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samples  of  possible  damage  patterns  and  to  assess 
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PREFACE 


This  report  describes  the  Airbase  Damage  Assessment  (AIDA)  computer 
model,  designed  for  rapid  examination  of  the  results  of  conventional  - 
air  attacks  on  complex  targets.  A user's  guide  and  program  listing  are 
included. 

AIDA  was  developed  by  Rand  for  use  in  a study,  conducted  at  Rq 
USAFE,  of  conventional  air  attack  effectiveness.  The  model  incorporates 
both  a Monte  Carlo  mode  and  a deterministic  (or  expected-value)  mode  of 
operation.  Approximations  used  in  the  expected-value  computations, 
while  differing  from  those  used  by  the  Joint  Munitions  Effectiveness 
Manual  (JMEM) , are  equivalent  in  precision  to  those  of  the  JMEM  hand 
methods.  A large  body  of  publications  furnishes  background  for  the 
methods  used  and  for  input  values  required  by  AIDA. 

AIDA  has  a variety  of  possible  applications.  It  can  be  used  as  an 
aid  in  planning  effective  attacks  against  complex  target  sets,  or — by 
testing  alternative  attack  headings,  aim  points,  and  the  like — in  examin- 
ing the  tradeoffs  between  damage  to  primary  and  secondary  targets.  It 
should  also  prove  useful  in  target  (e.g.,  airbase)  protection  studies, 
both  to  provide  realistic  and  detailed  samples  of  possible  damage  pat- 
terns and  to  assess  different  protection  options. 

AIDA  has  been  discussed  with,  and  made  available  to,  the  Operations, 
Intelligence,  Plans,  Logistics,  Engineering,  and  Communications  staffs 
at  Hq  USAFE;  to  the  staff  of  the  Assistant  Chief  of  Staff  for  Studies 
and  Analysis  at  Hq  USAF;  and  to  a number  of  other  DoD  and  NATO  organiza- 
tions. This  report  is  being  published  to  provide  a record  of  the  model 
and  to  make  it  available  to  a wider  audience.  The  computer  program  is 
available  from  The  Rand  Corporation. 

This  work  was  conducted  under  the  Project  RAND  research  project 
entitled  "Rand  Ramstein  Activity." 
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SUMMARY 
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This  report  describes  a new  computer  model  that  permits  examina- 
tion of  conventional  bombing  attacks  on  complex  targets — e.g.,  on  an 
airfield.  A complete  user's  guide  is  included.  Damage  statistics  for 
up  to  250  individual  targets  can  be  quickly  assessed  for  attacks  involv- 
ing as  many  as  50  delivery  passes  and  10  types  of  weapons.  Both  point- 
impact  weapons  (such  as  general-purpose  (GP)  bombs  and  precision-guided 
munitions  (PGMs))  and  area  weapons — cluster  bomb  units  (CBUs) — can  be 
handled,  and  targets  may  be  grouped  into  20  different  vulnerability 
categories  to  distinguish  different  levels  of  weapon  effectiveness. 

If  the  user  is  concerned  only  with  the  expected  numbers  of  hits  with 
point-impact  weapons,  and  is  not  interested  in  either  CBU  weapons  or 
the  coverage  and  damage  variations  expected  with  point-impact  weapons, 
a special,  more  efficient  expected-va.'  .e  mode  is  provided. 

In  its  basic  mode,  AIDA  determines  the  actual  impact  points  (pat- 
tern centroids  for  CBUs)  by  Monte  Carlo  procedures — i.e.,  by  random 
selections  from  the  appropriate  error  distributions.  GP  bombs  and  PGMs 
that  impact  within  a specified  distance  of  a target  are  classed  as  hits, 
and  the  results  include  the  total  number  of  hits  on  each  target  and  the 
cumulative  probability  of  kill.  For  CBU  munitions  the  program  assesses 
the  fraction  of  each  target  covered  by  each  pattern,  and  the  results 
include  the  fractional  coverage  from  all  patterns  and  a cumulative 
probability  of  kill  for  each  target.  In  addition  to  these  results  for 
the  complete  attack,  the  attack  can  be  repeated  automatically  for  sev- 
eral trials  to  provide  statistics  on  the  average  damage  levels  to  be 
expected.  In  the  special  expected-value  mode,  average  hit  densities 
are  determined  directly,  without  recourse  to  Monte  Carlo  procedures. 

In  the  basic  mode,  up  to  5 targets  may  be  designated  as  runways 
or  taxiways  suitable  for  aircraft  operations,  and  the  model  will  examine 
these  to  see  if  an  area  of  specified  size  is  available  for  such  opera- 
tions; if  not,  the  minimum  number  of  craters  that  would  need  to  be  re- 
paired to  obtain  an  area  of  that  size  is  determined. 

The  AIDA  program  is  designed  so  that  many  cases  can  be  examined 
conveniently  during  each  run — a feature  that  will  prove  useful,  for 
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example,  when  one  desires  to  quickly  examine  several  attack  options - 
Computer  output  includes  both  the  input  data  and  the  attack  results 
that  have  been  specified.  Most  model  features  are  illustrated  with 
a sample  problem. 
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I.  MODEL  DESCRIPTION 
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The  AIDA  (Airbase  Damage  Assessment-)  computer  model  permits  exami- 
nation of  bombing  attacks  on  a complex  set  of  targets — e.g.,  on  an  air- 
field. In  AIDA's  basic  mode  of  operation,  the  actual  bomb  impact  points 

* 

are  obtained  by  Monte  Carlo  procedures,  and  the  attack  can  be  repeated 
for  several  trials  to  provide  statistical  estimates  of  the  average  dam- 
age and  variability  of  that  damage  for  each  of  the  many  targets.  Al- 
ternatively, an  expected-value  mode  is  offered  when  only  point -impact 
weapons  are  employed  in  the  attack  and  when  the  user  is  interested  only 
in  the  expected  numbers  of  hits.  With  a variant  of  the  expected-value 
mode,  one  may  also  generate  hit-density  patterns  for  complex  attacks 
without  specifying  an  actual  target  system.  Both  modes  of  operation 
and  several  different  sets  of  problems  may  be  treated  by  successive 
cases  during  a single  computer  run.  The  several  features  available 
with  AIDA  are  illustrated  in  Sec.  IV. 

-j- 

In  AIDA  the  target  system  may  be  composed  of  up  to  250  separate 
targets,  for  example,  shelters,  hangers,  maintenance  buildings,  runways, 
taxiways — even  pipelines.  The  complete  attack  may  consist  of  up  to  50 
distinct  weapon-delivery  passes.  Each  target  is  a rectangle  of  specific 
size  and  orientation  and  an  attack  pass  is  defined  by  the  expected  proba- 
bility of  arrival,  a heading,  and  the  aim  point,  delivery  accuracy,  and 
dispersion  for  a stick  of  weapons.  Targets  may  be  grouped  into  10  or 
20  different  vulnerability  categories  and  there  may  be  up  to  10  differ- 
ent kinds  of  weapons  dropped  in  an  attack. 

This  computer  model  has  a variety  of  possible  applications.  Since 
an  entire  airbase  may  be  represented  rather  accurately,  AIDA  can  be  used 
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That  is,  the  actual  mean  point  of  impact  and  the  actual  impact 
points  are  determined  by  random  variates  drawn  from  the  appropriate 
error  distributons. 
t 

Most  array  dimensions  may  be  changed  fairly  readily,  and  they  have 
been  in  several  instances.  Some  of  the  dimensions  mentioned  in  this 
report  differ  from  those  outlined  in  an  earlier  version  of  AIDA;  the 
changes  were  made  to  maximize  program  utility  while  limiting  core  stor- 
age requirements.  Users  will  find  instructions  for  changing  program 
dimensions  at  the  beginning  of  AIDA’s  MAIN  subroutine. 
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to  study  how  one  might  best  design  an  attack  to  damage  both  primary  and 

secondary  targets.  Whereas  many  bomb-damage  rids  deal  only  with  a single 

target  (or  set  of  like  targets)  and  provide  no  direct  evidence  as  to  the 

collateral  damage  to  be  expected,  AIDA  could  be  used  to  help  design  an 

attack  that  would  get  the  maximum  benefit  from  those  bombs  that  miss  the 

primary  target.  AIDA  should  also  prove  usef  ' in  a variety  of  studies 

of  base  protection  and  repair,  in  that  it  could  pro\ide  a more  realistic 

and  detailed  picture  of  possible  damage  than  is  customarily  available. 

Yet  another  area  in  which  AIDA  may  prove  useful  is  as  the  first  step  in 

an  analysis  of  the  sortie-generation  capabilities  of  a damaged  base. 

Hopefully  the  results  of  this  model  can  be  used  as  input  co  a sortie- 
* 

generation  model,  and  the  combination  may  then  be  able  to  provide  a 
much  better  understanding  of  the  effects  of  attacks  on  aircraft  opera- 
tions than  exists  today.  In  addition  to  analyzing  problems  involving 
airfield  complexes,  AIDA  could  just  as  easily  be  used  to  assess  attack 
options  and  expected  damage  levels  for  attacks  against  the  complex  tar- 
get arrays  found  at  field  headquarters,  SAM  sites,  supply  depots,  etc. 

Although  AIDA  is  straightforward  and  simple  mathematically,  its 
substantial  flexibility  should  prove  useful  in  a wide  range  of  studies. 
For  the  reader  who  wants  more  detail  on  the  methodological  treatment 
than  is  discussed  in  the  body  of  this  report,  a program  listing  is  in- 
cluded in  Appendix  C. 

MONTE  CARLO  MODE 

In  this  mode  of  AIDA,  weapons  may  be  of  two  basic  types:  point- 

impact  weapons  (such  as  general-purpose  (GP)  bombs  and  precision-guided 
munitions  (PGMs)  or  area  weapons  (such  as  cluster  bomb  units  (CBUs)).  A 
weapon  reliability  may  be  associated  with  each  kind  of  weapon.  For 
each  kind  of  point-impact  weapon  an  effective  miss  distance  (EMD)  may 
be  specified  for  each  target  type  (i.e.,  that  miss  distance  at  which 
the  weapon  is  effective  and  an  impact  is  to  be  categorized  as  a hit). 


One  that  can  account  for  the  effects  of  various  levels  of  damage 
to  the  various  elements  of  a base. 
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Approximations  used  in  the  AIDA  computations,  while  different 
from  those  used  in  the  Joint  Munitions  Effectiveness  Manual  (JMEM),  are 
about  equivalent  in  precision  to  those  of  the  JMEM  hand  methods. 
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When  this  is  done,  target  coverage  is  computed  as  that  fraction  of  the 
target  area  that  is  covered  by  a circle  having  a radius  of  EMD  and 
centered  at  the  impact  point.  Additionally,  the  user  may  specify  a 
value  for  target  kill  probability,  given  a hit  (as  defined  above),  or 
he  may  specify  a different  radius  than  the  EMD  and  that  will  also  be 
used  for  computing  a value  of  coverage.  This  option  may  be  specified 
individually  for  each  target  type  and  weapon  type,  as  with  the  EMD. 

For  example,  with  this  feature  one  can  use  the  EMD  to  represent  the 
radius  of  the  mean  area  of  effectiveness  (MAE)  for  severe  structural 
damage  to  a building,  and  the  alternate  value  to  represent  the  radius 
of  the  MAE  for  severe  damage  to  the  building's  contents,  given  a hit 


on  the  structure.  The  final  results  will  include  the  cumulative  cover- 


age fractions  for  each  target  for  all  point-impact  weapons,  correspond- 
ing to  both  the  EMD  and  to  the  optional  factor,  computed  according  to 


the  m l” tier: 


FC  = 1 - 


(1  - C (i»  , 

Wjl 


all  hits 


where  in  the  first  case,  Cw  ^(i)  is  the  target  coverage  based  on  the 


EMD  of  the  i th  hit  on  the  target  of  type  t by  a weapon  of  type  w;  in 
the  second,  it  is  assessed  according  to  the  user  specified  option. 

If  weapons  are  specified  as  CBUs,  the  model  first  computes  the 
fraction  of  the  area  of  each  target  that  is  covered  by  the  rectangular 
bomblet  pattern.  The  total  fractional  coverage  of  a target  for  all 
passes  is  that  fraction  of  the  area  that  has  been  covered  by  one  or 
more  patterns.  If  a probability  of  kill  is  associated  with  "coverage" 
by  CBU  patterns,  the  model  will  generate  the  total  probability  of  kill, 
taking  into  account  the  actual  position  of  each  of  the  CBU  patterns 
that  covered  any  portion  of  the  target,  according  to  the  relation: 


PK  = 1 


SlUK, 


N(a,i)(  / 
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where 


a = the  attack  number, 
i = a point  on  the  target  grid, 

N(a,i)  = the  number  of  times  point  i was  "covered"  during 
attack  a, 

M - the  total  number  of  CBU  attacks, 

T = the  total  number  of  target  grid  positions,  and 


= the  probability  of  kill  of  a portion  of  a target 


w,t 


of  type  t that  is  "covered"  by  a bomblet  pattern 


from  a weapon  of  type  w. 


The  results  for  each  trial  include  the  number  of  hits  by  point- 
impact  weapons  and  the  fractional  coverage  by  CBUs  for  each  target  as 
well  as  the  point-impact  weapon  coverage  (FC)  and  CBU  kill  probability 
(PK) ; in  addition,  for  the  targets  that  the  user  has  specified  (for  a 
maximum  of  20  targets  other  than  the  runways  and  taxiways),  the  impact 
points  and  weapon  types  are  printed  for  up  to  25  weapons  per  target. 

For  multiple  trials,  the  results  for  each  target  include  the  fraction 
of  trials  with  at  least  one  hit,  the  average  number  of  hits  and  aver- 
age CBU  coverage,  the  standard  deviation  of  these  two  measures,  and 
the  average  values  of  FC  and  PK  for  the  several  trials.  A full  descrip- 
tion of  the  output  options  is  presented  in  the  next  section  and  is  sum- 
marized in  Table  1,  p.  9. 

£ 

The  user  may  also  specify  that  certain  (up  to  5)  of  the  (rectangu- 
lar) targets  are  actually  runways  or  taxiways  that  are  suitable  for  air- 
craft operations.  The  model  will  then  test  to  see  if  such  operations 
are  possible  from  these  areas;  i.e.,  tests  are  made  to  see  if  the  min- 
imum clear  length  and  width  required  for  operations  are  available  after 
the  attack.  In  testing  for  runway  availability,  only  point-impact  weap- 
ons are  considered,  and  the  crater  radius  is  taken  as  the  EMD.  Up  to 
250  hits  can  be  stored  and  examined  for  each  such  target  (the  locations 
will  be  listed  along  with  hi t data  on  other  targets) . If  the  runway 
does  not  meet  the  minimum  requirements,  the  user  may  request  an  assess- 
ment of  how  few  craters  would  need  to  be  repaired  to  meet  those  require- 
ments. The  user  may  also  request  an  approximate  computer  plot  of  the 
impact  points  for  each  runway. 


By  d' "ining  as  type  #1  targets. 
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Coverage  by  CBU  patterns  is  not  taken  into  account  in  the  runway 
* 

assessments;  examination  of  the  effects  of  CBUs  on  the  runway  would 
require  that  the  user  either  make  a visual  estimate  using  the  approxi- 
mate computer  plot  or  plot  the  impact  points  and  the  resulting  bomblet 
patterns  more  accurately  on  a plan  of  the  airbase. 


EXPECTED-VALUE  MODE 

This  mode  operates  with  target  and  attack  descriptions  that  are 
identical  to  those  used  with  the  Monte  Carlo  mode  of  operation.  How- 
ever, only  point-impact  weapons  may  be  included  in  the  attack.  Com- 

f 

putationally  this  mode  derives  an  average  value  of  the  hit  density 
for  each  target  and  for  each  attack.  These  are  combined  to  provide 
the  total  expected  number  of  hits  for  all  the  attacks.  Althougn  this 
mode  provides  no  evidence  regarding  the  statistical  variations  that 
must  he  expected  in  actual  hit  patterns,  it  does  p~ovide  a quick,  ef- 
ficient means  of  assessing  the  expected  values.  Since,  for  many  prob- 
lems, a rather  large  number  of  Monte  Carlo  trials  is  required  in  order 
to  get  a reasonably  accurate  estimate  of  average  values,  one  can  some- 
times use  both  the  Monte  Carlo  and  expected-value  modes  to  advantage. 

The  Monte  Carlo  mode  can  first  be  used  with  a limited  number  of  trials 
to  provide  gross  estimates  cf  target  coverage  and  of  the  variability 
of  hits;  the  expected-value  mode  can  then  speedily  provide  a reasonably 
accurate  estimate  of  the  average  number  of  hits.  As  will  be  illustrated 


m 


Except  that  the  centroids  are  shown  on  the  computer  plots  when 
they  fall  within  the  area  included  in  the  plot. 
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The  point  value  of  the  hit  density  is  first  determined  at  several 
points  on  each  target  for  the  specified  aiming  errors  and  dispersion; 
these  values  are  computed  for  each  of  several  positions  along  the  in- 
tended line  of  bomb  impacts.  The  average  hit  density  is  taken  as  the 
arithmetic  average  of  these  point  values  for  the  several  target  points 
and  the  several  positions  along  the  bomb  impact  line.  If  the  target 
dimensions  are  small  (i.e.,  less  than  one-quarter  the  RMS  of  the  range 
and  deflection  errors  projected  parallel  to  each  target  edge),  only  the 
target  corners  are  used  for  averaging.  For  larger  targets,  a uniform 
grid  of  internal  points  is  established  with  spacing  no  greater  than 
one-quarter  the  RMS  noted  above.  The  expected  number  of  hits  is  com- 
puted as  the  product  of  the  average  hit  density  and  the  area  of  the 
target  (including  a border  as  wide  as  the  EMD) . 
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in  Sec.  IV,  this  two-step  calculation  would  require  only  two  more  cards 
than  would  either  calculation  alone. 

A special  feature  of  the  expected-value  mode  permits  the  user  to 
quickly  generate  a hit-density  grid.  This  feature  is  controlled  by 
special  target  cards  that  may  be  used  either  alone,  or  in  conjunction 
with  normal  target  cards.  Each  such  card  generates  a 17  x 17  grid  of 
hit-density  values  measured  over  a square  of  specified  dimensions  (the 
dimension  should  be  a multiple  of  16) . If  no  dimension  is  specified, 
hit  densities  are  provided  at  250  ft  intervals  over  a 4000  x 4000  foot 
square.  The  southwest  corner  for  each  grid  is  placed  at  that  position 
specified  on  the  special  target  card.  Such  cards  are  identified  by 
specifying  target  type  #21 — an  entry  that  acts  as  a control  signal  with- 
in AIDA. 

AIDA  OPERATION 

AIDA  offers  several  features  designed  to  simplify  its  operation  and 
to  permit  a series  of  cases  to  be  analyzed  during  a single  computer  run. 
Most  are  illustrated  in  Sec.  IV.  The  first  feature  permits  a multi- 
aircraft attacK  against  the  same  objective  to  be  specified  simply;  when 
two  or  more  attacks  have  common  parameters  (i.e.,  heading,  desired  mean 
point  of  impact  (DMPI),  CEP,  dispersion,  arrival  probability),  a single 
entry  will  generate  the  additional  attacks.  Other  convenience  features 
derive  from  use  of  the  REDO  card  (see  the  next  section).  When  this  card 
is  encountered  it  acts  as  a terminator  card,  terminating  the  input  for 
one  case  and  announcing  that  there  will  be  a subsequent  case.  If  de- 
sired, an  entirely  new  set  of  data  may  be  input  following  a REDO  card. 
Alternatively,  the  subsequent  case  may  be  a modification  of  the  preced- 
ing case;  a few  simple  inputs  permit  the  user  a wide  range  of  modifica- 
tion options.  One  may  change  AIDA's  mode  of  operation  or  any  of  the 
control  variables.  In  addition,  targets  may  be  added  to  the  prior  list; 
or  some  may  bo  removed  and  then  others  added.  The  same  can  be  done 
with  attacks.  For  example,  this  feature  could  prove  particularly  use- 
ful in  investigating  a series  of  alternative  attack  options  on  a given 
airfield. 
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AIDA's  features  provide  substantial  flexibility.  For  example, 
since  target  location  is  not  restricted,  it  is  possible  to  have  two 
identical  targets  at  the  same  location;  by  assigning  these  targets 
different  type  numbers,  one  can  assess  results  for  two  weapon  effects. 
This  would  be  useful,  for  example,  in  determining  expected  personnel 
and  materiel  losses  (targets  with  different  MAEs)  in  open  parking 
areas.  The  kill  probability  options  available  with  impact  weapons 
complement  this  technique.  Consider,  for  example,  two  identically 
located  buildings  (of  different  target  types),  each  with  EMDs  fixed 
by  the  MAE  for  a structural  kill.  If  one  also  inputs  an  effective  kill 
radius  based  on  equipment  damage  for  one  of  the  buildings  and  the  kill 
radius  for  personnel  with  the  other,  the  Monte  Carlo  results  will  in- 
clude estimates  of  the  structural  damage,  as  well  as  of  the  materiel 
damage  and  personnel  losses  resulting  from  a structural  hit. 
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II.  INPUT  PROCEDURES 


Up  to  seven  basic  types  of  cards  can  be  used  in  operating  AI1)A, 
although  only  three  are  required.  Four  types  describe  the  target  and 
. attack  characteristics  and  three  others  are  used  to  control  operations. 
For  each  basic  card  the  type  i9  identified  in  the  first  four  columns 
(left-adjusted): 

TGT  target  data  card 

ATT  attacker  data  card 

ATT2  alternate  attacker  data  card  (optional) 

EMD  weapon  card  (optional) 

CONT  control  card  (optional) 

REDO  controls  sequential  cases  (optional) 

END  last  card 

For  two  of  these  basic  card  types — ATT2  and  EHD — supplementary  cards 

are  used.  For  the  ATT2  card  a following  card  (with  additional  data, 

as  will  be  described)  is  mandatory.  One  to  three  suppleme.i£ary  cards 

may  be  employed  with  the  EMD  card,  as  will  be  described. 

There  may  be  as  many  as  250  TGT  cards,  50  ATT  or  ATT2  cards,  and 

10  EMD  cards.  For  a given  case  there  may  be  at  most  1 CONT  card.  The 

* 

order  of  the  cards  is  immaterial  except  that  a REDO  card  or  an  END 
card  must  be  used  to  signify  the  completion  of  input  for  a given  case. 
The  targets  and  attackers  are  numbered,  internally,  in  the  sequence 
in  which  their  descriptions  are  read  in;  each  target  may  also  have  an 
alphanumeric  designator  ('■>.. g. , building  number).  A detailed  descrip- 
tion of  how  data  are  to  be  entered  on  each  type  of  card  is  presented 
in  Appendix  A.  The  input  data  are  printed  for  each  case;  Table  1 out- 
lines the  output  options  for  the  results.  The  program  will  function 
in  its  simplest  form  with  only  TGT  cards,  ATI'  or  ATT2  cards,  and  an 
END  card. 

£ 

The  only  exception  is  that  if  the  special  ’’target  type  921"  cards 
are  used  in  conjunction  with  the  special  hit-density  grid  feature,  one 
of  these  cards  must  be  the  last  target  card. 
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Table  1 

OUTPUT  CONTROL 


NPRINT° 

Control 

Value 


For 

Each  Trial 

All 

All  Hits 

Impact 

(and  Target 

Stored 

Hit 

Runway 

Points 

Corners) 

Hit  Data 

Summary 

Results 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

x 

X 

X 

X 

x 

xb 

c 

X 

Multiple 

Trial 


*Enter  in  Columns  23,  24  of  CONT  card. 

^Compact  listing  of  hits  and  required  repairs  for  runways/taxiways . 
'Compact  listing  of  hits  on  each  target. 


A control  (CONT)  card  will  be  required  if  advantage  is  to  be  taken 
of  more  than  the  most  basic  of  AIDA’s  features;  without  this  card  AIDA 
examines  only  one  Monte  Carlo  sample  of  the  attack  and  provides  the 
actual  numbers  of  hits  on  all  targets  and  the  stored  hit  locations  for 
specified  targets.  Specifically,  a CONT  card  is  needed  if  (1)  more  than 
one  trial  is  required,  (2)  an  alternative  output  mode  is  desired,  (3)  a 
different  mode  of  operation  is  desired  (e.g.,  the  expected-value  mode), 
or  (4)  the  runway  availability  features  are  to  be  exercised.  As  ex- 
plained more  fully  on  pp.  38-39,  this  card  is  used  to  specify  the  number 
of  trials,  the  mode  of  operations,  the  output  formats,  the  dimensions  of 
the  minimum  runway  surface,  whether  or  not  the  user  wants  the  minimum 
repair  requirement  assessed,  and  the  distances  that  the  "minimum  runway 
rectangle"  will  be  shifted  laterally  and  longitudinal 1y  in  testing  for 
an  available  area. 
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TARGETS 

• Reference  (westernmost)  comer 
Ll  Northeasterly  heading  boundary 
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9 Orientation  angle 


X-axis 

ATTACK  PASSES 
* DMPI 

f>  Attack  heading 
SL  Stick  length 
f — I"  Nominal  bomb  impacts 


Fig.  1 — Target  and  attack  layout  for  AIDA 
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Several  special  features  should  be  noted  with  respect  to  the  EMD 
cards.  To  begin  with,  these  cards  are  optional;  if  the  EMD  card  is 
not  used  for  one  or  more  weapon  types,  these  weapons  are  assumed  to  be 
point-impact  weapons  that  must  impact  within  the  target  outline  to  be 
counted  as  "hits."  If  more  of  the  features  are  to  be  used  for  a par- 
ticular type  of  weapon,  an  EMD  card  is  input  with  the  number  of  the 
weapon  entered  in  Columns  11  and  12.  If  weapon  reliability  is  less 
than  unity,  a 1 is  placed  in  Column  6 of  the  EMD  card  and  the  relia- 
bility is  entered  in  Columns  7 to  12  of  a special  card  placed  immed- 
iately following  the  EMD  card.  If  the  user  wants  to  consider  20  types 
of  targets,  rather  than  only  10,  a 1 is  placed  in  Column  5,  and  data 
for  the  second  set  of  10  target  types  are  entered  on  a supplementary 
card  (2  cards  if  a 1 also  appears  in  Column  6) . The  other  entries  for 
the  EMD  card  and  the  special  following  cards  differ,  depending  on 
whether  it  is  a CBU-type  munition  or  a point-impact  weapon. 

CBU  munitions  are  denoted  by  a negative  number  in  Columns  13  to 
18  on  the  EMD  card.  The  absolute  value  of  that  number  is  taken  as  the 
CBU  pattern  dimension  along  the  flight  direction;  pattern  width  is 
specified  in  Columns  19  to  24.  If  tne  user  wants  to  associate  a prob- 
ability of  kill  with  CBU  coverage  on  some  or  all  of  the  target  types, 
that  fact  is  also  denoted  by  the  integer  1 in  Column  6 of  the  EMD  card. 
When  this  is  done  the  program  interprets  entries  on  the  following  card 

(also  used  for  weapon  reliability)  as  p,s  for  that  weapon  against  the 

* 

desired  target  types. 

Other  special  features  are  available  for  use  with  point- impact 
weapons.  If  a weapon  can  effectively  damage  a target  when  it  actually 
falls  near  but  outside  the  target  outline,  the  EMD  for  a "hit"  can  be 
entered  in  Columns  13  through  72  (in  10  fields  of  6 columns)  of  the  EMD 
card  (and  auxiliary  card)  for  the  10  (or  20)  target  types.  The  appro- 
priate entry  in  many  instances  would  be  the  radius  of  a circle  whose 
area  is  equal  to  the  MAE  (as  presented  in  JMEM  manuals)  for  the  corre- 
sponding target-weapon  combination.  In  the  case  of  hits  on  runways  or 
taxiways,  however,  the  appropriate  entry  is  crater  radius;  when  AIDA 

* 

In  ten  6-column  fields  from  Columns  13  to  72. 
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checks  for  the  availability  of  a minimum  runway,  each  reliable  impact 
is  assumed  to  have  a crater  radius  equal  to  EMD. 

For  each  "hit"  with  a point-impact  weapon  an  estimate  is  made  of 

the  fraction  of  the  target  area  that  is  covered  by  a circle  of  radius 
* 

EMD.  The  user  may  also  select  one  of  the  additional  options  (see 

Fig.  2) ; using  the  supplementary  card  he  either  may  input  a p^  for  a 

particular  weapon-target  type  or  specify  a radius  different  than  EMD 

for  defining  another  "effects"  circle  to  be  used  in  computing  an  addi- 

•f* 

tional  coverage  measure.  In  either  case  (p^  or  radius)  the  value  is 
entered  on  the  supplementary  card  in  the  field  corresponding  to  the 
appropriate  target  type. 

A REDO  card  terminates  the  input  for  one  case  and  initiates  a new 
case.  If  no  additional  entries  are  made  on  the  REDO  card,  the  targets 
and  attacks  for  the  next  case  will  include  all  those  entered  in  the 
preceding  case.  If  some,  but  not  all,  of  the  preceding  targets  and/or 
attacks  are  to  be  included,  the  number  of  targets  that  are  to  be  re- 
tained is  entered  in  Columns  7 to  12  and  the  number  of  attacks  to  be 
retained  is  entered  in  Columns  13  to  18.  The  numbers  to  be  retained 
are  selected  from  the  beginning  of  the  ordered  lists  generated  in  the 
input  process.  If  a negative  entry  (e.g.,  -1)  is  made  in  either  or 
both  of  the  target  or  attack  fields,  none  of  the  prior  targets  and/or 
attacks  will  be  considered  in  the  new  case. 

After  specifying  which — if  any — of  the  prior  inputs  are  to  be  re- 
tained, additional  targets  and  attacks  may  be  added  with  appropriate 
TGT  and  ATT  cards.  If  a 1 is  entered  in  Column  24  of  the  REDO  card, 
the  target  list  and/or  the  attack  and  weapons  lists  will  not  be  in- 
cluded in  the  input  listing,  if  no  change  has  been  made  in  those  lists 
from  the  prior  use. 

Control  conditions  will  remain  as  in  the  preceding  case  unless  a 
new  CONT  card  is  entered  with  appropriate  data.  An  END  card  must  be 
used,  at  the  end  of  the  last  case,  to  complete  the  card  deck. 


As  noted  earlier,  the  cumulative  coverage  of  a target  is  given 
by  1 minus  the  product  of  the  probability  of  noncoverage  to  all  hits; 
this  is  equivalent  to  saying  that  in  estimating  the  target  area  cov- 
ered by  a particular  weapon  the  effect  of  prior  hits  is  neglected. 

+ 

A value  greater  than  1 is  interpreted  as  a radius. 


i>  w*!"^  * ^ v^V^ff5 


X^KiC- r * - jTfg ^'?^^.;-<^>»g^»^i»iX^^gJiiBtfaaMt^  'jjtifi\L 


III.  PROGRAM  NOTES 
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AIDA  is  written  in  FORTRAN  IV  and  first  operated  on  the  Honeywell 

6060  used  with  Hq  USAFE's  Worldwide  Military  Command  and  Control  System 

(WWMCCS);  subsequently,  it  has  been  run  on  the  IBM  360/370  series.  AIDA 

is  easily  transferable  to  other  computers  because  there  are  no  system- 
ic 

unique  features  employed.  The  program  is  composed  of  approximately 
1950  card  images  organized  into  a MAIN  routine  and  17  subroutines.  As 
presently  dimensioned,  AIDA  operation  requires  36K  words  of  core.  The 
organization  of  the  various  subroutines  and  their  functions  are  indi- 
cated in  Fig.  3;  this  structural  outline  of  AIDA  should  also  prove  use- 
ful for  anyone  who  wants  to  overlay  the  program  so  as  to  reduce  core 
storage  requirements.  Potential  users  will  find  that  the  source  code 
is  moderately  well  annotated  with  comment  cards.  A full  listing  of 
the  definitions  for  the  key  variables  and  arrays  xtfill  be  found  in 
Appendix  B. 

AIDA's  CPU  requirements  obviously  can  vary  widely  depending  on 
the  complexity  of  the  problem  (i.e.,  numbers  of  targets  and  attacks 
and  use  of  the  special  runway  features)  and  upon  the  required  number  of 
replications.  For  problems  involving  about  150  targets,  750  bombs,  and 
20  replications  of  the  attack,  computing  times  have  run  between  1-1/2 
to  2 minutes  on  the  Honeywell  machine.  To  limit  these  requirements  the 
process  used  to  test  for  hits  with  point-impact  weapons  has  been  de- 
signed to  reduce  the  TARGETS  X BOMBS  dimensionality  problem  somewhat, 
by  ordering  targets  in  a particular  manner  and  by  only  checking  those 
targets  "near"  each  bomb  impact  point.  Nevertheless  there  is  a very 
substantial  amount  of  checking  and  computation  required. 

The  computation  problem  is  especially  severe  when  CBU  munitions 
are  used.  To  estimate  the  joint  coverage  of  all  CBU  patterns,  the  tar- 
gets each  have  a uniform  grid  of  points  superimposed,  and  each  of  those 
points  is  checked  to  see  if  it  lies  within  one  or  more  of  the  patterns. 


Minor  format  changes  may  be  required  in  a few  instances;  name 
changes  will  probably  be  necessary  for  the  random  number  generator. 
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MAIN  controls  program  operations 
and  initializes  numerous  arrays 

Reads,  organizes,  and  lists 
input  data 

Performs  bomb  trajectory 
calculations  as  required 

Computes  and  stores  target 
corner  locations 

Computes  and  prints  hit- 
density  grid  on  demand 

Computes  and  prints  expected 
number  of  hits  for  point-im- 
pact weapons 

Orders  targets  by  location 
of  westernmost  corner 
Groups  targets  into  zones 

Determines  actual  MP1  and 
weapon-impact  points 


INPUT 


JMEMO 


TGTDIM 


TGTORD 

tgtzon’ 


GAUSS 


MAIN 


GRIDEN 


EXPHIT 
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t 

1 if 

■X  Jag 
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Checks  for  hits  with  point 
impact  weapons 

Checks  for  coverage  by 
CBU-type  weapons 

Controls  runway  tests 

Checks  for  runway  mini  mums 
anc!  determines  minimum  repairs 
Plots  impact  points  on/near 
runways  and  main  taxi ways 

Prints  results  from  individual 
trials 


_TESTHT_ 

HITTGT 


CHECKR 


RUNWAY 


PLOTHT 


PRINT 


More 

trials? 


Computes  and  prints  the 
statistical  results  for 
multiple  trials 


STATIS 


More  cases? 


Fig. 3 — Organization  and  function  of  the  AIDA  ^routines 
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The  fraction  of  the  target  grid  points  that  is  covered  by  one  or  more 

CBU  patterns  is  used  as  the  estimate  of  fractional  coverage.  For  small 

targets  only  a 4 * 4 grid  is  used,  but  for  the  very  large  targets  (di- 

* 

mensions  in  excess  of  1000  x 1000)  a 16  x 16  grid  is  used;  quite 
obviously  the  computational  requirements  are  affected  by  the  presence 
of  large  targets  and  the  use  of  CBUs,  but  not,  fortunately,  in  direct 
proportion  to  the  number  of  grid  points. 

The  procedure  used  for  checking  the  availability  of  an  adequate 
section  of  runway  also  involves  substantial  processing.  A rectangle 
of  the  required  dimensions  is  positioned  first  at  a corner  of  the  run- 
way and  tested  to  see  if  there  is  a hit  within  the  rectangle.  If  there 

•j* 

is,  the  rectangle  is  moved  5 ft  laterally  and  rechecked;  this  is  re- 
peated across  the  runway  or  until  an  open  area  has  been  found.  If  none 

i 

is  found  the  rectangle  is  shifted  250  ft  along  the  runway  and  the  pro- 
cess is  repeated.  If  the  dimensions  of  the  minimum  usable  runway  sec- 
tion are  small,  or  if  the  runway  is  large  and/or  there  are  many  hits,  a 
very  substantial  amount  of  processing  can  be  required.  These  require- 
ments are  further  increased  when  an  assessment  of  the  minimum  repair 
requirements  is  requested.  Furthermore,  if  there  are  other  runways, 
or  taxiways  suitable  for  emergency  flight  operations  (signified  by  entry 
as  a target  type  #1) , the  entire  process  must  be  repeated  for  each  sur- 
face (for  a maximum  of  5) . Whenever  the  problem  can  be  effectively 
reduced  to  one  of  establishing  the  availability  of  the  minimum  clear 
width  (e.g.,  when  a 9C00  ft  runway  is  being  cut  at  the  center  to  deny 
a 6000  ft  clear  length) , the  minimum  clear  length  should  be  set  equal 
to  the  runway  length  to  avoid  going  through  the  lateral  shift  sequence 
unnecessarily  (11  times  in  the  example). 


It  will  be  noted  that  hits  may  not  be  recorded  when  a target  is 
very  large  and  the  CBU  pattern  so  small,  proportionately,  that  the 
pattern  can  fall  between  grid  points. 

+ 

These  are  the  default  values;  the  user  may  specify  a different 


value . 
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IV.  SAMPLE  PROBLEM 


The  layout  of  Base  XYZ  is  shown  in  Fig.  4.  This  base  consists  of 
a 150  ft  x 8000  ft  main  runway,  several  taxiways,  a parking  ramp,  eight 
support  facilities,  and  a housing  area.  To  examine  the  effectiveness 
of  bombing  attacks  against  this  target  complex  with  AIDA,  one  first 
must  describe  the  target  elements  and  the  attack  in  a common  coordinate 
system.  Targets  are  defined  by  their  westernmost  corner,  their  size, 
and  their  orientation;  the  attack  heading  and  the  desired  mean  point 
of  impact  fix  the  attacks.  For  this  illustration  four  medium  bombers 
will  each  drop  25  bombs  in  an  effort  to  cut  the  runway  at  two  points; 
two  others  are  targeted  on  the  operations  building  near  the  main  taxi- 
way, and  one  will  aim  at  the  electronics  shop.  In  addition,  one 
fighter-bomber  is  assigned  to  dive-bomb  each  of  the  main  aircraft  main- 
tenance buildings,  B1  and  B2,  and  one  will  drop  a stick  of  5 CBUs  on 
the  housing  area. 

INPUT 

Figures  5a  and  5b  reproduce  the  card  images  needed  to  describe 
this  sample  problem  and  to  control  four  distinctly  different  assessments. 
These  assessments  are  iesignated  as  "cases"  and  illustrate  the  following: 

Case  1.  Hit  statistics  for  all  targets,  using  5 trials. 

Case  2.  Detailed  results  of  a single  trial. 

Case  3.  Results  using  the  expected-value  mode,  with  a sample 
hit-density  grid. 

Case  4.  Summary  results  of  the  runway  attack,  using  25  trials. 

For  clarity,  the  control  (CONT)  card,  the  target  (TGT)  cards,  the  weapon 
(EMD)  cards,  and  the  attack  (ATT  and  ATT2)  cards  are  listed  in  order; 
the  ordering  of  cards  for  a given  case  is  actually  immaterial  except 
that  the  EMD  card  pairs  and  the  ATT2  card  pairs  must  be  together.  A 
REDO  card  or  an  END  card  defines  the  end  of  the  input  data  for  a case. 
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A careful  review  of  the  various  annotations  in  Figs.  5a  and  5b 
along  with  a reading  of  Appendix  A,  hopefully  will  lead  to  a full  and 
rapid  understanding  of  the  various  input  requirements  and  optional 
features  provided  by  AIDA.  As  can  be  seen  in  Case  1,  the  first  card 
directs  that  the  Monte  Carlo  mode  be  used  and  that  a statistical  sum- 
mary of  five  replications  of  this  attack  be  printed;  no  assessment  of 
runway  availability  is  requested  for  this  case.  The  specifications 
for  the  targets  are  straightforward.  The  weapon  data  are  somewhat  more 
complex;  two  cards  are  used  in  describing  each  weapon’s  effectiveness 
against  the  10  types  of  target  (four  cards  would  be  required  for  11 
to  20  target  types).  The  first  card  lists  the  EMD  for  10  target  types 
and  the  second  card  lists  weapon  reliability  and  the  optional  effective- 
ness descriptor  (an  alternate  radius,  or  a p^) . In  this  example  weap- 
on type  //I  has  an  effective  (crater)  radius  against  paved  surfaces 
(target  types  #1,  #2,  and  #3)  of  22  ft;  against  target  types  #4  and  #5 
a hit  will  be  assessed  if  the  impact  is  within  40  ft  of  these  structures. 
The  fractional  damage  on  these  structures,  given  a hit,  is  to  be  based 
on  a 50  ft  effects  radius  for  target  type  #4,  and  a 75  ft  radius  for 
target  type  #5.  The  third  weapon  type  is  a CBU  and  the  special  input 
procedures  are  illustrated. 

The  10  attacks  are  listed  last.  The  first  (and  fourth  through 
tenth)  attack  uses  the  simpler  ATT  specification;  the  REP,  DEP,  disper- 
sion, and  stick  length  are  given  directly  in  feet,  measured  in  the 
ground  plane.  For  the  second  and  third  attacks — the  fighter-bombers 
directed  to  attack  the  main  maintenance  hangars  B1  and  B2 — these  fac- 
tors are  to  be  computed,  based  on  the  aircraft's  attack  characteristics; 
both  will  release  their  last  bomb  at  a 5000  ft  altitude  in  a 45  deg  dive 
at  450  kn  with  an  intervalometer  setting  of  100  ms.  Aiming  accuracy  and 
dispersion  in  the  plane  normal  to  the  trajectory  are  30  mils  and  5 mils, 
respectively.  For  three  of  the  attacks,  two  aircraft  will  attempt  to 
fly  the  identical  path,  so  only  one  card  is  needed  to  describe  each 
pair. 

The  remaining  eight  cards  are  all  that  are  needed  to  define  and 
control  three  additional  cases.  Case  2 simply  calls  for  a single  Monte 
Carlo  attack,  but  with  full  printout  and  with  an  examination  of  the 
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availability  of  an  undamaged  50  ft  x 4000  ft  section  of  either  the 
runway  or  main  taxiway  for  aircraft  operations. 

Case  3 requests  an  estimate  of  the  expected  numbers  of  hits  using 
the  expected-value  mode;  since  this  mode  will  not  function  for  CBU-type 
munitions,  only  the  first  nine  attacks  are  retained  (the  housing  area 
targets  for  the  tenth  attack  are  also  dropped).  The  target  card  that 
is  added  in  this  case  directs  the  additional  computation  of  a 1600  ft 
x 1600  ft  hit-density  grid,  located  as  noted. 

Case  4 focuses  on  the  availability  of  a minimum  surface  for  air- 
craft operations.  Only  the  runway  and  main  taxiway  are  retained  as 
targets;  however,  all  attacks  are  to  be  considered,  since  any  crater 
must  be  taken  into  account,  whether  it  represents  a hit  on  the  intended 
target,  or  collateral  damage  from  some  different  attack.  Twenty-five 
attack  trials  are  to  be  run  with  the  Monte  Carlo  mode,  repair  require- 
ments are  to  be  assessed,  and  the  crial-to-trial  runway  results  are  to 
be  printed  along  with  the  statistical  summary  for  the  25  trials. 


OUTPUT 

The  very  first  output  (Fig.  6),  presented  even  before  the  main 
summary  of  the  input  data,  lists  the  input/output  for  any  trajectory 
calculations  that  were  required;  in  this  case  this  was  necessary  for 
attacks  2 and  3.  The  basic  input  listing  follows;  as  will  be  noted, 
the  factors  determined  with  the  trajectory  calculations  are  listed  in 
their  normal  location  in  this  format.  It  will  also  be  noted  that  the 
targets  and  the  attacks  were  assigned  numbers  in  the  order  in  which 
they  were  located  in  the  input  deck;  the  pairs  of  identical  attacks 
each  were  assigned  two  numbers. 

The  target  damage  statistics  for  five  Monte  Carlo  replications  of 
the  attack  in  Case  1 are  shown  in  Fi''.  7.  The  various  annotations  will 
help  to  clarify  the  nature  of  the  various  statistics  that  are  provided. 
As  can  be  noted  in  conjunction  with  the  type  #10  targets,  "hits"  are 
not  assessed  with  CBU  weapons;  the  results  are  in  terms  of  coverage 
and  fractional  kill  (PK) . Perhaps  the  most  significant — and  typical — 
observation  that  might  be  drawn  from  these  sample  statistics  is  the 
very  substantial  uncertainty  that  must  be  attached  to  any  measures  of 
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Target,  attack,  and  trajectory  data.  Cases 
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the  "average  level  of  damage";  as  can  be  seen,  the  standard  deviation 
is  frequently  as  large,  or  larger,  than  the  mean. 

Case  2 called  for  a full  printout  of  one  trial,  but  with  the  input 
data  suppressed.  These  results  are  shown  in  Figs.  8 through  10.  The 
first  results  presented  are  the  hit  patterns  on  the  runway  (Fig.  8)  and 
main  taxiway  (Fig.  9)  as  well  as  statements  as  to  their  status.  As  will 
be  noted,  all  hits  within  an  EMD  of  the  target  are  included,  since  a 
crater  will  affect  the  surface  up  to  a distance  equal  to  the  crater 
radius.  In  the  case  of  the  taxiway  against  which  no  attack  was  spe- 
cifically planned,  it  is  apparent  that  either  an  attack  aimed  at  the 
operations  building  or  one  of  those  aimed  at  buildings  B1  or  B2  was 
wide  and  to  the  right  of  the  aim  point.  A careful  study  of  the  target 
hit  summary  (Fig.  10)  indicates  that  it  was  one  of  the  attacks  on  the 
operations  building,  since  the  hits  are  with  type  //I  weapons.  Us-’ng 
the  EMD  of  40  ft  for  structural  damage  and  a 75  ft  radius  for  the  con- 
tents of  type  #5  buildings,  the  results  show  that  this  attack  will  de- 
stroy an  estimated  15.7  percent  of  the  structure  and  damage  contents 
on  57.5  percent  of  the  floor  space  of  the  avionics  maintenance  facility 
(B3) ; on  the  other  hand,  the  AGE  facility  (B5),  which  is  similarly 
located  with  respect  to  the  intended  aim  point  near  the  operations  build- 
ing, was  missed  entirely  as  was  the  intended  target.  One  of  the  main 
aircraft  maintenance  hangers  (B2)  was  hit.  Other  dive-bombing  attacks 
were  both  wide  and  to  the  right,  placing  bombs  on  the  main  parking  ramp. 
For  the  housing  areas , there  was  no  damage  to  the  exposed  targets  in 
the  three  areas  in  this  particular  trial.  The  actual  hit  locations 
on  building  B2  and  on  the  operating  surfaces  are  shown  at  the  bottom 
of  Fig.  10.  Of  the  four  targets  for  which  hits  were  to  be  retained 
(#6,  it 7,  #11,  and  #13),  cnly  it 7 was  hit  in  this  case. 

Case  3 used  the  expected-value  mode  to  provide  estimates  of  the 
expected  numbers  of  hits  and  to  compute  a hit-density  grid.  The  latter 
is  controlled  by  specif ving  a "type  #21"  target;  as  will  be  noted  in 
Fig.  11,  the  grid  dimension  is  printed  in  the  position  normally  reserved 
for  target  angle.  Any  grids  that  have  been  specified  are  computed  first 
and  the  sample  requested  in  this  case  is  shown  in  Fig.  12;  as  will  be 
noted,  the  results  are  presented  in  terms  of  the  expected  numbers  of 
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hits  per  10,000  sq  ft,  or,  in  effect,  within  the  walls  of  a 100  * 100 
ft  building  centered  at  the  grid  position. 

Figure  13  presents  the  expected  hits  for  the  regular  target  array. 
As  a comparison  will  quickly  indicate,  the  particular  results  for  the 
single  trial  in  the  second  case  are  substantially  different  fron  the 
average  result  that  is  to  be  expected.  Furthermore,  as  comparisons 
with  the  statistical  averages  presented  in  Case  1 will  confirm,  a 
rather  large  number  of  trials  are  required  if  one  is  to  obtain  reason- 
ably accurate  estimates  of  a mean  in  problems  of  the  sort  examined  here 
using  Monte  carlo  techniques.  It  is  for  reasons  such  as  this  that  it 
was  suggested  earlier  that  the  two  modes  can  be  used  together  benefi- 
cially so  as  to  more  efficiently  obtain  both  a reasonably  good  estimate 
of  the  mean  as  well  as  useful  evidence  on  the  variability  about  that 
mean. 


Case  4 examines  the  availability  of  a minimum  operating  surface 
for  aircraft.  Since  this  type  of  problem  is  rather  complex  and  tends 
to  demand  relatively  large  numbers  of  trials,  it  will  often  be  prefer- 
able to  separate  this  examination  from  others,  as  we  have  here  in  Cases 
1 and  4,  to  avoid  unnecessarily  processing  data  relating  to  attacks  and 
targets  that  are  not  of  interest.  For  this  reason  all  targets  but  the 
two  of  interest  have  been  dropped  (see  Fig.  14),  and  all  attacks  that 
could  not  have  any  effect  should  also  be  dropped;  for  this  case  we  re- 
tained all  the  attacks,  since  none  are  extremely  distant  from  both 
surfaces. 

The  first  result  in  Case  4 is  the  trial-by-trial  record  (top  of 
Fig.  15)  of  the  total  number  of  hits  and  minimum  numbers  of  repairs 
required  for  both  of  the  targets;  there  is  no  entry  only  if  there  are 
no  hits.  This  record  is  especially  useful  if  one  wishes  to  examine 
the  distribution  of  attack  results  in  more  detail  than  can  be  done  using 
the  overall  statistical  results  them 'elves . The  latter  are  presented  at 
the  bottom  of  Fig.  15.  As  can  be  seen  by  comparing  these  results  for 
the  taxiway  with  those  in  the  previous  case,  even  25  trials  is  inadequate 
to  assure  a reasonable  estimate  of  the  mean  (which  is  actually  nearly 
twice  that  which  has  been  obtained  in  this  sample  of  25  trials). 


Expected  target  damage/  Case 
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Total  computation  time  for  these  four  cases  was  just  under  53  sec 
on  USAFE's  WWMCCS  Honeywell  6060  computer.  If  not  overlaid,  36K  words 
of  core  are  required  for  AIDA  as  presently  dimensioned,  that  is,  for 
250  targets,  50  attacks,  20  target  types,  and  10  weapon  types. 
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Appendix  A 


DETAILED  DESCRIPTION  OF  AIDA  INPUT 


The  basic  types  of  input  cards  employed  with  AIDA  are  as  noted 


below: 


control  card 

target  card;  one  per  target 

attack  card;  one  per  weapon  delivery  pass  (or  group  of 
identical  passes) 
alternate  attack  card 

effective  miss  distance  card;  one  for  each  weapon  type 
controls  sequential  cases 
terminates  overall  computation 


;r> 


The  ATT2  card  is  actually  two  cards  in  sequence  and  the  EMD  card  may 
have  up  to  three  supplementary  cards.  A detailed  description  of  the 
entries  for  each  type  of  card  is  presented  on  the  pages  that  follow. 

The  general  arrangement  of  data  on  all  basic  card  types  is  similar; 
the  card  type-name  is  placed  (left-adjusted)  in  the  first  four  columns 
and  the  data  are  listed  in  eleven  6-column  fields  between  Columns  7 and 
72.  All  data  are  read  with  a F6.0  format;  i.e.,  they  are  to  be  real 
numbers.  If  a whole  number  is  to  be  input,  it  may  be  entered  (right- 
adjusted)  in  the  field  without  a decimal  point;  the  decimal  point  is 
necessary  otherwise.  Columns  5 and  6 on  the  ATT,  ATT2,  and  EMD  cards 
are  also  used,  as  will  be  described,  and  the  name  of  the  target  complex 
being  studied  and  a name  for  each  target  may  be  included  in  Columns  73 
through  80  of  the  CONT  and  TGT  cards,  respectively;  any  alphanumeric 
names  are  acceptable. 

* 

All  linear  dimensions  should  be  in  consistent  units  (e.g.,  feet) 
and  the  target  orientation  and  the  attack  heading  entries  should  be  in 
degrees . 

If  ATT2  cards  are  to  be  used,  all  linear  dimensions  must  be  in 
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The  CONT  card  controls  the  mode  of  operation,  the  choice  of  random 
number  generator,  the  number  of  trials  (attack  replications),  and  print- 
out options;  specifies  the  minimum  clear  length  (MCL)  and  minimum  clear 
width  (MCW)  for  runway  attack  effectiveness  calculations;  and  controls 
the  runway  repair  assessment. 
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Columns  Data  Entry 

1-4  CONT 

11-12  When  0,  the  seed  for  the  random  number  generator  is  the 

same  for  all  runs.  If  greater  than  0,  the  seed  is  changed 
from  run  to  run;  if  equal  to  -1,  the  random  number  generator 
is  locked  out.  If  equal  to  -2,  the  expected-value  mode  of 
operation  replaces  the  Monte  Carlo  mode. 

13-18  Desired  number  of  replications.  Default  is  1. 

23-24  Controls  printout  options  as  follows.  If  entry  is 

5 Prints  multiple  trial  statistics  plus  a condensed 

listing  of  hits  by  trial 

4 Prints  multiple  trial  statistics  plus  a condensed 
listing  of  runway  status  by  trial 
3 Prints  multiple  trial  statistics  only 
2 Above  plus  runway  results  for  each  trial 
1 All  above  plus  hit  summary  for  each  trial 

0 All  above  plus  stored  hit  data  for  each  trial 

-1  All  above  plus  all  hits  and  target  corners 
-2  All  above  plus  all  impact  points 
30  Controls  printout  of  intermediate  information  for  piogram 

test  purposes;  should  normally  be  0.  If  set  to  greater 
than  7,  the  random  number  generator  is  locked  out.  See  the 
program  source  listing  for  the  effect  of  other  values. 

31-36  MCL  for  aircraft  operations.  (Used  to  test  if  the  runways 

are  open.) 

37-42  MCW  for  aircraft  operations.  (Used  to  test  if  the  runways 

are  open.) 


it  ft 


55-60 


61-66 


73-80 


Data  Entry 

When  entry  is  1,  runway  results  will  include  the  minimum 
number  of  craters  to  be  repaired  for  the  runway  to  meet 
the  MCL  and  HCW  criteria. 

When  the  entry  is  1,  a plot  of  all  impact  points  will  be 
included  for  all  closed  runways  (if,  also,  the  printout 
option  entry  in  Columns  23  and  24  is  less  than  3) ; when 
the  entry  is  2,  impact  plots  are  provided  for  each  runway 
whether  or  not  it  is  closed. 

The  distance  that  the  "minimum  runway  rectangle"  is  to  be 
shifted  laterally  in  checking  for  an  adequate  section;  the 
default  value  is  5. 

The  distance  along  the  runway  that  the  minimum  runway 
rectangle  is  to  be  shifted  in  checking  for  an  adequate 
section;  the  default  value  is  250. 

A name  can  be  entered  here  for  the  entire  target  complex 
and  it  will  appear  in  the  heading  of  the  output  listing. 
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Each  TGT  card  designates  the  location,  size,  and  orientation  of 
a rectangular  target. 


Columns 


Data  Entr'v 


:■! 


m 


■"  Ss’"'§i 


13-18 


19-24 


25-30 


31-36 


41-42 


The  X-coordinate  of  the  westernmost  corner  of  the  target. 

The  Y-coordinate  of  the  westernmost  corner  of  the  target. 

If  a target  boundary  runs  exactly  north-south,  the  X and 
Y coordinates  of  the  southwestern  corner  should  be  specified. 
Target  dimension  along  the  boundary  running  northeast  (or 
north)  from  the  X and  Y coordinates  of  the  reference  corner 
specified  in  the  two  previous  fields. 

Target  dimension  along  the  boundary  running  southeast  (or 
east)  from  the  reference  corner. 

Heading  in  degrees  of  the  northeast  (or  north)  heading 
boundary  of  the  target  (along  the  dimension  specified  in 
Columns  19  to  24).  (Meaning  varies  for  target  type  #21; 
see  below. ) 

Target  type.  Targets  may  be  grouped  into  up  to  10  (or 
20)  different  categories  with  like  vulnerabilities.  This 
entry  is  used  in  conjunction  with  the  effective  miss  dis- 
tance on  the  EMD  card.  Target  type  #1  is  restricted  to 
runways  and  taxiways  that  may  be  used  for  flight  operations; 
there  will  be  no  more  than  5 targets  of  this  type.  Enter- 
ing a 21  for  target  type  actually  acts  as  a signal  (but 
only  in  conjunction  with  the  expected-value  mode)  directing 
that  a 17  x 17  grid  of  hit-density  values  be  tabulated  over 
a square,  the  southwest  corner  of  which  is  entered  in 
Columns  7 to  12  and  13  to  18.  In  this  case,  entries  in 
the  third,  fourth,  and  seventh  fields  have  no  meaning. 

Unless  a different  value  is  entered  in  Columns  31  to  36 
(preferably  a number  divisible  by  16) , the  default  dimen- 
sion of  the  square  is  4000,  for  a grid  increment  of  250. 
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Columns  Data  Entry 

There  may  be  one  or  ir  jre  target  type  if 21  cards,  and  they 
may  be  intermingled  .th  normal  target  cards;  however,  when 
present,  one  of  the  type  if 21  cards  r/rust  be  the  last  target 
card  entered  for  a case. 

48  If  greater  than  0,  all  hit  locations  will  be  saved  (and 

printed  when  entry  in  Column  24  of  the  CONT  card  is  0 or 
less) . 

73-80  A name  or  number  for  the  target  (any  alphanumeric)  may  be 

entered  here.  This  name  as  well  as  the  sequence  number 
that  is  assigned  automatically  will  appear  for  target 
identification  in  the  output  listing. 


Ir 


. V4 


&&%**•»&* 

~~~ — — saas — 


,’,.  *’•.-.  P ■'  , , 
w*  « ; ' raUfrV  ^.-vw?v>  g* A 


The  ATT  card  specifies  the  parameters  of  each  weapon-delivery 
pass.  Inputs  required  are  the  attack  heading  (measured  from  north  in 
the  coordinate  system  used  to  specify  the  targets),  the  desired  mean 
point  of  impact  (DMPI)  for  a single  weapon  or  for  the  middle  of  a stick 
of  weapons,  the  aiming  error  expressed  as  REP  and  DEP,  the  ballistic 
error  of  the  individual  weapons,  the  number  of  weapons  to  be  delivered 
in  the  pass,  the  stick  length,  and  the  weapon  type  (related  to  the  ef- 
fective miss  distance  on  the  EMD  card) . 


Columns 


Data  Entr-s 
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5-6  Total  number  of  passes  with  the  following  characteristics; 

default  = 1. 

10-12  Attack  heading  in  degrees  from  north. 

13-18  The  X-coordinate  of  the  DMPI  of  a sirgle  weapon  or  the 

middle  of  a stick  of  weapons. 

19-24  The  Y-coordinate  of  the  DMPI  as  above. 

25-30  The  REP. 

31-36  The  DEP. 

37-42  Ballistic  dispersion  In  range  of  individual  weapons  (R-DISP) . 

43-48  Ballistic  dispersion  in  deflection  of  individual  weapons 

(D-DISP) . Default  value  is  R-DISP. 

49-54  The  number  of  weapons  in  the  stick. 

55-60  The  length  of  the  stick  (the  distance  between  the  first  and 

last  weapon  of  the  stick  in  the  absence  of  dispersion) . 

61-66  The  weapon  type  (used  in  effectiveness  calculations  together 

with  EMD  and  target  type) . An  entry  is  required  (an  integer 
from  1 to  10);  otherwise  hits  will  not  be  recorded. 

67-72  Probability  of  arrival  at  target;  default  = 1.0. 


AT  1’2 


The  ATT2  card  should  be  used  in  place  or  the  ATT  card  when  the  user 
wishes  assistance  with  trajectory  calculations.  When  this  card  is  used 
the  user  expresses  the  3ttack  in  terms  of  speed,  altitude,  dive  angle, 
intervalometer  settings,  etc.,  and  a special  subroutine  converts  these 
inputs  to  those  demanded  on  the  ATT  card.  The  conversion  procedure  is 
the  JMEM/AS  Open  find  Method  Zero  as  outlined  in  the  Users  Manual  for 
JMEM/AS  Open-End  Methods , WANG  Labs.,  Inc.,  Tewksbury.  Mass.,  August  1974. 

Both  ATT  and  ATT2  type  cards  may  be  used  in  the  same  run;  the  order 
of  entry  is  of  no  importance.  When  ATT2  cards  are  used  the  input  data 
will  be  reproduced  as  submitted,  as  well  as  being  tabulated  in  the  nor- 
mal manner,  after  conversion. 

Data  input  with  the  ATT2  procedure  require  two  cards.  The  first 

card  is  labeled  ATT2  in  the  first  4 columns  and  has  input  similar  to 

that  on  an  ATT  card  (all  fields  are  read  with  a F6.0  format);  a second 
unlabeled  card  is  mandatory  following  each  ATT2  card.  The  format  for 

both  cards  follows.  When  these  cards  are  used,  all  linear  dimensions 

in  the  input  data  will  be  in  feet. 


Columns 


Data  Entry 
ATT  2 

Total  number  of  passes  with  the  following  characteristics; 
default  = 1. 

Attack  heading  in  degrees  from  north. 

The  X-coordinate  of  the  DMPI  of  a single  weapon  or  the 
middle  of  a stick  of  weapons. 

The  Y-coordinate  of  the  DMPI  as  above. 

The  CEP  in  the  normal  plane  in  mils,  or,  if  DEP  is  specified, 
a constant  which,  when  divided  by  the  sin  of  the  impact 
angle,  gives  the  REP,  in  mils. 

The  DEP  in  mils  (if  omitted,  CEP  controls). 

Ballistic  dispersion  in  mils . 

The  number  of  weapons  in  the  stick. 

The  weapon  type. 

Probability  of  arrival  at  target;  default  = 1.0. 


555 


The  data  format  for  the  second  card  of  each  ATT2  pair  is  as  noted 
below  (this  card  is  used  with  a 6F6.0,  3F6.3  format).  Typical  ballis- 
tic data  required  for  thi * card  are  noted  in  Table  A-l. 
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Columns  Data  Entry 

7-12  Aircraft  velocity  (kn) . 

13-18  Release  altitude  of  last  bomb  (ft). 

19-24  Dive  angle,  at  release  (deg) . 

25-30  Terminal  velocity  of  weapon  (cluster)  or  first  leg  of  a 

high-drag  bomb  (ft/sec)  (VT1  in  JMEM) . 

31-36  Terminal  velocity  of  a cluster  bomblet  or  a high-drag  bomb 

(ft/sec)  (VT2  in  JMEM). 

37-42  ProbabJe  error  in  estimating  and  correcting  for  wind  ef- 

fects (ft/sec). 

43-48  Cluster  opening  time  or  fin  opening  time  for  a high-drag 

bomb  (ms),  or  cluster/fin  opening  altitude  (ft).  (a  decimal 
point  is  mandatory  when  altitude  is  input.)  (TD  or  in 
JMEM.) 

49-54  Intervalometer  setting  (ms) . 

55-60  Dispenser  intervalometer  setting  (ms)  (0  for  clusters). 


Illustrative  values  are  noted  on  Table  A-l. 
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Table  A-l 

TYPICAL  BALLISTIC  PARAMETERS 


Weapon 

, 

VT 

(fps) 

Ballistic  Parameter 

VT, 

(fps) 

TD  or  Hf 

Mk-81  Mod  1 

185C 

0 

0 

Mk-81  SE 

1100 

208 

300  ms 

Mk-82  Mod  1 

1900 

0 

0 

Mk-82  SE 

1200 

240 

350  ms 

Mk-83 

2500 

0 

0 

Mk-84 

2850 

0 

0 

M-117  Unretarded 

1950 

o 

0 

M-117  Retarded 

900 

168 

300  ms 

M-118 

2450 

0 

0 

AN-M64A1 

1600 

0 

0 

AN-M65A1 

2000 

0 

0 

Mk-36  DST 

1200 

240 

330  ms 

CBU-38 

*50 

o 

0 

CBU-52B/B 

1000 

230 

Variable  Altitude  (ft) 

CBU-58/B 

950 

1 

! 215 

Variable  Altitude  (ft) 

SOURCE:  //set’s 

WANG  Labs.,  Inc., 


Manual  for  JMEM/AS  Open-End  Methods - 
Tewksbury,  Mass.,  August  1974. 


The  EMD  card  is  optional  and  provides  information  regarding  weapon 
performance  against  the  various  types  of  targets.  The  entries  for  this 
card  are  different  for  point-impact  weapons  and  for  CBU-type  munitions. 
For  point-impact  weapons,  a hit  is  assessed  for  any  impact  within  a 
distance  of  EMD  from  the  target.  For  CBU  munitions,  the  EMD  card  is 
used  to  specify  the  dimensions  of  the  rectangular  bomblet  pattern. 

The  methods  for  expressing  weapon  coverage  also  differ  for  the  two 
types  of  munitions.  For  point-impact  weapons  the  EMD  is  also  used  as 
the  weapon  kill  radius,  and  coverage  is  determined  as  that  fraction  of 
the  target  area  chat  is  covered  by  a circle  of  that  radius. 

For  point-impact  weapons  (GP  bombs  or  PGMs)  the  entries  are: 


Columns 


Data  Entry 


Enter  1 if  data  are  to  be  entered  for  20  target  types. 


Enter  1 if  data  on  weapon  reliability,  p^,  or  effective 


11-12 


13-18 


19-24 


67-72 


kill  radius  for  this  weapon  type,  are  to  be  entered  (on 
the  following  card) . 

Weapon  type  (used  in  conjunction  witn  Columns  61  to  66  on 
ATT  card) . 

EMD  for  point-impact  weapons  versus  target  type  it  1. 

EMD  versus  target  type  it 2. 

EMD  versus  target  type  it  10. 


If  the  weapons  are  CBU-type  munitions,  use  the  following  entries 


on  the  EMD  card. 


Columns  Data  Entrv 


Enter  1 if  data  are  to  be  entered  for  20  target  types. 


When  more  than  i.0  target  types  are  involved,  the  EMD  data  and,  if 
specified,  the  supplemental  coverage  data  (see  p.  48),  f .•*.  target  types 
it 1 1 through  it 20  are  entered  in  ten  6-column  fields  from  Columns  13  to 
72  on  cards  that  immediately  follow  the  EMD  card  (and  supplementary  card) . 
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Columns 


Data  Entrv 


11-12 


13-18 

19-24 


Enter  1 if  data  on  weapon  reliability  and/or  on  kill  prob- 
abilities are  to  be  entered  for  any  target  type  on  the  fol- 
lowing card. 

Weapon  type  (used  in  conjunction  with  Columns  61  to  66  on 
ATT  card) . 

Enter  CBU  pattern  length  as  a negative  entry. 

Enter  CBU  pattern  width  as  a positive  entry. 
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SUPPLEMENTAL  CARD  FOR  WEAPON  RELIABILITY  AND  COVERAGE  FACTORS 

If  a 1 is  entered  in  Column  6 of  an  EMD  card,  a supplemental  card 
must  be  included  next  with  the  weapon  reliability  and  a set  of  entries 
for  the  several  target  types.  Note  that  this  card  is  not  identified, 
but  one  must  follow  each  EMD  card  that  has  an  entry  in  Column  6.  If 
a 1 is  entered  in  Column  5 of  an  EMD  card,  as  well  as  in  Column  6,  a 
second  supplementary  card  is  required  for  target  types  #11  through  # 20; 
this  card  is  the  fourth  of  four. 

All  entries  on  these  cards  are  optional;  the  default  value  for 
reliability  is  1.0.  If  an  entry  is  made  in  any  of  the  last  10  (20) 
fields  and  it  is  not  greater  than  unity,  it  is  taken  as  the  user  esti- 
mate of  the  p^  for  that  particular  weapon-target  combination  for  either 
point-impact  weapons  or  CBU  munitions.  For  point-impact  weapons,  an 
entry  that  exceeds  unity  is  taken  as  an  additional  kill  radius  and 
another  coverage  fraction  is  determined  as  that  fraction  of  the  target 
area  that  is  covered  by  a circle  of  that  radius,  given  a hit  within  EMD 
of  the  target.  Thus,  when  there  are  entries  on  the  supplemental  card 
for  certain  target  types,  coverage  fractions  are  computed  both  for  the 
corresponding  value  of  EMD  as  well  as  for  the  value  on  the  supplemental 
card. 


Columns  Data  Entry 

* 

7-12  Reliability  of  this  weapon  type;  default  = 1.0. 

"t* 

13-18  p^  or  kill  radius  for  this  weapon  versus  target  type  #1. 

19-24  p^  or  kill  radius  for  this  weapon  versus  target  type  #2. 


67-72  p^  of  this  type  weapon  versus  target  type  #10. 

Entries  for  target  types  #11  through  #20  on  a second  supplemental 
card  will  be  in  the  ten  6-column  fields  between  Column  13  and  Column  72. 


Since  these  entries  are  read  with  an  F6.0  format,  the  decimal 
point  rust  be  included. 

1" 

Only  for  point-impact  weapons. 


REDO 


The  REDO  card  is  used  to  terminate  the  input  for  one  case  and 
initiate  a new  case  with  some  or  all  of  the  previous  inputs,  as  de- 
scribed earlier. 

Columns  Data  Entry 

1-4  REDO 

7-12  Number  of  prior  targets  to  be  retained.  All  will  be  re- 

tained if  there  is  no  entry.  Use  a negative  entry  if  none 
are  to  be  retained. 

13-18  Number  of  prior  attacks  to  be  retained.  All  will  be  re- 

tained if  there  is  no  entry.  Use  a negative  entry  if  none 
are  to  be  retained . 

39-24  An  entry  of  unity  suppresses  the  input  listings  for  targets 

and/or  for  attacks  and  weapons  if  no  changes  have  been  made 
in  these  data  sets  from  the  prior  case. 
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END 


An  END  card  must  be  included  at  the  end  of  all  data  entry  cards 


Columns 


Data  Entry 


Lateral  distance  the  minimum  runway  rectangle  is  shifted 
in  checking  for  an  adequate  section. 

ITRIAL  Number  of  the  current  trial. 

KCBU  Switch;  set  to  unity  if  any  weapons  are  CBUs. 

KPTI  Switch;  set  to  unity  if  any  weapons  are  the  point-impact 

type. 

KTEST  Index  controlling  variety  of  debugging  printout  options. 

LIST  Switch;  when  set  to  unity,  target  and/or  attack  input  lists 

are  suppressed  when  unchanged. 

MCL  Minimum  adequate  length  for  required  runway. 

MCR  Switch;  set  to  unity  when  runway  availability  is  to  be 

checked. 

MCVJ  Minimum  adequate  width  for  required  runway. 

MODE  Index  controlling  mode  of  operation  (see  pp.  9 and  38). 

MTT  Largest  target  type  number  in  the  target  array. 

NA  Total  number  of  weapon-delivery  passes . 

NAM  Maximum  permissible  number  of  weapon-delivery  passes. 

V 

ND  Number  of  types  of-  weapons  in  overall  attack. 

NHITD  Switch;  set  to  unity  when  the  expected-value  mode  is 

specified. 

NJMEM  Number  of  weapon-delivery  passes  that  required  trajectory 

calculations . 
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NPLOT 

NPRINT 

NREDO 

NREP 

NSAVE1 

NSAVE2 

NST 

NSTAT 

NT 

PTM 

NTRIAL 


Switch;  set  to  1 or  2 if  runway  impact  plots  are  desired. 

Index  controlling  results  output  (see  pp.  9 and  38). 

Switch;  set  to  unity  if  an  additional  case  is  specified. 

Switch;  set  to  unity  when  repair  requirements  are  to  be 
assessed. 

Number  of  targets  to  be  retained  for  a subsequent  case. 

Number  of  weapon-delivery  passes  to  be  retained  for  a sub- 
sequent case . 

Maximum  number  of  targets  for  which  hits  can  be  stored. 

Cumulative  number  of  trials  in  which  the  minimum  runway 
was  available. 

Total  number  of  targets. 

Maximum  permissible  number  of  targets. 

Total  number  of  trials  specified. 
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KEY  ARRAYS 

AMD(I,J,K) 

Weapon  effectiveness  data. 

I 

Weapon  type. 

J 

Target  type. 

K = 1 

Effective  miss  distance. 

2 

Effective  damage  radius  or  probability  of  kill. 

ATT(I, J) 

Strong  array  for  weapon-delivery  data. 

I 

Weapon-delivery  pass  number;  numbered  internally  in 
order  of  entry. 

J = 1 

Heading  (deg) . 

2 

X coordinate  of  desired  mean  point  of  impact. 

3 

Y coordinate  of  DMPI. 

4 

Range  error  probable  of  DMPI. 

5 

Deflection  error  probable  of  DMPI. 

6 

Dispersion  in  range  (ground  plane) . 

7 

Number  of  weapons  released  in  pass . 

8 

Length  of  stick  (in  ground  plane). 

9 

Weapon  type. 

10 

Dispersion  in  deflection. 

11 

Probability  attacker  arrives  at  target. 

CBUHT  ( J , K) 

Impact  coordinates  of  the  centroid  of  the  J th  CBU 
pattern. 

K = 1 

X coordinate. 

2 

Y coordinate. 

COV(L) 

Fraction  of  target  L covered  by  one  or  more  CBU  patterns. 

HIT (I, J,K) 


I 

J = 1 
2 
3 
K 

HITR(I,J,K) 

I,J,K 

IR(N) 

IZONE(K,J) 

J = 1 
2 

MH1T(K) 

MSTAT(J) 

J = 1 
2 

3-8 

MTYPE(I) 

NCBU(L) 


Storage  array  for  hit  locations  on  specified  targets. 

1th  of  those  targets  for  which  hit  data  are  to  be 
stored. 

X coordinate. 

Y coordinate. 

Weapon  type. 

Number  of  hit  on  the  1th  target. 

Storage  array  for  hit  locations  on  type  #1  targets 
(i.e.,  runways  and  taxiways). 

See  HIT(I, J,K) . 

Switch;  set  to  unity  if  the  N th  weapon-delivery  attacker 
fails  to  reach  target. 

Denotes  which  of  the  ordered  targets  (see  TO)  fall  in 
the  Kth  target  zone. 

Lowest  numbered  target  in  the  K th  zone. 

Highest  numbered  target  in  the  K th  zone. 

Target  number  of  the  K th  target  for  which  hit  location 
data  are  to  be  stored. 

Storage  array  for  accumulating  trial  results  of  runway 
availability  tests. 

Minimum  number  of  repairs  required  to  open  a minimum 
runway . 

Square  of  J = 1,  above. 

Not  used. 

Index  that  specifies  whether  or  not  supplementary  data 
are  to  follow  the  EMD  card  for  weapon  type  I. 

Number  of  CBU  weapon  patterns  that  cover  all  or  part 
of  target  L. 


NHIT(L) 

NRW(I) 

P(L,K) 

K = 1 


STAT (L , J ) 
L 

J = 1 
2 

3 

4 

5 

6 

7 

8 

STAT2(I, J) 

i 

I 

J = 1 
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Number  of  hits  on  target  L;  by  both  point-impact  and 
CBU  weapons. 

Target  number  of  the  1th  runway  entered . 


Expected  fraction  of  target  L that  is  covered  by  the 
effects  of  point-impact  weapons;  or  probability  of 
kill  of  target  L due  to  point-impact  weapons. 

Probability  of  kill  of  target  L due  to  CBU  weapon 
patterns . 

Storage  array  for  accumulating  trial  results. 

Target  number. 

Number  of  hits  by  point-impact  weapons. 

Square  of  J = 1,  above. 

Trials  with  at  least  one  hit. 

Fractional  coverage  by  CBU  weapons. 

Square  of  J = 4,  above. 

Expected  fractional  coverage  by  point-impact  weapon; 
or  probability  of  kill.  ’ 

Probability  of  kill  by  CBU-type  weapons. 

Unused. 

Storage  array  for  accumulating  trial  results  for  targets 
of  a given  type. 

Target  type. 

Fraction  of  the  targets  of  type  I that  received  at 
least  one  hit. 

Square  of  J = 1,  above. 

Storage  array  for  target  data. 
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Target  number;  numbered  internally  in  order  of  entry. 


H 

It 

r? 

X 

coordinate 

of  westernmost 

2 

Y 

coordinate 

of 

corner 

#i. 

3 

X 

coordinate 

of  corner 

#2. 

4 

Y 

coordinate 

of 

corner 

n. 

5 

X 

coordinate 

of 

corner 

#3- 

6 

Y 

coordinate 

of 

corner 

». 

7 

X 

coordinate 

of 

corner 

H. 

OO 

Y 

coordinate 

of 

corner 

#4. 

9 Heading  of  northeast  target  leg. 

10  Target  type. 

11  Switch;  hits  stored  when  reset  to  unity. 

12  Dimension  oi  NE  target  leg. 

13  Dinension  of  SE  target  leg. 

T0(I,J)  Target  order  array  in  which  targets  are  ordered  accord- 

i-  g to  increasing  values  of  the  sum  of  the  coordinates 
of  the  western  corner. 

I Ith  target  in  tho  ordered  array. 

j*l  Value  ot  (X+Y)  for  five  Ith  ordered  target. 

2 Number  of  the  target  as  initially  entered. 

VPNREL(l)  Reliability  of  weapon  type  I. 
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Appendix  C 

PROGRAM  LISTING 


1. 

c 

MflN-4nA  41R.3AS-  OM. 

05  ASSESSMENT  90  3:L 

km 

Cn*,'ON  /<  SJAYO/  tit  ( 250. 131  * M T f 53 . 11 > , A*0( 10, 20 ,2 > . T0<250,2) , 

3. 

XI  ZJN5(  50,?  l.NHIT  ,2501  ,**H(  T (20)«hIT(20,3,24),YRW(  SI  ,H1 TP( 5, 3,250) 

4. 

X , O{25o,3l,C1Vl240) 

,*TYPi  til)),  N'A  96(250,21  »WPNRcL ( 10 1 , NCB’Jl  250 ) 

5. 

cn^n'./nT/vT  ,9t,N0, 

NTM,KT£jT,«C8,9C2,9Cl,400e,NPFlMT,NAM,NST,MTT 

t . 

C3K5»  ■»*./*»*  TS/mV  UL. 

1TR,21,9ST1T,ST»T(250,81,STAT2(20,5),MSTAT(8I 

7. 

C0“''3N  /CPV|T ,«./  i=- 

P.NPL'-  T.l'ifc,INL,NSAVvl,NSAVe2.UST,,iJ9?9 

B. 

C09MJH  • HlT‘9  / 9:'!T0,  M!;>J 

9. 

CCMM3N  /CBUHIT/  C.H'Jn*l  200,21  , l<(50),  KC8>J,  KPTl 

10. 

NT  * 0 

11. 

r 

*****  >1T“  ts  TH*  •tAX(yt)''  NUKBV8  £F  TARGETS  *««•■ 

1?. 

/■ 

uo&YS:  rr.T,  :o,  nhit,  p,  ;jv,  9A>*c,  ncbu.  STat 

13. 

NTM  =250 

14. 

t 

*****  YA«  IS  TH.:  MAX(W 

U*  MUPpip  ,'.F  attacks 

i . 

c 

'iRS  AY  St  ATT, 

It 

16. 

MA”  *40 

17. 

C 

8 f*  TARCcTS  nna  „H  (CM  HtTS  CAN  BE  ST3F63. 

18. 

c 

'.RHYS:  HIT,  *h;t 

10. 

VST  * ?0 

20. 

r 

<*•**  t;  C h/  NG f any  Of 

THF  PRJ  CE'- 1 VO  3 1 “-NS  I CMS,  MAKE  THE  APPROPRIATE 

21. 

r 

•»•**  CHANGES  19  TMC  1 

8 u AYS  A 13  THcN  ChANGc  THs  LIMITING  VALUE. 

22  . 

NA  » 0 

23. 

NO  * 0 

24. 

«C°  * 3 

25. 

O 

r— 

ii 

•-* 

24. 

4T2ST  = C 

27. 

NRtOl  * 0 

28. 

PO  5 1 * 1.10 

29. 

“TYOiin  * 0 

30. 

KP'IVLIII  * 1.0 

31. 

no  5 .1  * i,?9 

32. 

or  5 K * 1,2 

33. 

5 

apc(i  ,j,k j * 3.o 

34. 

10 

ITRI4L  T o 

35. 

00  12  1*1.5 

36. 

12 

NRw<  ()*0 

37. 

CAU  INPUT 

38. 

IF  l-.HITr  .=  ).  n 

GO  Tf  40 

39. 

IF  INT'-IAl  ,LT.  21 

GC  TC  25 

40. 

NSTA*  * 0 

41. 

00  15  1 * 1.4.T 

42. 

00  15  N * 1,5 

43. 

15 

STATU, *i)  * j.o 

44. 

no  20  v = l ,23 

45. 

DO  20  \ * 1 ,5 

46. 

20 

ST*T2  (u,4I1  = 0.3 

47. 

25 

00  30  1*1.  NST 

48. 

30 

**M  I T C 1 1 * 0 

49. 

40 

CONTI NU5 

50. 

DO  45  V = 1.8 

51. 

45 

**STAT(M  * 0 

52. 

CAU  TGTDI9 

53. 

«F  (MHJTr  .fO.  01 

GO  TC  60 

54. 

IF  (TGTINT.IJ)  ,fO. 

21.1  CAU  GAIDEN 

55. 

DO  50  l * l,NT 

5t . 

COVU  ) * 0.0 

57. 

P(T » 1 ) * 0.0 

5f  . 

P(l , 2)  * 0.0 

59. 

50 

PU,31  * 0.0 

60. 

CAU  5XPHIT 
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61. 

CO  T'*  230 

62. 

40 

IF  (KPU  .r 0.  0)  Tu  100 

6?. 

CALL  TCtrm 

64. 

CALL  T57.TJ* 

65. 

100 

CONTINUE 

66. 

(TRIAL  * ITAUL  6 l 

67. 

90  196  It  ' !•  NST 

66. 

TO  105  12  * 1.3 

65. 

PO  135  13  * 1.  25 

70. 

105 

MITtI  1*12*131  * 0.0 

71. 

no  no  r*i*  nt 

72. 

COV( I ) = 0.0 

73. 

6(1,11  * o.o 

74. 

PII,A)  * 0.0 

75. 

PIT, 3)  * 0.0 

76  . 

NCBU(I)  * 0 

77. 

no 

NHI T( I 1 *0 

71. 

DO  115  11  * l,  5 

79. 

no  114  12  * 1.  3 

60. 

TO  115  13  » 1*  250 

.31. 

115 

MIT4< 11,12,131  * 0.0 

82. 

CALL  BOMB 

83. 

IF  (KCMJ  . = 0.  11  CALL  LBU 

84. 

IF  («C5  .60.  0)  V'  »**  130 

85. 

CAt  L CHrCFF 

8<  . 

130 

CONTI  Ml? 

67. 

IF  (NTP1AL  .LT.  ?>  Gr  TO  I/O 

86. 

DO  140  I * l,  NT 

86. 

A!C  » NHITltl  - MCttUllI 

90.  . 

STATU, 11  * STATU,  U *410 

91. 

STATU,  2)  * STA  T( ; . 2 1 A !T»A1P 

92. 

IF  ( AID  .r.T.  0.01  STATU, 31  * STATU 

93. 

STATU, 4)  = STATU,  4)  * CCVUl 

9*. 

STATU, 51  * STATU, 51  ♦ CCV(U*CTV(I) 

95. 

STATU, 61  * STaT(.,4|  ♦ P(t,U 

96. 

STATU, 7}  * STATU, 7)  ♦ P1I.2I 

97. 

STATU, 81  * STAUI.HI  ♦ P(l,il 

9e. 

140 

CPNUWJF 

99. 

f>0  160  H * 1,20 

100. 

PM  » 0 

101. 

« 0 

102. 

AIDl  > 0.0 

103. 

AI03  « 9.0 

104. 

A1C4  * 0.0 

1C5. 

AI04  « 0.0 

106  . 

00  150  1 * 1.3T 

107. 

IF  (TGTl  1,131  .NF.  «>  CC  TO  153 

ioe. 

MM  = NN  ♦ 1 

105. 

NA19  - MMITd)  - MCBUIU 

110. 

AREA  x TGTU,  ;2)«TfiT(  1,131 

111. 

4IDI  x A! Cl  ♦ A°.cA 

112. 

AI03  = AI03  ♦ A^sAx P (1,11 

113. 

A 104  * 4104  ♦ l <t  AxP  1 1,21 

114. 

AIC5  » A 105  ♦ ARiAxPIl  ,31 

115. 

IF  INAIC  .GT.  0)  M»  x PM  ♦ l 

116. 

150 

CONTINUE 

117. 

IF  (HA:  .FO.  01  GO  TO  160 

118. 

Air  x NN 

115. 

STA  * KM/ AID 

120. 

STAT2 (p» U x 5T(T2(K,1)  ♦ STA 

121. 

STAT2IK.2I  x $T  AT 2 ( k , 2 ) ♦ STA»STA 
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t22. 

ST  A T 1 ( *■• , 3 ) = ST  AT2I  v , 3 1 

♦ AID3/AIP1 

123. 

3TA 7 2 (M  • 4 1 = RTAT2IM.4) 

♦ AI04/A1D1 

124. 

ST AT 3 (w ,5  ) = ST  AT 2 ( m , 5 1 

♦ A1DR/AIDI 

125. 

1 60 

CONT I NUE 

126. 

1 70 

COf.’T  I N’,R 

127. 

IF  ( INT  ,r.T.  11  GC 

TP  1 BO 

120. 

CALL  PRINT 

120. 

110 

;'Or:Ti  nos 

i AC. 

! e ( N P , i : : T . i: . 5)  G"> 

10  190 

131  . 

DD  115  L = 1 , NT 

. 3 c . 

18K 

WRT  '1 (6,1001)  I T B I AL,  L 

, NHIT(L) 

133. 

190 

rfi,.r  , 

134. 

TF  ( ! Tf  I AL  .L  f . MV-  ! *L  I 

GO  TO  IDO 

13  5. 

IF  ( siTS  I AL  . ST.  1)  CALL 

STAT)  S 

136. 

r jo 

IF  (VB[.jr'  .£1.  1)  GO  Tt  10 

137. 

STOP 

138. 

1001 

FORMAT! • * , 1 T P 1 AL  * , 1 4 , • 

TGT ’ » 1 4 , * 

139  . 

END 

140. 

SUBROUTINE  !N°UT 

14)  . 

(NT cOFD  '4  LA8FL,  AN,  NAPE l , Kt'A-i, 

142. 

COMMON  /".F'AYi/  TGT  ( 250  , 

13),  ATT ( 50  , 

143. 

XI 20N?( 50 , 2 ) , NH IT ( 2^0 1 , HM 1 T (20) »HtT(2D 

H-r  T s • . I 


X ,P( 250, 3) ,f Jv (250) ,MTYPE(  ID ) , NAME ( 2 50 » > ) .WPNREL  ( 10 J ,NCBU(250) 
CCMM1\/I'jT/*|Tf,MA,‘J'),NT’<,KT6ST  .^CR.^Cw.mcl.m-JDuNPP.  INT , NAM, NST , MT T 
COPM'N/STAT$/MT.<iAL  , IT-MiL,NSTAT,  ST  AT  C 250  , ft ) , ST  4 T 2 ( 20 , 5 ) , *ST  AT  ( fl  » 
COMMCN  /CC’ITKL/  NFf P.NPlCT, JNW, INI,  NSAV!  l,N3AVF2,LIST,NJMeM 

comm:n  /cbuhit/  cBHHr( 200 ,21,  iV,(5o>.  kchu,  kpti 

CCMMCN  / HIT  OK  / WHITP,  VJfcDG 
DIMENSION  1A6£L(6),  OATAdll 

PAT  ft  LABEL  /‘TOT  * , • ATT  • , • ft  TT2  • , 1 ; MD  • , • CONT  • , • P EDO  * / 

NJM^m  =0 

l i s r i » o 

LI5T2  * 0 


'l  155. 

IF 

(NRtOP  . 

SO.  3) 

NH  ITT  * 0 ‘ 

I,  156. 

IF 

(NR EDO  . 

c).  0) 

00  TO  2 

5 157. 

C 

SEE 

N |T£  AT 

LABEL  *45 

, 

j 158. 

IF 

(NSAVE1 

.r,T.  j ) 

NT  * NSAYC1 

1 14J. 

IF 

( *J  SAVE  2 

.OT.  ()) 

NA  = NSAVC2 

i 160. 

IF 

( N SAVE  1 

.LT.  01 

NT  = 0 

! 161. 

IF 

(NSAVE2 

.LT.  0) 

NA  = 0 

162. 

IF 

( N SAVE  1 

.20.  0) 

LISTl  * LIST 

163. 

IF 

( NS. AVI  2 

.10.  0) 

LI  ST  2 = LIST 

nsa/ei  = 0 

NSAVR2  =0 

CONTINUE 

NKEOP  = 0 

READ  (5,101)  AN,  NT  Y PE  » ( DATA ( I ),l *1 , 111 , NAME  l ,NAME2 
Ni)  ENTRY  IS  RcOUI»ED  IN  CCLJMNS  5 AND  5.  IF  AN  INTEGtR  IS  FOUND 
ON  AN  ATTACK  CAPO,  THE  ATTACK  WILL  Bc  REPEATED  SO  THAT  THERE  WILL 
BE  THAT  TOTAL  NUMH£K  OF  ATTACKS  WITH  THE  STATED  CHARACTERISTICS. 

(ONE  ATTACK  IS  ASSUMED  IF  THERE  IS  NO  ENTRY.) 

IF  AN  ENTRY  IS  NOTED  ON  AN  AMD  CARD,  RELIABILITY  AND/OR  PROBABILITY 
OF  KILL  DATA  rfILL  Bt  EXPECTED  ON  THE  NEXT  CARD. 


175. 

IF 

( AN 

.EC. 

LABEL! 1 ) ) 

GO 

TO 

10 

I7t  . 

IF 

< AN 

. . 0. 

L AB  EL ( 2 ) ) 

GO 

TO 

20  '•'> 

1.7. 

IF 

(AN 

• «T  C« 

LAP  El ( 3 1 ) 

GO 

TO 

26 

s 

170. 

IF. 

(AN 

. EO. 

LABEL (4) ) 

GO 

TO 

30 

179. 

IF 

(AN 

.EC. 

LAB<=1  (5)  ) 

GO 

TO 

40 

3 

180. 

IF 

(AN 

.EO. 

L ABEL (6  J ) 

GO 

TO 

45 

3 

181. 

GO 

TO  50 

• 

NT  = NT  * 1 


10 
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183. 

USTl  « 0 

194. 

IF  ( NT  .GT.  ‘iTM ) f,C  Tl  120 

165. 

GO  12  ! * 1,2 

136. 

12 

TG7  {NT, 11  =3‘T4  (1 | 

197  . 

00  13  1 * 3,4 

186. 

13 

TGTC’CT,  !»91  * OATAtll 

189. 

I>0  14  1 * 5,7 

190. 

1* 

TGTINT.l  ♦ 41  * DATA  111 

191. 

NAME  ("IT  , l ) » NAM31 

192. 

NINE  (NT,  21  * .\‘V*52 

193. 

GO  TO  6 

1'4. 

20 

N4  * NA  ♦ 1 

195. 

l 1 $T2  « 0 

196. 

IF  I Nft  .GT.  NANI  f,D  70  130 

197. 

00  22  1 * 1,6 

196. 

22 

4TT(NA,i;  « CATAUI 

195. 

ATT  ( '14,101  * OA  T A ( 7 1 

200. 

IF  lftATA<71  .iO.  0.01  / TT  (NA,  10)  * ATT  INA,61 

?01. 

IF  (OAT AC  11  > ,C7.  0.1  CATAUU  » 1.0 

702. 

ATT  (NA,  111  * 0ATAC111 

203. 

nc  24  t t 7,9 

204. 

24 

ATT IN A, 1 1 * 7ATAC1+11 

205. 

NTT«*r  « >AX0CCNTYP*--1|,0) 

206. 

IF  CNTYCS  .50.  0»  GC  TP  6 

207. 

GO  TO  20 

208. 

2* 

NA  c ‘tA  > l 

205. 

L1ST2  * 0 

210. 

IF  « NA  .GT,  NA*4)  GO  Tt  130 

211. 

IF  (OATAClli  .63.  0.01  0ATAC111  - 1.0 

212. 

CAU  JM5M0C  NJNCN,  OAT A1 

213. 

00  27  !«  1,11 

214. 

27 

ATTINA.P  « OATAIII 

215. 

2S 

NTY0C  * MfxOJINTY«»=-l),OI 

216. 

IF  (NTYPt  .CO.  0)  GO  TO  6 

217. 

NA  • NA  ♦ 1 

218. 

DO  29  1 « 1.11 

214. 

29 

ATTCKA.t)  « ATTHnA-11,1) 

220. 

GO  Tn  2P 

721. 

30 

NO  a NO  * 1 

227. 

LIST?  * 0 

223. 

9 • OAT AC  11 

224. 

IF  CM  .GT . 101  GO  TO  140 

225. 

HTYFC  (Ml  « NTYP? 

226. 

no  31  N * 2,11 

227. 

31 

AMD(K,'e-l,ll  = DATA  ('ll 

228. 

IF  CINTYPE  .FO.  01  .CR.  CNTYPfc  .SO.  10)1  GO  TO  34 

22  9. 

REA0C5,  1141  VCP  , (nAT  Ai  1 ) , l«1.10> 

230. 

HPA^CUNl  a W# 

231. 

IF  (wk  .to.  0.01  MPN»5LCN1  « 1.0 

232. 

00  3?  '1  * 1,10 

233. 

32 

AH0CM.N.21  . DATACM 

234. 

34 

CONTINUE 

235. 

C th 

3 FOUOWNG  8 STATEMENTS  PEP*»it  20  TARGET  TYPES  MHEN  NEEDED 

236. 

IF  I "TYPE  .IT.  101  GO  TO  6 

237. 

AfcAO  (5,1141  (OITA ( II,  { « 1,111 

238. 

00  36  N "•  11,20 

239. 

36 

AM0(M,N,11  ■ OATACN-91  ^ 

240. 

IF  (NTYPf  .CT.  Ill  GO  TO  6 

241. 

AEADC5, 114)  (OATACII,  1 « 1,111 

24:  . 

00  38  N • 11,20 

243. 

38 

3hOCN,N,2}  « OATACN-9) 

\ 


yjM&'S-fttS** 
J?  ' 


ytVAW;-  **■  W1-^’  3c£*$g>*  $**$  <<*  “»  '*  ^sv^'^ajjf  >, 


„ V-'. 

'"'su *i«i.  j,  ;.‘  .-£.<'  j-^fr^^i-.sri  — ^ 


*sipn»".M,!r-4-.^»> 


' 


« 


244. 

GO  T • 6 

1 

245. 

40 

MODS  = )£  TA  ( 1 ) 

I 

246. 

If  (MQOS.EO.OJ  Cail  °STAPT(7) 

<jZ 

$ 

247. 

NTRIAL  - DATA!?) 

4 

248. 

NPRINT  = PAT  A { 3 ) 

Lj 

249. 

KTEST  = 0 AT  4 ( 4 ) 

* 

250. 

IF  (NTRIil  .IT.  2)  NTRIAl  = 1 

} 

251. 

VCP.  = 0 

% 

252. 

NH I TO  = 0 

252. 

IF  (MpDe  .SO.  -2!  ,'IH  1 TO  = 1 

l 

254. 

IF  IMHITO  .so.  1)  CATA < 5 ) = 0.0 

i 

255. 

IF  { OAT  A I 5)  .GT.  0.01  MCR  = 1 

k 

256. 

MCL  = 04TAI51 

* 

257. 

MCW  = TATAI6) 

< 

258. 

NREP  = TATA  171 

A 

259. 

NPLOT  = TATA (8 ) 

0 

260. 

INW  = OA  T A ( 9 ) 

< 

261. 

INL  = PATAdOI 

t 

262. 

IF  (INW  .SO.  0)  IOW  = 5 

j 

263. 

1=  (IM  .80.  0)  INL  = 250 

ft 

264. 

NBASS1  = NAMcl 

26*. 

NB  A S 2 2 - NAME2 

266. 

GO  TO  6 

; 

267. 

45 

NRED*1  = 1 

« 

266. 

C 

THE  F 1 3 ST  TW0  TATA  cNL  dES  CO  THE  REOG  CARP  m4y  85  LSED  TO  SPECIFY 

1 

269. 

C 

THE  NUMBcRS  Dc  AR1TR  TAPPETS  A\‘C  ATTACKS,  R ES°t  CT  I VEL  Y , TO  BE 

2 7C. 

C 

INCLUOEC  IN  THE  NSW  CALCULATION.  THOSE  NUMBERS  WILL  BE  SELECTSO 

271. 

c 

IN  RANK  3P0ER  ck DM  THOSE  IN»UT  PREVIOUSLY;  IN  NO  CASE  MA Y THE 

272. 

c 

NUMRF=  8=  LARGER  THAN  THE  NJMSEP  AVAILABLE.  IF  NC  NUMBER  IS 

273. 

c 

ENTERED,  ALL  =810“  TARGETS  AND  ATTACKS  WILL  BE  INCLUDED.  IF  A 

„ 

274. 

c 

NEGATIVE  NUMB?5  IS  ENTERED  (EG  -l),  NONE  Oc  THE  PRIOR  MEMBERS 

275. 

c 

OF  THAT  SFT  WILL  BE  TREATED. 

276. 

NSAVc  1 = DATAd  ) 

* 

277. 

NSAVE2  = DA  T A ( 2 ) 

278. 

f 

IF  THE  THIRO  c NT  BY  (COL  24)  IS  SET  TO  UNITY,  THE  TARGET  LIST 

279. 

r1 

V. 

AND/OR  THE  ATTACK/WcADON  l I ST S WILL  BE  SUPPRESSED  IF 

’ 

280. 

C 

THEY  HAVE  NOT  BREN  CHANGED. 

281. 

LIST  = 0ATM3) 

282. 

50 

CONTINUE 

283. 

WRITS  (6,111)  NTRIAL,  NPPINT,  MODS,  MCL,  INL,  MCW,  INW, 

284. 

X NRS °,  NPLOT,  KTEST 

285. 

IF  < ( NBAS  £1  .EO.  0)  .AND.  ( NBASE2  .EQ.  0))  GO  TC  55 

286. 

WRITE  (6,100)  NBASSl,  NBASS2 

287. 

56 

CONTINUS 

288. 

IF  <1  IST1  .ED.  1!  GO  Tr  65 

289. 

WRITE  (6,102) 

290. 

m.tt  = 0 

291. 

00  60  I = 1 »NT 

292. 

IF  ( (TGT (1,10)  .GT.  VTT)  .AND.  (TGT(I.IO)  .NE.  21.)) 

293. 

X MTT  = TGT ( I, 10) 

„ 

294. 

60 

WRITS  (6,112)  I,(TGT(I,J)»J=1»2)»(TGT(I,J)»J=12,13)» 

295. 

X (TGT( I , J) »J=9,11) , NAME (1,1) ,NAME ( 1,2) 

296. 

65 

IF  (LIST2  .EO.  1)  GO  TO  95 

297. 

WRITE  (6,104) 

298. 

MWPN  = 0 

299. 

KPT  I =0 

300. 

DO  70  I = 1,  NA 

301. 

IF  ( A TT ( I ,9 ) .GT.  MWPN)  MWPN  = ATT(1,9) 

302. 

IF  ( AMD (ATT (I, 9), 1,1)  .Gt . 0)  KPTI  = 1 

‘ 

303. 

70 

WRITS  (6d03)  l,  (ATT(I,J),J=l,6),ATT(I,10),(ATT(I,J)tJ=7,9) 

304. 

X , ATM!  LI) 

4 

' i,  /-,*  ;.  * ln»  ,v^,1*«f'Virt^iv*>  ■' 


T'tJ<fA^  *&ffiR*$f  *£**-?«  rt~t- o& . 


■~t&?Y**i^r<i  « jjtft  j 


j1  ’wmtwoaaw*!*  B^araT 


|x 


- «k*£l 

■-j  j||g] 


w- 

3&F  • 


IF  CIO  ,F3.  Oi  GC  t0  15 
WRITE  (6,1051  (I,  1=1, ’3) 
no  90  1=1,  MW PM 

WR1TE16, 106)  It  WPNRFL(I) ,(AMP(I, J,l) ,J=1,10) 

IF  ((MTYPE(t)  • EO.  0)  .OF.  (MTYPECI)  .EO.  101)  GO  TO  80 
WRITE  (6,1071  (AMO  t I » J * 2 f , J = 1,10) 

CONTINUE 

IF  ( MTYnF ( I ) .CT.  10)  GC  TO  90 
WRITEI6.115)  (4M0(I,J,1),  J=ll,20) 

IF  ( MTYPc ( I ) .63.  10)  GO  TO  90 
WRITE (6, 107 ) ( AMO ( I , J , 2 ) , J=ll,20) 

CONTINUE 

CONTINUE 

TARGET  TYPE  *1  IS  RES'RVcC  FOP  RUNWAYS  ANO  TAXIWAYS  (OR  OTHER 
LARGE  TARGETS  IF  MC°=0>  ANO  HIT  STORAGE  IS  PROVIOcT  FOR  250  HITS 
BUT  FOR  A MAXIMUM  OF  f TV?  TARGETS  0C  TYPE  <fl. 

NTX  = 0 

00  99  I = 1.  NT 

IF  (TGT(I,10)  .NE.  1.)  GO  TO  99 

NTX  = NTX  + 1 

IF  (NTX  .GT.  5)  GO  Tn  150 

NRW(NTX)  --  I 

CONTINUE 

RETURN 

WRITE  (6,108) 

STO° 

WRITE  (6,109) 

STOP 

WRITE  (6,110) 

STOP 

WRIT?  (6,116) 

STOP 

FORMAT ( • ',  15X, •**«**  BASE  COMPLEX  NAME  - ',244,*  ****«',/,/) 

FOFMAT(  A4,I2,11F6.0,2A'») 

FORMAT  (•  •,  20X, ‘TARGET  DATA*,/,*  NUMBER  X-DIM  Y-DIM 

X NE  LIMB  SE  LIMB  ANGLE  TGT  TYPE  STORE  BLOG  NO‘,/,/1 
FORMAT! * ‘,I4,2X,I0F10.0,F10.3> 

FORMAT('  ',/,  20X, ‘ATTACK  DATA',/,'  NUMRER  HOG  X-OMPI 

X Y-OMPI  REP  CEP  R-DISP  D-OISP  NO  WPNS  LENG 

XTH  WRN  TYPE  ARP  I V AL ’ , // J 

FCRMAT('1‘,  3 5 X , 'MISS  DISTANCES  ALLOWED  FOR  EFFECTIVE  HITS',/, 
X50X, ‘TARGET  TYPES*,/ ,20X,10(6X,I2,2X),/,'  WPN  TYPE  WPN  REL  •, 
X/,/,/) 

FORMAT ( ' ‘,I6,7X,F5.3,10F1.0.0) 

FORMAT ( * »,  2 IX,  10F 10.3) 

FORMAT ( • ♦,/,'*****  TOO  MANY  TARGETS  HAVE  BEEN  SPECIFIED  *♦**•) 
FORMAT!'  ',/,'♦****  TCO  MANY  ATTACKS  HAVE  BEEN  SFECIFIED  ****•) 
FORMAT!'  ',/,'*****  TCO  MANY  TYPES  OF  WEAPONS  HAVE  BEEN  •, 

X 'SPECIFIED  »****•) 

FORMAT  (' 1'  , 1 30( •*•)♦/» 

X 40X,'*****AIDA*****  AN  AIRBASE  DAMAGE  ASSESSMENT  MODEL',/, 

X 40X, 'DEVELOPED  8Y  THE  RAMSTEIN  OFFICE  OF  THE  RANC  CORPORATION',/ 
X , 130 ('*'  ),/, ' NO  OF  TRIALS  ',  13,'  NPRINT  ',12,'  MOOE  ',12, 

X ' MCL  ',I6,'(',I4,'»  MCW  ',14,'C,  12,')  MIN  REPAIR  ', 

X II,'  PLOT  HITS  ',11,'  TEST  '»I2»/»1 30 (•*')»/,/) 

FORMAT  (•  ',  I4,4X,7F10.0,2X,2A4> 

FORMAT  ( 6X » 1 1F6.0I 
FORMAT  (•  • ,I8X,  10F10.0) 

FORMAT  (•  ',/,'*****  TOC  MANY  RUNWAYS/TAXIWAYS  HAVE  BEEN', 

X • SPECIFIED  ***-•  ) 

END 
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SU6R0UT1NE  TGTPIM 

COMMON  /ARRAYS/  TGT(250,13I,  ATT  ( 50  , i 1 ) , AMD(10,  20  ,2  ) . T0(250,2), 
Xt  Z0NE(50,2>,NHIT(?50)  ,MHI7  120 ) ,HI  T 12  0, 3 *25  )',NF  W(  5)  ,HITR(  5,3,250) 

X , P ( 250*  3) ,C0V(250) .MTYPEt 10 J , NAME (2 50, 2 > , WPNREl ( 10) ,NCBU( 250) 
COMMON  / 1 NT  /NT  , NA ,ND»  NTM ,KT  EST ,MCR , MCW»MCL  * MQDL ,NPR  INT  * NAM , NST  »MTT 
COMMON/ STATS/NTRIAL, 1TF I AL.NSTAT,  ST  AT  ( 250 , 8 > , ST  AT  2(  20, 5 ) , MSTAT(  8) 
COMMON  / HI  TON  / NHtTD,  NR EDO 
IF  (NPRINT  .IT.  0)  WRITE  (6,104) 

10  20  I = I,  NT 

TGTDI M COMPUTES  AND  STORES  THE  LOCATION  OF  THE  QTHE«  CORNERS. 

LI  = TGT ( I,U) 

L2  = TGT( I, 13) 

THETA  = TGT  (I,  9>/57.3 


S = 

SIN  (THETA) 

C = 

COS  (THETA) 

US 

= L ITS 

HC 

= L 1*C 

L2S 

= L2*S 

L2C 

= L2Y  C 

TGT 

(1,3)  = TGT 

(I, 

i) 

4- 

LIS 

TGT 

(1,4)  = TGT 

( I. 

2) 

4- 

L1C 

TGT 

(1,5)  = TGT 

( I. 

3) 

4- 

L/C 

TGT 

(1,6)  = TGT 

(I, 

4) 

- 

L2S 

TGT 

(1,7)  = T GT 

( I, 

5) 

- 

LIS 

TGT 

(1,8)  = TGT 

(I. 

6) 

- 

L1C 

IF 

((KTEST  .GT. 

2) 

.OR 

• 

( NPR 

X WRITE  (6,102)  I,(TGT(t,K),  K=l,8) 

CONTINUE 

tF  ( NH1T0  .EQ.  1)  GO  TC  50 
NR  = 0 

00  30  I = 1,  NT 

FOR  SPECI CI EO  TARGETS  THE  TARGET  NUMBER  IS  STORED  IN  MH IT  FOR 
LATEF  REFERENCE. 

IF  l (TGT (1,11)  .LT.  1.)  .OR.  (TGT (I  , 10 ) .EQ.  1.))  GO  TO  30 
NR  = NR  ♦ 1 

IF  (NR  .GT.  NST)  GO  TC  80 
MHIT(NP)  = I 

IF  ( (KTSST  .GT.  41  .AND.  (ITRIAL  .LT.  2)) 

X WRITE  (6,101)  NR,  MHIT(Nfi) 

CONTINUE 

IF  (NR  .EQ.  NST)  GC  TO  50 
NR1  = NRH 

DO  40  1 = NR l,  NST 

MHIT(I)  = 0 
CONTINUE 

IF  ((NPRINT  .LT.  0)  .OR.  (KTEST  .GT.  2))  WRITE(6,104) 


GO  TC  80 


WRITE(6,104) 


RETURN 

WRITE  (6,103) 

STOP 

FORMAT ( * * , 1 MH I T ( * , 12,  •)  = *,12) 

FORMAT! ’ •, ‘TARGET  CORNER  : TGT  « » , 14 ,4 t 4X,F6. 0 , IX , F6. 0 > ) 

FORMAT (’O’,  « COMPUTATION  STOPPED:  HIT  DATA  SPACE  •, 

X ‘REQUIRED  FOR  MORE  THAN  "NST"  TARGETS’) 

FORMAT  (»l«) 

END 

SUBROUTINE  TGTORD 

COMMON  /ARRAYS/  TGT(250,13),  ATT ( 50 ,1 1 ) , AMD! 10, 20 ,2) , T0(250,2), 

XI  ZONE  ( 50,2), NH IT  (250)  ,MHIT(20>  ,H|  T (20,3 ,25  > ,NP.W{  5)  ,HITR(  5,3,250) 

X ,P(250,3),C0V(250),MTYPE( 10 ), NAME (250,2) , WPNREL ( 10) ,NCBU(250) 

COMMON/I  NT /NT, NA, NO, NTM, KTEST, MCR,MCW,MCL, MODE, NPRINT , NAM, NST ,MTT 
THIS  ROUTINE  CREATES  AN  ARRAY  IN  WHICH  THE  TARGET  NUMBERS  ARE 


0 


i \ 

0 4 


I 

1 


fetf 


42  7. 

C 

DRDEREC  ACCORDING  TO  INCREASING  VALUES  of  IXfY)  OF  THE  INDEX  CORNER 

vi 

42  8. 

C 

OUTPUT  VALUES  : T3(I.1»  IS  THE  (X+Yl  AND  TO ( I , 2 1 IS  THE  •ORIGINAL* 

.3 

429. 

c 

TARGET  NUMBER;  I.F.  ITS  POSITION  IN  THE  INPUT  LIST 

Ja 

* 

430. 

DO  10  1=1, NT 

$ 

431. 

c 

INITIALIZES  TO  (TARGET  CRCER)  ARRAY 

i 

432. 

IF  ( TGT  (1,9)  .GT.  45.)  GO  TO  5 

5 

433. 

T0(I,1)=TGT(1,1)«-TGT(I,2) 

434. 

GO  TO  8 

fi 

435. 

5 

Ton,  1I=TGT(1,7)+TGT(I,8) 

& 

436. 

8 

TCi  i,?)  = I 

f 

437. 

10 

CONTINUE 

j 

438. 

CU  20  J =2 ,NT 

439. 

NTEST  =0 

440. 

CO  15  K = 2 »NT 

V 

441. 

c 

REORGANIZES  THE  TJ  ARRAY  INTO  INCREASING  VALUES  OF  THE  INDEX  CORNER 

< 

442. 

I = NT-K«-2 

i 

443. 

IF  { T0( 1,11  .G£.  TO ( I -1 , 1 1 ) GO  TO  15 

i 

444. 

NTFST=1 

445. 

T=TOt 1-1,1) 

j 

446. 

TN=TO ( I -1  ,2) 

447. 

TO ( I -1, 1) =T0( 1,1) 

448. 

T0( 1-1,2) =T0( I, 2) 

i. 

449. 

TO ( I , 1 ) =T 

450. 

TO ( I , 2) =TN 

451. 

15 

CONTINUE 

J* 

452. 

IF  ( NTtST  .GO.  0)  GO  TC  25 

i 

453. 

20 

CONTINUE 

‘ 

454. 

25 

CONTINUE 

455. 

IF  (KTEST  .LT.  3)  GO  TO'  40 

456. 

00  30  I = 1,  NT 

4 

457. 

NTO  = TO (1,2) 

458. 

30 

IF  (KTFST  -GT.  5)  ViR  I Tf  (6,101!  I , NTO 

459. 

40 

CONTINUE 

460. 

NT1  = NT  +■  l 

461. 

DO  50  I = NT1,  NT M 

462. 

DO  50  J = 1,2 

463. 

50 

TO ( I , J ) = 0.0 

464. 

RETURN 

465. 

101 

FORMAT  (•  •,'  RANK  »,!?,•  TARGET  * ',  13) 

, 

466. 

END 

’ *' 

467. 

SUBROUTINE  TGT Z ON 

468. 

COMMON  /ARRAYS/  TGTI250.13),  ATT ( 50 , 1 1) , AMOJ 10,  20 , 2) , T0(250,2), 

469. 

XI  ZONE (50,2  ),NH  IT (250)  ,MHIT (20) ,HI T ( 20, 3 ,25 > ,NRW(  5 ) ,H  I TR(  5,3,250) 

470. 

X , P ( 250, 3 ) ,C0v(250),MTYPE(10),NAM£(250,2), WPNREL ( 10) ,NCBU(250) 

471. 

COMMON /I  NT /NT , NA ,ND, NTM ,KT  EST ,MCR ,MCN,MCL , MODE ,NPR INT ,NAM, NST ,MTT 

- 

472. 

DIMENSION  I Z( 250) 

473. 

c 

TGTZON  IDENTIFIES  TARGET  L^CAT  I C*1  IN  TERMS  OF  ITS  IZONE  SO  THAT  THE 

474. 

V. 

SUBSEQUENT  SEARCH  PROCESS  CAN  BE  REDUCED.  CONSIDER  THE  ENTIRE 

475. 

c 

TARGET  AREA  MAPPED  BY  LINES  CF  CONSTANT  (X4-Y).  ALL  TARGETS  WITH  • 

“ 

476. 

c 

•INDEX  CORNER*  (X+Yl  FALLING  INTO  THE  K TH  500  FOOT  SEGMENT  OF 

‘ 

477. 

c 

(X+Y)  ARE  IN  THE  K TH  ZONE . THE  ORDERED  INDEX  NUMBER  FOR  THE 

478. 

c 

TARGET  WITH  THE  LOWEST  (X+Y)  IN  TnE  ZONE  IS  I ZONE < K ,1 ! ; THAT  WITH 

479. 

c 

THE  HIGHEST,  IS  I ZONE (K, 2).  IF  THERE  ARE  NO  TARGETS  IN  A ZONE,  THE 

480. 

c 

IZONE  VALUES  ARE  BOTH  EQUAL  TO  THE  INDEX  NUMBER  OF  THE  LAST  TARGET 

481. 

c 

AS  I ZONE! K-l , 2 ) ). 

482. 

DO  10  1 = 1, NT 

483. 

A A = TQ( I , 1 ) /500. 

t 

484. 

10 

um  = AA 

485. 

I ZONE (1,1)  = 0 

; 

486. 

IZONE (1,2)  = 0 

487. 

IF  (12(1)  .NE.  0)  GO  TO  14 

+ 

l ^ 

IS 

fete 


"j 


- . / '* 


\.  ' > ^ * 
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488. 

K = 1 

48  9. 

I ZONED  ,1)  = 1 

490. 

1 ZONE (1,2)  = 1 

491. 

00  12  I « 2, NT 

492. 

IF  IIKII  .GT.  0)  SO  TO  22 

493. 

U 

I ZONE  (1,2)  =■  I 

494. 

C 

ALU  HITS  IN  7 ON 6 ¥ i Hr^c. 

495. 

GO  TO  30 

496  . 

14 

00  16  K = 2,50 

497. 

IF  (mil  .EQ.  (K-l) > GO  TO  18 

498. 

I ZONc  (K  , t ) , o 

499. 

16 

I ZONE (K , 2 ) « 0 

500. 

18 

CONTI NUF 

501. 

I ZONE (K , 1 ) * 1 

502. 

IZ0N5 IK-1 ,2 ) = 1 

503. 

I ZONE  IK  , 2 ) * l 

504. 

C 

AT  THIS  POINT  < IS  ZONE  OF  FIRST  HIT 

505. 

00  20  I « 2, NT 

506. 

;f  i iz: i j .gt.  (K-i 1 1 go  to  2? 

507. 

20 

1 ZONE (K  « 2 ) = I 

508. 

r 

IN  TR&NSe5P  TO  '22'  K IS  FIRST  0CCUPI=0 

509. 

c 

IN  ( K + 1 ) ZON*. 

510. 

22 

CONTINUE 

511. 

N = I 

512. 

00  28  I = N , NT 

513. 

c 

SKIP  T->  26  IF  HIT  IN  ZONE  CF  PRIOR  HIT 

514. 

IF  (IZII)  .EO.  12(1-1)1  GO  TO  26 

515. 

24 

K = K + l 

516. 

I ZONE (K , 1 ) * 1-1 

517. 

I20NE1K.2)  * I-l 

518. 

C 

IF  NO  HITS  IN  ZONE  INCREMENT  ZONE 

519. 

IF  (IZ(I)  .GT.  (K-l))  GC  TO  24 

520. 

I ZONF (K , 1 ) = 1 

521. 

1Z0NFIK.2)  = I 

522. 

GO  TO  28 

523. 

C 

INCREMENT  upper  hit  IN  zoni 

524. 

26 

I ZCNE  (K  , 2 ) = I 

525. 

28 

CONTI NUE 

526. 

30 

CONTINUE 

527. 

IF  ( IKH)  ,GT,  50)  GO  TO  36 

528. 

C 

FILL  ALL  EXCFSS  ZONES 

529. 

K1  = K ♦ 1 

530. 

00  32  L = K 1 , 50 

531. 

I ZONE  (L  , 1 ) = NT  ♦ 1 

532. 

32 

I ZCNE  IL, 2)  = NT  1 

533. 

36 

CONTINUE 

534. 

IF  (KTEST  .LT.  3)  GO  TO  50 

535. 

WRITE  (6,101) 

536. 

DO  40  K = l,  50 

537. 

40 

WRITE  (6,102)  K,  I ZONE (K , 1 ) , I ZONE( K, 

538. 

X TC( I ZONE (K , 2 ) , 2 ) 

539. 

50 

CONTINUE 

540. 

RETURN 

541. 

101 

FORMAT  Cl',  ‘TARGETS  BY  ZONE',/, 

542. 

X 1 ZONE  LOWFR  UPPER  (LOWER  UPPER 

543. 

102 

FORMAT  ('  ',  14,4X , 13 ,4X ,I3,7X,F4.0, 

544. 

END 

545. 

SUBROUTINE  BOMB 

546. 

COMMON  /ARRAYS/  TGT(250,13),  ATT(50,1 

547. 

X I ZONE (50,2),NhIT(250) ,MHIT 120)»HIT(20» 

548. 

X *P( 250, 3) ,COV ( 2 50) , MTYPE ( 10) , NAME (2  50 

ZON?  ANO  I IS  FIRST  HIT 


T0(250,2), 


Hi  ’Sm® 


--  v .Vsr;. 
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549  . 

CGMRON/INT/NT.NA ,N0,  NTM , KT5S7  , NCR , KC.  w,MCL  , MQ9 

•550. 

CORHON  /CBUHI  T / C8UH-<  200,2)  , IRS50),  KCBU,  ! 

551  . 

IF  (KTtST  • GT  • 2)  «T<ITf  (6,102i 

552. 

NCBUHT  = 0 

553. 

KC  BU  =0 

554 . 

DO  40  1=1  ,MA 

555. 

IR  ( I ) = 0 

556. 

NW  = ATT ( 1,9) 

557. 

FMO  = AMOINW, 1,1) 

558. 

IF  ((MOOS  .IT.  0)  .09.  (ATT (1,11)  .60.  1.0)1 

559. 

C 

TF  (MOOI  . EO . 0)  GO  TO  1 

560. 

C 

RN  = RANDT(1.) 

561. 

C 

GO  n 2 

562. 

1 

RN  = P.AN(  1) 

563. 

2 

IF  (S')  . Ifc  . 4TT(  1,11))  GO  TO  3 

564. 

IP  ( I ) = 1 

565. 

GO  T t 40 

566. 

3 

CONTINUE 

567  • 

IF  I6KD  .LT.  0.0)  KCBU  = 1 

566. 

0X=0 

569. 

nv=o 

570. 

NS=ATT (1,7) 

571. 

PHl = ATT ( I ,11/57.3 

572. 

S=$ I N ( PH I ) 

573. 

C=COS(PHI  ) 

574. 

SIGR=:.4B3*ATTf 1,4) 

575. 

SIG0=1.4B3=ATT< 1,5) 

576. 

CALL  GAUSS ( S I GR  ,R  EF  R ) 

577. 

CALL  GAUSSISIOO.OERR) 

57e. 

AGZX=ATT( 1,2) ♦RERR-S»D(PF*C 

579. 

AGZY= ATT  ( I,  3M-RERR*C-Di.RF*S 

580. 

BOGZX=AGZX-S*-ATT(  1,81/2. 

581. 

BDGZY=AGZY-C*ATT( l,8)/2. 

582. 

IF  (NS  .LT.  ?)  GO  TO  10 

583. 

D*ATT  U * 8 )/ « NS-1 ) 

584. 

cx=s*o 

585. 

0Y=C'0 

586. 

10 

CONTI  NUS 

587. 

StGX  = A T T ( ! ,6) 

588. 

SIGY  = ATT (1,10) 

589. 

DO  20  M=1,NS 

59C  . 

IF  (SMC  .LT.  0.0)  NCRUHT  = NCBUHT  + 1 

591. 

IF  (NCB'JHT  .GT.  200)  GO  TC  60 

592. 

IF  ( MODE  .LT.  0)  GO  TC  13 

593. 

C 

IF  (NODE  .EO.  0)  GO  TO  11 

594. 

C 

RN  * RANDT ! 1 . I 

595. 

C 

GO  TO  12 

596. 

11 

RN  = RAN(l) 

597. 

12 

IF  (RN  .GT.  WPNREL(NW))  GO  TO  17 

598. 

13 

CONTINUE 

599. 

CALC  GAUSS (S 1 G X, X) 

600. 

CALL  GAUSS (S IGY, Y I 

601. 

BAGZX=BOGZX«-X 

602. 

BAGZY=BDGZY«-Y 

603. 

IF  ( ( KTFST  .GT.  2)  .OR.  (NPRINT  .LT.  -1)) 

604. 

X WRITE  (6,101)  I,  ",  BAGZX,  BAGZY 

605. 

IF  (5MD  .GE.  0.0)  GO  TO  16 

606. 

C3UHT (NCBUHT , 1 ) = BAGZX 

607. 

CBUHT (NCBUHT ,2)  = BAGZY 

608. 

GO  TO  18 

609. 

16 

IF  ( IBAGZX+BAGZY)  .LT.  -500.)  GO  TO  18 

v}&VatSZ&3*& 


BOMB  #',!<,,«  X-OIM  * ,F8.0, 


INCEX  = I 

CALL  TE  STHT 1 1 NDSX.BA GZX, E AGZ Y ) 

GO  TO  18 

IF  IcMD  .LT.  O.O)  CBUHT (NCBUHT , 1 ) = -10000. 

3DGZX=B0GZX«-DX 

BOGZY  =BOG  ZY4-0Y 

CONTINUE 

CONTINUE 

RETURN 

WR l Tc  (6 , 1 03  ) 

STOP 

FORMAT!'  •,/, 'ATTACK  rf',14,'  BOMB  *',!4,'  x-OIM  *,F8.0, 

X • Y-HM  •,  =8.0) 

FORM-T  ( * 1« ,301 • »• ) , • BOMB  IMPACT  4N0  HIT  DATA  ',  30 (•*•>,/,/ ) 
FORMAT!  'O'*  ' THr  CBUHT  ARRAY  MUST  BE  L-NLAkGED  TO 
X • AC  CO  MO DATE  MORE  CBU  WEAPONS') 

ENO 

SUBRJUTINt  TFSTHT( I ,BX,BV) 

COM«‘N  /ARRAYS/  TGT(?50,13),  ATT  ( 50  , U I . A'OIIO,  20 ,2 1 . TOI250.21. 
X ! ZONE ( 50. 2 I , NH JT  1 256 ) »MHI 7 (20) »H1 T ( 20, 3,25  I »NB  W( 5 I ,HITR t 5, 3, 250 ) 

X »P(250»3|»C0V(250),mtYP£( 10) .NAME (2  SO,?) , WPNREL  1 10) .NCRUI250) 
COMmTN/!NT/NT,NA  ,M!)»  NTM.KTFST,  MCR  , MO  M , MCL  » MOOe  »NPF  INT  » NAM.NST  ,MTT 
XY  = BX  t BY 
NN=XY/500. 

K = MflxOINN'-l.l) 

K = "INOIK.AO) 

U = 0 

if  (K  .to.  n s:;  TO  10 

LL  = IZONE  IK-1,  1) 

IF  (K  .£0.2)  GO  TO  20 

IF  I I IZONEIK-l, 1)  .EC.  I/CNEIK-2, 2)  ) .AM'. 

X ( I ZONE  I K -2 , 2 ) .NE.  l)i  LI  = IZDNEIK.l) 

GO  TO  20 

LL  = IZONE  (1,  1) 

CONTINUE 

LU  = IZONE  (<K*l),  2) 

IF  (KTEST  .GT.  3)  WRIT!  (6,102)  I,  LL , LU 
00  100  IL  = LL,  LU 

CONSIDER  ALL  TARGETS  BETWEEN  THE  LIMITS  OF  LL  AND  LU 
L = TCI IL,?) 

IF  (L  .Lc.  0)  GO  TO  100 

IF  I I TGT  I L,  12  ) TGT  IL  , l 3 ) ) .GT.  500.)  GO  TO  100 
U = 1.A1A*AM0(ATTII  ,0)  .TC-TIL,  10)  , 1) 

IF  KTGTIL.l)  - Cl  .GT.  BX)  GO  TO  100 

IF  I I TGT l L, A)  ♦ 0)  .LT.  BY)  GO  TO  100 

IF  I I TGT I L, 5 ) «■  D ) .LT.  BX)  GO  TO  100 

IF  I I TGT  I L, 8)  - D)  .GT.  BY)  GO  TO  100 

IF  IKTEST  .GT.  '♦)  WRITE  I 6,  101 ) I , L,  BX.BY 
CALL  HITTGTII ,L, BX.BY) 

CONTINUE 

00  120  L = 1,  NT 

IF  1 1 TGT I L, 12 ) ♦ TGT  I L , 1 3 ) ) .LF.  500.)  GO  TO  120 


GO  TO  100 


GO  TO  120 


IF  IKTEST  .GT.  A)  WRITE  (6,1015  I, L, BX.BY 
IL  = L 

CALL  HITTGTII ,IL,BX,BY> 

CONTINUE 

CONTINUE 

RETURN 

FOFMAT  (•  • ,10X,*TESTHIT:  ATTACK  »,I3,‘  TGT  *,IA,  • ) 

X ■ , F6.0,  • Y-OIM  • , T6.0) 

FORMAT!*  »,20X, 'ATT  ACKER  ',13,'  TARGET  RANK  LIMITS  ',216,/) 


-o8- 


SUBROUTIHE  HITTGT(t,L,BX,RY) 

COMMON  '7 RO AYS/  T GT ( 250, ' 3 1 ♦ ATT ( 50 , 1 1 ) , A“D( 10, 20 , 2 ) . T0(250,2), 
XI  ZONE  I 50,2 ),NH IT <2501  .MHITC20) ."I T t 20, 3 , 2 5 ) ,NRW(  5 ) ,HI TR( 5 , 3 » 250 ) 

X *P ( 250, 3 1 fCOV  (250), M TYPE!  10 ) , NAME ( 2 50 , 2 ) , WPNREL { 10) ,NCBU( 250) 
COMMON /I  NT /NT  , NA,N0*N7M,KTFST  » MC?  , MCW, MCL  » MODE  »NPRINT,NAM,NST,MTT 


677. 

01 

MANSION  GNU) 

,F(4)  , 

, I S t A 

■ ),D2( 

4) 

67P . 

= 0 

675. 

T 

= 0 

680. 

DO 

= AMD 

( ATT ( ! , 

9) ,TGT(L». 

.0) ,1) 

681. 

s 

= SIM 

(TGT  (L 

, 9)757.3) 

682. 

C 

= COS 

(TGT  (L 

,9)75 

7.3) 

682. 

CO 

“ s. 

684. 

IF 

( DO  . 

EO.  0.0 

> GO  T0  10 

685. 

T 

= (S  + 

C!  *DP 

686. 

n 

= (S  - 

c)*rr 

687. 

10  XI 

= TOT 

(L  , 1 ) 

- T 

68E  . 

V 1 

= TGT 

( L , 2 ) 

+■  n 

685  . 

X2 

= TGT 

(L,  5) 

f 0 

690. 

Y2 

= TGT 

(L,4) 

+ T 

691  . 

X3 

= TGT 

( L , 5 ) 

* T 

692. 

Y3 

= TGT 

(L  ,6  ) 

- D 

693. 

• X4 

= TGT 

(L,7) 

- D 

694. 

Y4 

= TGT 

(L.fl) 

- T 

695. 

IF 

(tax  . 

LT.  XI) 

.DP  . 

(BX 

.GT. 

X 3 ) ) 

696. 

IF 

((BY  . 

LT.  Y 4 ) 

• DR. 

(PY 

.GT. 

”2 ) ) 

697. 

IF 

(TGT 

(L,  9) 

.£0.  0.) 

GO  TO 

60 

GO  TO  100 
GO  TO  100 


T = S/C 
C = l./T 


700. 

IF 

((BX  .LT. 

X2  ) 

.AND. 

(BY 

.GT. 

( Y 1 *C* ( BX-X1 1 1 ) ) 

GO 

TO 

100 

701. 

IF 

((P.X  .G’. 

X2  ) 

.AND. 

( E Y 

.GT. 

( Y2-TM3X-X2)  ))  ) 

GO 

TO 

100 

702. 

IF 

((BX  .GT. 

X4) 

♦ AND. 

(PY 

.LT. 

(Y4+CMBX-X4)  ) ) ) 

GD 

TO 

100 

703. 

IF 

< ( PX  .LT. 

X4» 

.ANC. 

IFY 

.LT. 

( Yl-TMBX-Xt)  ))  ) 

GO 

TO 

100 

(NPPINT  .IT.  0))  WRITE(6,101H,BX,BY 


(MCR  .ME.  0))  GO  TC  110 


60  CONTINUE 

IF  ( (KTEST  .GT.  2)  .OR. 

NHITd  ) = NHIT(L)  ♦ 1 
NW  = ATT! 1,9) 

NTGT  = TGT(L, 10) 

IF  ((NTGT  .EO.  1)  .AND. 

PP  = AMO ( NW , NTGT , 2 ) 

RESULTS  INCLUDE  AN  EST I “ATE  OF  THAT  FRACTION  OF  THE  TARGET  AREA  THAT 
IS  COVERED  BY  A CIRCLE  OF  RADIUS  »EMD'.  IF  A VALUE,  NOT  GREATER 
THAN  ONE,  IS  INPUT  WITH  THE  SUPPLEMENTARY  EMD  CARD,  A HIT  WILL  BE 
ASSUMED  TO  ACHIEVE  THAT  FRACTIONAL  KILL  OF  THE  TARGET.  IF  A VALUE 
GREATER  THAN  ONE  IS  SPECIFIED,  THAT  VALUE  WILL  BE  USED  IN  THE 
MANNER  AS  THE  EMD  TO  COMPUTE  AN  ESTIMATE  C1F  THE  TARGET  FRACTION 
THAT  IS  COVERED. 

COMPUTE  DISTANCES  NORMAL  TO  THE  FOUR  SIDES  OF  THE  TARGET. 

GO  TO  75 


719. 

IF  (TGT( 

L ,9 ) .EO.  0. 

0) 

720. 

YD  = 

Y1 

+ C*(BX-Xl) 

- BY 

721. 

DN<  1 ) 

= 

S*YD 

722. 

YD  = 

Y 2 

- T*( 6X-X2) 

- BY 

723. 

ON  ( 2 ) 

- 

CO*  YD 

724. 

YD  = 

BY 

- Y4  - C* (B  X 

» 

X 

725. 

DN(  3) 

= 

S*YD 

726. 

YD  = 

BY 

- Y1  + T* (B  X 

- XI 

727. 

ON  ( 4 ) 

= 

CO*  YD 

728. 

GO  TO 

80 

729. 

75 

CONTINUE 

730. 

ON  ( 1 ) 

= 

8X  - XI 

73 1. 

ON  ( 2 ) 

= 

Y2  - PY 

ir-->""  ■ 


5>’jS»i.  --f-  ..».jajfefes»fif« Su>i 


i I 

a tp 


-3 


ON t 3 ) = X3  - GX 

ON  1 4 ) - BY  - Y1 

CONTINUE 

TOT  = 0.0 

DO  8 3 II  = 1,4 

IS(N)  = 0 

02  ( N ) = CD-ON(N) 

IF  ( D 2 ( N ) . LE . 0.0)  GO  TO  83 
TOT  = TOT  + 0?( N1 *02 (N ) 

CONTINUE 

IF  (TOT  .GT.  00*0")  GO  TO  100 


1 

743  . 

RL  = 00 

1 

744. 

IRAO  - 0 

745. 

81 

CONTINUE 

* 

746. 

AL  = 3. 141594  RL  *PL 

1 

747. 

Nt 

II 

2 

>*• 

CO 

O 

C 

c 

748. 

F (N ) = 0.0 

f 

749. 

R = 00  - 0NIN1 

i 

750. 

IF  <(R  .GT.  0.0) 

• ANC.  (R  .GE.  RL)) 

? 

751. 

IF  <F(N)  .EQ.  1.0)  GC  TO  84 

j; 

752. 

IF  ( ( F .LT.  0.0) 

.ANC.  ( -R  .GT . RL ) ) 

f 

753. 

IF  ( R .LT.  0.0) 

R = -R 

\ 

754. 

IF  <R  .EQ.  0.0) 

GO  TO  82 

755. 

Z = R/RL 

? 

756. 

RHO  = 2 . * AT  AN  ( ( 1. 

/( Z*Z )-l  .)*♦(  .5)  ) 

t 

757. 

A = RL*Rl*<RHO  - 

SIN (RHO ) >/2. 

J* 

758. 

IF  (DN(N)  .GT.  00)  A = AL  - A 

£ 

759. 

F ( N ) = 1.  - A/AL 

760. 

GO  TO  84 

i 

761. 

82 

F(N)  = ,5 

762. 

84 

CONTINUE 

i 

763. 

FX  = 1.  - F ( 1 ) - 

FI  3) 

< 

764. 

FY  = 1.  - F ( 2 > - 

F ( 4 ) 

k 

765. 

-IF  (FX  .LT.  0.0) 

FX  = 0.0 

\ 

766. 

IF  (FY  .LT.  0.0) 

FY  = 0.0 

2 

767. 

IF  (KTEST  . GT .4  ) 

WRITE  (6,102)  FX, 

768. 

C 

NOTE  THAT  THE  USE  OF  FX  ANr  FY  PROVIDES 

769. 

FAC  = FX*FY*AL/(TGT(L»12)*TGT(L,13) ) 

VL 

770. 

PS  = 1.  - AMINK1 

.,  FAC) 

771. 

IF  ( IRAO  .EQ.  1) 

GO  TO  91 

772. 

90 

P(  L , 1 1 = 1.  - (l. 

-P ( L» I ) ) *PS 

773, 

IF  (PP  .LE.  1) 

GO  TO  92 

* 

774. 

IRAO  = 1 

775. 

RL  = PP 

z 

776. 

GO  TO  81 

1 

/ 

777. 

91 

P(  L , 3 ) * 1.  - (l. 

- P(L,3))*PS 

w*1 

778. 

GO  TO  93 

1 

779. 

92 

IF  (PP  .EQ.  0.0) 

GO  TO  93 

Yn 

V 

780. 

P( L , 3 ) = 1.  - (1. 

-P(L,3) )*<1.  - PP) 

F ( N ) = 


GO  TO  84 


CONTINUE 

IF  (KTEST  .GT.  31  WRITE 
IF  ( T GT ( L » 1 1 J .LT.  I.)  G 
IF  (NHIT(L)  .GT.  25}  GO 
00  95  J = It  NST 
IF  (MHIT(J)  .EQ.  0)  GO  TO 
IF  (MHIT(J)  .NE.  LI  GO  T 
HIT  (J,1,NHIT(L))=BX 
HIT  CJ,2,NHIT(L))=BY 
HIT  (J, 3,NHIT(L1)=ATT(I,9) 
GO  TO  100 
CONTINUE 


WRITE  (6,10011 
) GO  TO  100 
GO  TO  100 

GO  TO  100 
GO  TO  95 


P (L,  1 1 » P(L , 3 1 


ij 


i‘*  ',''~'  J;''  ,jfal^;|^  ■if  ) '' 


“'>*  ■“".  -T’*t 


^ ■-». p'*-k ^-v ^ , j.  ^ w. 


ag*S3ay£3>vc«gr?^^  g> 


»,  1SX.  •**▼*«  HHTGT 
•,  • FX  • *F6.3  , • 
•.*  TGT  H ',13,* 


NMIT  ( • , J3  , • ) * , JP 6. 0) 

FY  *,F6.3> 

PK1  ' * F8. 3 1 • PK3  • ,F8. 3) 


CONTINUE 

IF  (NTGT  .NE.  1)  RETURN 
CONT I NU5 

00  120  J * 1,5 
IF  (NRW(J)  .tO.  0)  GO  TO  130 
IF  (NRW(J)  .Nr.  L)  GO  TO  120 
HIT»( J,1,NH!T(I  )!  * OX 
H!TR(J,2,NHIT(L))  = RY 
HlTR(  J,3,NH!T(U)  = 4TT(l,g) 

CONTINUE 

CONTINUE 

RETURN 

F0FM4TP  »,  IPX,  «**▼*«  HITTGT  NHIT  ( • , 13  , • ) * , ?P  6. 0) 

FORMATC  ',  ' FX  ' *F6.3 1 • FY  ',F6.3> 

L FORMAT ( • •,«  TGT  H ',13,*  PK1  *,F8.3,*  PK3  *,F8.3) 

FNO 

SUBROUTINE  GAUSS  ( S ♦ V> 

COMMON/) NT /NT  ,NA  ,ND,NTM,KTEST,MCR,MCW,MCL,MOOfc  ,NPR  INT , NAM ,NST ,MT  T 
IF  IIKTEST  .GT.  71  . TR.  {MODE  .LT.  0)1  GO  TO  50 
A=0.0 

if  (mcds  .so.  o)  g:  TO  20 

DO  10  1=1,12 

V = RAN0TI1.I 
A = 4 f Y 

GO  TO  40 
00  30  1=1,12 
Y=RAN(I) 

A = 4 + Y 

V = ( A-6 . 01 *5 
RETURN 

CONT I NUC 

V = 0.0 
RETURN 
END 

SUBROUTINE  CHECKR 

COMMON  /ARRAYS/  TGTI  250.13),  ATT ( 50 , 11 ) , AMD( 10, 20 ,2 1 , T0!250,2). 
XI  ZONE  I 50, 21, NHIT 1 250) ,MHI T (20 ) ,Hl T (2 0, 3 , 25 1 , NRW( 5 ) ,HITR( 5, 3, 250 > 

X »P(250,3) ,COV 1 2 50 ) , M TYPE ( 10) , NAME ( 2 50, 2) , WPNREL ( 10) ,NCBU(250) 
CCMMON/INT/NT,NA,NO,NTM,KTEST,MCR,MCW,MCL,MODE,NPRINT,NAM,NST,MTT 
COMMON/STATS/NTRIAL, ITR I A t , NST AT , ST AT  I 250 , «) ,STAT2(20,5) , MST  AT ( 8) 
COMMON  /CONTRl/  NREP.NRLCT, INW, INL , NS AVE1 , NSAVE 2, LI  ST ,N JMEM 
NC  = 0 
NN=0 

LHGLFS  = 10000 
00  40  MRW  =1,5 

CYCLE  THRU  AS  MANY  AS  8 RUNWAY/TAX  I WAYS . 

IRW  = NRW(MRW) 

EXIT  IF  NO  TARGET  NUM35R  (IRW)  FOUND. 

IF  (IRW  .EQ.  0)  GO  TO  50 
NN=NN+1 

IF  (NHIT(IPW)  .EO.  0)  GO  TO  40 
IF  (KTEST  .GT.  4)  WRITE  (6,102)  IRW 
INDEX  = MRW 

CALL  RUNWAY  (INDEX,  IRW,  ICONO,  NHOLESI 
RUNWAY  SU8P0UTI NE  RETURNS  ICONO  = 0 IF  RUNWAY  HAS  PEQUIRED  SPACE; 
ICONO  = 1 IF  NOT. 


IF  (ICONO  .EO.  1)  NC=NC*1 

IF  (NHOLES  .LT,  LHOLES ) LHOLES  = NHOLES 

IF  (NPRINT  .EO.  4)  WRIT6< 6, 103 ) ITRI AL , I RW, NHI T ( IR W) , NHOLES 

CQNTINUC 


V- 


wmsm 


,« ..-,  -.  KW.K'iBitswaE-irS'yi 


IF  (NC  .EO.NN)  GO 
IF  (NPRINT  .LT.  3 ) 
NSTAT  = NSTAT  ♦ 1 
CONTINUE 

IF  ( (NFEF  .SO.  1)  . 
RETURN 

MSTAT(l)  = MSTAT(l) 
WST AT  (2  ) s MSTUI2) 


TO  60 

WRITE  (6,101) 


MSTAT(l)  = MSTAT(l) 
mSTAT(2)  s WST ' T( 2) 
RETUrN 

FORMAT!  • • AT  LE, 

FORMAT ( • CHECK 

FORMAT! • •.•TRIAL*, 

2ND 

SUBROUTINE  RUNWAY  ( 
COMMON  /ARRAYS/  TGT 
X I ZONE ( 50 ,2 ) « NH  l T (260 
X ,P(250,3) ,C0V(250), 
COMMON/ ST  ATS /NT R i AL 
common/i mt/nt,na,no, 

COMMON  /HITS/  XN( 2 
COMMON  /C0NTRL/  NRi 
DIMENSION  NTS  ST ( 250 


AND.  (NC  .80.  NN>)  GO  TO  70 

♦ I.HCLES 

♦ LHOL2S*LH0LES 

ST  CNE  RUNWAY  IS  AVAUA8LE') 
TARGET  *•  , !'>) 

14,*  TGT • » 14  > • HITS', 14,* 


REPAIRS', 13) 


MRW,  IRW,  I COMO,  NHPLCS) 

(250,13),  ATT ( 50 , 1 1 ) , AMD ( iO, 20,2 ) , T0(250,2), 

) ,MHIT(20) ,HIT(20,  3,25),NFW( 5)  ,HI TR( 5 , 3 , 250) 
MTY°E(10), NAME (2 50, 2) ,WPNREL( 1C) »NCBU ( 250 ) 

, ITPIAL,NSTAT,STAT(250,8) ,STAT2(20,5) ,MSTAT( 8) 
NTM.KTEST  , MCR.MCW.MCI  . MODI:  ,NPR  INT  , NAM  ,NST  ,MTT 
50) ,YN(250> ,NZt?50) 

P.NPLl T, INW,INL,NSAVcl,NSAVc2,LIST,NJMEM 
),  YH ( 2H0 ,2  I 


CHECKS  FOR  THE  EXIS^ANCE  CP  THE  SPECIFIED  RUNWAY  MIMMUMS  (MCL  X 
MCw)  ON  EACH  RUNWAY  AND  DESIGNATED  TAXIWAY  (TYPE  #1  TARGETS), 
STOPS  SEARCHING  A GIVEN  RUNWAY  WHENEVER  PEQUlRt.MENT  IS  SaTIS.  IED. 


TH=TGT(IPW,9)/S7.3 

ESTABLISH  ORIGIN  (XO,  YO)  FOR  A fECTANGULAR  COORDINATE  SYSTEM  WITH 
THE  X-AXIS  ON  THF  MOPE  SOUTHERLY  FOGS  OF  THE  RUNWAY. 

NH I =0 
NHOLES  =1000 
00  5 N = l,  250 
NZ(N)  = 0 

IF  ( TGT  ( IRW,  12  ) .C,T.  TGT  ( I RW,  13  ) ) GO  TO  10 
NO  IR  = l 

XO=TGT ( I PW,  1 ) 

YO=TGTl IRW,  2) 

LTH  = TGT (I RW , 1 3 ) 

WID  = TGT (IRW, 12) 

GO  TO  20 
NDIR=2 

XO=TGT( IRW, 7) 

YC=TGT( I R W, 8 ) 

LTH  = TGT I I RW , 1 2 ) 

WID  = TGT (IRW, 13) 

CONTINUE 

IF  (KTEST  .GT.  4)  WRITE  < 6, 1004) I RW ,X0, YO ,LTH ,W I D.MCL ,MCW 
IF  (MCW  .GT.  WID)  GO  TO  320 
NHITl  = NHIT(IRW) 

DO  50  1=1,  NHITl 

IF  (AMD(HITRIMRW,3, I ), 1; 1 ) .LT.  0.)  GO  TO  50 
NTW  = HITRIMRW, 3,1) 

GO  TO  60 
CONTINUE 
CONTINUE 
NON  = 1 

EMD  = AMD (NTH, 1,1) 

DO  70  I = l,  NHITl 


WRITE  (6,  10O4)IRW,XO,YO,LTH,WID,MCL,MCW 
GO  TO  320 


GO  TO  50 
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915. 

if  tnmiMPw.3, n 

• EO.  NT W ) GO  T9  70 

916. 

IF  (AMC(HITP(MRW,3 

.1  ),l,l)  .LT.  0.)  GO  Tfl  70 

917. 

NON  = 0 

918. 

EMO  = 0.0 

919. 

GO  T3  so 

920. 

70 

CONTINUE 

921. 

80 

continue 

922. 

00  140  I = 1,  'JHIT1 

923. 

C 

TRANSFORM  wjt  C '‘HR.)  I NATES  10  RUNWAY  C 30R0IN  AT  ES  . 

924. 

NHI  = NHI  1 

925. 

XB  = HITMvp.w.l.n 

926. 

YB  = HI  TP (MRW.2,1 ) 

927. 

IF  (AMDIHITRIMRV, } 

, I . ,1,1)  .LT.  0.)  N7II)  = 1 

928. 

IF  ( TH  .EQ.  0.0  > 

GO  TC  1 10 

929. 

XX=X9-X0 

930. 

YY= Y8-Y0 

931. 

R = ( XX*XX«-YY*YYI* 

-10.5) 

932. 

IF  ( KTt  ST  7) 

WRITE  (6,1010)  R , XX  , Y Y 

933. 

YZ  = YY/= 

934  • 

TH>  = AT  AN ( YY/ XX  I 

935. 

IF  (XX  .LT.  0.0) 

THl  = TH1  + ?.  1 <,1 4 

936. 

TH2  ~ TH ) + TH 

937. 

le  (NO!*  .EQ.  2) 

TH2  = TH2  - 1.5706 

938. 

XNiI)=P«C0S(TH2) 

939. 

YN( I ) =R  *S IN ( TW2 ) 

940. 

IF  (KTeST  .GT.  5) 

WRI Tc 16, 1009)1, XNI I ) , YNII) 

941. 

oG  T.3  130 

942. 

110 

IF  ( NCI  R .FO.  2) 

GO  TC  120 

943. 

XNII)=X3-X0 

944. 

YNI I ) =YB- Yn 

945. 

GO  TO  130 

946. 

120 

XNI I ) = Y3-Y0 

947. 

YNI I I =X0-X3 

949. 

130 

IF  ( I .GT.  249) 

GJ  TO  150 

949. 

140 

CONTINUE 

950. 

GO  TJ  160 

951. 

150 

WRIT*  (6,1001)  IRW,  1TR1AL 

952. 

160 

CONTINUE 

953. 

NH  = NHI 

954. 

C 

IF  NPL9T  .SO.  2 RUNWAY  IMPACTS  A°E  PLOTTED  FOR  ALL  CONDITIONS 

955. 

C 

IF  N2L°T  .EC.  1 IMPACTS  CNLY  PLOTTED  WHEN  RUNWAY  IS  CLOSED. 

956. 

IF  (IMRW  .1:0.  1)  .AND.  INPLOT  ,tj.  2!)  WRITE  1 6,  10 12 ) I TR  I AL 

957. 

IF  INPLOT  .EO.  2) 

CALL  PL0THT(NH,IRW,LTH,wIO) 

958. 

IF  IKTFST  . GT  .6  ) 

WRITE  (6,1006)  NHI T 1 1 P W) ,NH 

959. 

00  170  I — 1 » NH 

960. 

YHI I ,1)  = YNI  I) 

961. 

170 

YH  ( I , 2)  = I 

962. 

IF  INH  .EO.  1)  GO  TO  190 

963. 

DO  180  J = 2» NH 

964. 

00  180  K=2,NH 

965. 

C 

ORDER  ALL  'MTS  FROM 

LOWEST  Y TO  HIGHEST.  YHII.l)  IS  THE 

966. 

C 

Y COORDINATE,  YHI 1,2)  THE  HIT  NUP6ER,  OF  THE  I TH  ORDERED  HIT 

967. 

I=NH-K+  2 

968. 

IF  I YHj l , 1 ) .GF.  YH (1-1,1))  GC  TO  180 

969. 

T*YH( 1-1,1) 

970. 

TN=YHI 1-1,2) 

971. 

YHII-1,1)=YHI 1,1) 

972. 

YHI T -1,2)=YHI 1,2) 

973. 

YH  1 1 , l ) =T 

974. 

YHI  I , 2)  -TN 

>75. 

180 

CONTINUE 

976. 

977. 
976. 

979. 

980. 

981. 

982. 

983. 
989. 
985. 
988. 

987. 

988. 

989. 

990. 

991. 

992. 

993. 
99*. 
993. 

996. 

997. 

998. 

999. 
1000. 
1001. 
1002. 
1003. 
100*. 
1005. 
lOOf  . 

1007. 

1008. 

1009. 

1010. 
1011. 
1012. 
1013. 
101*. 

1015. 

1016. 

1017. 

1018. 

1019. 

1020. 
1021. 
1022. 
1023. 
102*. 

1025. 

1026, 

1027. 

1028. 

1029. 

1030. 

1031. 

1032. 

1033. 
103*. 

1035. 

1036. 


CONTI  MU c 
XL  = 0.0 
XU  = MCL 

a = 0.0 

YU  = 

NYU  = l 

IF  (NON  .EO.  01  0.0  TO  210 

YL  = YL  - SMC 
YU  - YU  + CM' 

CONTINUE 
IHOL  F = 0 

00  250  I - NYt  , NH 

IF  C *-■  2 € I ) .SO.  1)  GO  TC'  250 
YT  = YH(I ,1) 

IF  (YON  . EO.  1)  GO  TC  230 
9 = AM0(HITR(MRW,3,YH(I,2n,l,l> 
YL  = YL  - R 
YU  = YU  + R 

IF  iVT  .LT.  YLI  GO  TO  2*0 
XT  = XN(YH< 1,21) 

IF  ((XT  .LT.  XU  .0«.  I XT  , GT . XI 
IF  (YT  .GT.  YU)  GO  TC  2fcO 

1 HOLE  = .HOLS  * L 

IF  (KTcST  .GT.  6)  WRIT*.  (6*  1008 ' 

NTE  ST ( I HOLS ) = I 

IF  (YFE-  .SO.  0)  GO  TO  260 

IF  (NCN  .F9.  II  GC  TO  250 

IF  (I  .SO.  NH)  GO  TO  260 

YL  = Yl  + o 

YU  - YU  - 9 

CONTINUS 

Conti cue 

IF  (IHCL/.  .SO.  0)  GO  TC  300 
IF  ( NCN  .SO.  1)  NYL  * NTcST(l) 
NHCLCS  = MINO(NHOLSS,  IHCLc) 

IP  (NON  .SO.  01  GO  TC  270 
YL  = YL  ♦ INW 
YU  = YU  IN-, 

IF  IYU  .GT.  (MIOtEM-UJ  GO  TO  2! 
GO  TO  220 


NH 

GO  TC'  250 
>0  TC  230 


GO  TO  2*0 

.0«.  (XT  , GT . XU)) 
GO  TC  2fc0 

WRIT-.  (6. 1008) 

GT  TO  260 
GC  TO  250 
30  TO  260 


GO  TO  2*0 


IHOLE,  YL.YU 


GO  TO  280 


YL  = YL  ' R ♦ INW 

YU  = YU  - R + INW 

IF  (YU  .GT.  WtO)  GO  TO  280 
IF  IKTSST  . GT.  *)  WPITE(6, 
GO  TO  220 
CONTINUE 
XL  = XL  ♦ INL 

XU  = YU  6 INL 

IF  (XU  .GT.  LTh I GC  TO  290 
GO  TO  200 

IF  INPRINT  .GE.  3)  GC  TO  29 
IF  ( ( MRW  .SQ.  1)  .ANO.  (NPLQT 
WRITE! 6, 1002 ) IRW 
IF  (NRFP  .GT.  0)  WRITE  (6,1 
IF  (NPL7T  .EO.  I)  CALI  PLOT 
iCONO  = 1 

IF  HNPLOT  .GT.  0)  .4N0.  (NPR 
RETURN 

IF  (NPRINT  .GT.  2)  GO  10  31 
IF  ( ( MRW  .EQ.  1 ) .ANO.tNPLOT  . 
WRITS  (6,1003)  IRW 


GO  TO  280 

WP1TE(6,1008)  NHHLES 


GC  TO  290 


GC  TO  295 


.NE.  2))  WRITF  ( 6 1 10 12 ) 1 TRI AL 


WRITE  (6,1011)  N HOLES 
CALL  PLOTHT(NH, IRW,LTH,WIDJ 

. ANO.  (NPRINT  .LT.  3)1  WRI TE ( 6 , 1000) 


GO  10  310 
(NPLOT  .NE.  2)) 


WRITE  16,1012)  I TRIAL 


1037. 

1033. 

1039. 

1040. 

1041. 

1042. 

1043. 

1044. 

1045. 

1046. 

1047. 

1048. 

1049. 

1050. 

105 1 . 

1052. 

1053. 

1054. 

1055. 

1056. 

1057. 

1058. 

1059. 

1060. 
1061. 
1062. 

1063. 

1064. 
1C6  5 . 
1066. 

1067. 

1068. 

1069. 

1070. 

1071. 

1072. 

1073. 
107*. 

1075. 

1076. 
1G77  . 

1078. 

1079. 

1080. 
1061. 
1082. 

1083. 

1084. 

1085. 

1086. 
1087. 
10S8 . 
1089. 
109C. 

1091. 

1092. 

1093. 

1094. 

1095. 

1096. 

1097. 
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310 


320 

1000 

1001 

x 

1002 

1003 

1004 

1005 

X 

1006 
1068 

1009 

1010 
1011 

X 

1 0i2 


I CONI  = 0 

IF  KNPLOT  .50.  2)  .AND.  (NPkINT  .LT.  3))  WR I TE ( 6 , 1000) 

NHGLFS  - 0 

R6TURN 

WR!  TE  ((• , 1 005  ) I»W 
STOP 

FORMAT! ' 1‘ ) 

FOKMAT ( 'O' ,'ONLY  FIRST 
• IN  TRIAL  <M,  14,/,/) 

FORMAT! <0. , ,RUNWtY  *•  ,13*' 

FPSWAT  ( <0*  ,' RUNWAY  ''•  * I 3 » ' 

FGR  M AT ( ' • * « RUNWAY  SPl CS 
FORMAT  {•0*, •*«***  TAP  Gi  T 
•FLIGHT  OPERATIONS’ ) 


:50  HITS  TESTED  FOR  TARGET  #',!4, 


IS  CLOSED',/) 

IS  OPEN',/) 

',  I4,4X,2F8.0,I8,Ftt.0,2I8! 
13, 1 IS  TOO  NARROW  FOR  ', 


u HI TS  TL  CHECK' , 216) 

TEST  POINT  #B  • ,21 10.2F10.0) 

TEST  POINT  #E',  14.2^10.0) 

R = *,F8.0,10X,2F10.0) 

HOLES  RUST  HF  REPAIRED  TO  NEE  T RUNWAY 


FORMAT (•  ', 

FORMAT!'  ', 

FORMAT!'  ', 

FORMAT!'  ',' 

FORMAT! • '.  14 
• MlNIMijMS'  ,/  ) 

FCRUAT  ! * 1 ' , • «*«**  TRIAL  «',I3,'  *»*«*',/) 

cNO 

SU8RSUTINE  PR INT 
INTEGER  *4  NAME 

COMMON  /ARP  AYS/  TGT  < 250, 13),  A" T( 50, 11), AMD (10, 20, 2),  T0(250,2), 
XIZ3NE!50,2),NHIT(250) ,MHIT!20) ,Ht T ( 20, 3 ,25 ) ,NRW( 5) ,HITR<  5,3,250? 

X ,P(250,3) ,GOV (?50),MTYPE( 10 ) .NAME  1250, 21  , WPNRcL ( 10) ,NC8U( 250) 
COMMOK/INT/NT,NA,NL',NTH,KTEST,MCR,NCW,MCL,MODE,NPR1NT,NAM,NST,MTT 
COHMON/STATS/NTRIAL, I TRIAL ,NSTAT,STAT(250,3I,STAT2(20,5) ,mSTAT(3> 
IF  (NTPIAL  .EG.  n GO  TO  1 


WRITS  (6,108) 
50  TO  3 
CONTINUE 
WRITE  (6,106) 
CONTINUE 
WRITS  (6,101) 
DO  10  m*1 ,MTT 
NN=  0 

DO  10  1=1, NT 


I TRIAL 


10 


T C 

4 » 

! TGT  C 

1,10)  . 

NE.  M 

) GC  TO  10 

NN= 

NN»  1 

IF 

(NN  . 

EO.  1 ) 

WR  I 

TE  (6,109)  S 

NA 10  = NHIT ( I ) 

- NC8U(  I ) 

WRI 

Tc  (6 » 

102)  I, 

NAI  0 

,COV ( I > ,P( I » 1 

NAME!  I 

,1),  NA*E(I, 

2) 

CONTINUE 

IF 

( NPRINT  .GT. 

0) 

GO  TC  30 

WRITE  (6 

, 103) 

DO 

20 

M - i , 

NST 

IF 

( MHIT 

(M)  .EO 

. 0) 

GO  TO  30 

NN= 

0 

NL 

= NHI 

T (MHI T( 

M)  ) 

IF 

(NL  . 

EO.  0) 

GO 

TC  20 

00 

15 

N=l,  NL 

NN= 

NN+1 

IF 

(NN  . 

SO.  26) 

GO 

TO  20 

IF 

(NN  . 

EO.  15 

WRITE 

(6,104)  MHI 

X=HIT(N, 

1 »NN) 

Y*H 

IITIH, 

2 ,NN) 

NWPN=HJT(M,3,NN) 

WRITE  (6,105)  X.Y.NWPN 


lose. 

15 

COMINUS 

1099. 

20 

CONTINUE 

1100. 

30 

CONTINUE 

1101. 

00  50  « = 1,  5 

1102. 

IF  (NFW(M)  .dQ.  0>  GO  TO  60 

1103. 

NN  = 0 

1104. 

NL  = NH IT (NRW  tM 1 1 

1105. 

IF  (NL  . f.  Q.  01  GO  TO  50 

1106. 

DO  40  N = 1,  Nl 

1107. 

MN  = NN  ♦ 1 

1108. 

IF  (NN  .£Q.  251)  GO  TO  50 

1109. 

IF  (NN  .60.  1)  WHITE  (6.104)  NRW(M) 

1110. 

X = MITR(M.I.NN) 

1111. 

Y = HITR(M  ,2,M) 

1112. 

NWPN  = H!TR(M,3,NN) 

1113. 

WRITE  (6,105)  X,  Y,  NWFN 

1114. 

40 

CONTINUE 

1115. 

50 

CONTINUE 

1116. 

60 

CONTINUE 

1117. 

RETURN 

1118. 

101 

FOPMATI ' 0 • * 

1119. 

X 

•TGT  NO.  C5‘J  ',10X, 

1120. 

X 

• 80MBS  CBU  BLDG • « / * IX , 

1121. 

X 

•NO.  HITS  COVERAGE  ’.10X, 

1122. 

X 

•EMO  OTHER  P K NO.') 

1123. 

102 

FORMS T t*  • ,I3,3X,I4,4X,F6.2,10X,2(3X,F5.3),6X,F5.3,7X,2A4) 

1124. 

103 

FORMAT! 'O' ,1SX, 'HIT  LOCATION  AND  WPN  TYPE  Fr)R  SELECTED  TARGETS* 

1125. 

X 

./,/) 

1126. 

104 

FGRMAT ( • *,•  TARGET  NUMBER ',14, • X-DIM  Y-DIM  *, 

1127. 

X 

•WPN  TYPE*,/) 

1128. 

105 

FORMAT ( • 21X.2F9.0.I 7) 

1129. 

106 

FORMAT!' 1' ,?0X, 'TARGET  HIT  SUMMARY',/) 

1130. 

108 

FORMAT! • 1* ,10X, 'TARGET  HIT  SUMMARY  TRIAL'. 15) 

1131. 

109 

FORMAT ('O',  13X,'«»  TARGET  TYPE  * ',13,'  «*',/) 

1132. 

END 

1133. 

SUBROUTINE  ST  AT  IS 

1134. 

COMMON  /ARRAYS/  TGT(250,13),  ATT I 50 , 11) , AMD( 10, 20 ,2 ) , T0(250,2), 

1135. 

XIZQNE(50,2),NHlT{250)  ,MHl  T ( 20 ) ,Hl  T (20, 3 ,25  1 ,NR  W(  5 1 ,H1  TM  5 .3 , 2501 

1136. 

X 

, P 1 250, 3) ,C0V(2  50),MTYPE( 10) ,NAM£ (2 50, 2) , WPNP tL 1 1C) ,NCBU( 250) 

1137. 

COMMON/!NT/NT,NA,NO,NTM,KTEST,MCR.MCW,MCL,MODE,NPRINT,NAM,NST,MTT 

1138. 

COMMON/ ST ATS/NTRIAL, ITR I AL. NST AT , ST  AT ( 250, 8) , ST AT2( 20,5) , MSTAT! 8) 

1139. 

COMMON  /CONTRL/  NR ED,NP1 OT, INW, I NL ,NSAV£ 1 ,NSAVE2,U ST.NJMEM 

1140. 

WRITE  (6,101)  NTRIAL 

1141  . 

AVGR6P  =0.0 

1142. 

DO  12  M = 1,  ITT 

1143. 

NN  = 0 

1144. 

SUM1  = 0.0 

1145. 

SUM 2 = 0.0 

1146. 

DO  10  I = l,  NT 

1147. 

IF  (TGT(l.lO)  .NE.  M)  GC  TO  10 

1148. 

CONTINUE 

1149. 

NN  = NN  1 

1150. 

IF  (NN  .EO.  1)  WRITE  (6,102)  M 

1151. 

FHIT  = (STATU,  3 ) /NTRIAL  > *100. 

1152. 

AH I T S = STATI 1,1) /NTRIAL 

1153. 

SUM!  * SUM1  «■  AHITS 

1154. 

TRIAL  = NTRIAL 

1155. 

SOH  = (STATU, 2)  - TRI AL* AHITS*AHIT $) / C TR I AL  - 1.) 

1156. 

SDH  = SDH**(0.5I 

1157. 

ACOV  = ST  AT  (1,4) /TRIAL 

1156. 

SUM2  * SUM?  + ACOV 

1159. 

1160. 
1161. 
1162. 

1163. 

1164. 

1165. 

1166. 
1167. 
1166. 

1169. 

1170. 

1171. 
1172  « 

1173. 

1174. 

1175. 

1176. 

1177. 
1176. 

1179. 

1180. 
1181. 
1182. 
1183. 
1164. 

1185. 

1186. 

1187. 

1188. 

1189. 

1190. 

1191. 

1192. 

1193. 

1194. 

1195. 

1196. 

1197. 
1196. 

1199. 

1200. 
1201. 
1202. 

1203. 

1204. 

1205. 

1206. 
1207. 
1206. 

1209. 

1210. 
1211. 
1212. 

1213. 

1214. 
1215* 
1216. 

1217. 

1218. 
1219. 


CO  TO  20 


SDC  = (STM*  I , S*  - TR I A l *4C0V»AC0V ) / ( TR  I AL  - i.l 

SDC  - <0O*  10.51 

PK1  * STATU  ,6)/TRt AL 

PK2  * ST  AT { 1 *71 /TP  1 Al 

PK3  * STATU  ,R)/TR  ! AL 

WRITc  (6.10’)  l,f-HIT,AHIlS.Sr>H,AC0V,SDC.PKl(PK3,PK2,NAMEU,n. 

X N A M F ( 1 , 2 • 

COM  1 Nil? 

IF  INK  .GT.  01  WRITE  (6,108)  SUMi , SUM2 

CONTINUE 

WRITE  (4,106) 

00  25  M x i,MTT 
00  15  I * 1 ,NT 

IF  ( TGT(I.IO)  .FO.  M)  CO  TO  20 

CONTINUE 

GO  TO  25 

FHIT  = STAT2(M, 11/TftIAL 

SDH  = ( STAT  2*  M , 2 > - TRIAt  *FH  I T*F  H I T ) / (TR1  AL  - 1.) 

FCCV3  = 100.*(STaT2(M,3)/TRIaL) 

FCCV4  = IOO.*(STAT2(M,4)/TMAL) 

FCGV5  = 1Q0.*(STAT2(M,5)/TRIAL) 

FHIT  = 100. *FHI T 

SDH  = n0.*(S0H**(0.5)  ) 

WRIT-  (3,107)  M,  FHIT,  SDH,  FC0V3 , FC0V5,  FCCV4 

CONTINU 

CONTINUE 

IF  (MCF  .£0.  0)  GO  TC  50 
STA  = NSTAT 

FOPEN  = ( S‘rA/TRiAL)*100. 

NClSO  = MTRIAL  - STA 

IF  ( NCLSO  . £0 . 0)  GO  TC  40 

CLSD  a NC ISO 

AVGRFP  * WSTATf  D/CLSD 

SDREP  - 0.0 

IF  (NCISC  .GT.  1) 

X SOPFP  = ( ( MSTAT ( 2) -CL  SD» AVGREP*A VGREP ) /(CLSD-l. ))»*(. 5) 

CONTINUE  - 

WRITfc  (6,104)  FOPEN 

IF  ( NRt 0 .EO.  1)  WRITE  (6,105)  AVGREP,  SUREP 

CONTINUE 

RETURN 

FORMAT  ( • 1* , 10X, 'TARGET  DAMAGE  STATISTICS  FOP', 14,'  TRIALS',/,/, 
X • TARGET  PERCENT  AVERAGE  HITS  STD.  DEV.  AVG.  C8U  STD.* 
X,'  DEV.  AVG.  8CMB  CCVER/Gc  C8U  BLDG',/, 

X • NUMBFk  ATTACKS  HIT  PER  ATTACK  OF  HITS  COVERAGE  COVE', 

X • RAGE  EMD  OTHER  PK  NO.',/,/) 

FORMAT  (•  ',/,  15X,  'TARGET  TYPE  * ',13 ,/) 

FORMATU  ', I6,6X,F6. I, 6X, F7 .2, 6X, F6 .2 »4X , F6. 2 ,4X , F3. 2 , 

X 2F9.3,3X,F9.3,2X,2A4) 

FORMATU  ',/,  • AT  LEAST  ONE  MINIMUM  RUNWAY  SECTION  WAS  OPEN  AFTE 
XR',  F6. I , ' PERCENT  OF  THE  ATTACKS',/) 

FORMAT!  ' ' , ' WHEN  ALL  RUNWAYS  WERE  CLOSED,  • , FA,,  1 , • ( « , F4. 1 , 

X')  HOLES  REQUIRED  REPAIR,  ON  THE  AVERAGE,  TO  PROVIDE', 
y • A MINIMUM  RUNWAY',/) 

FORMATU  «,/,/,lOX,' CAMAGE  STATISTICS  §Y  TARGET  TYPE*,/,/, 

X • AVERAGE',/, 

X • TARGET  PERCENT  STANDARD  COVERAGE  ',/, 

X • TYPE  HIT  DEVIATION  EMD  OTHER  CBU'/) 

FORMAT ( * • » 5X « 12,  2(6X,F5.1),  3(2X,F5.1)l 

FORMATU  • ,25X  » • » ,16X,' • ,/ ,26X  ,F6.2. 16X,  F6.2  ) 

END 
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/ 
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1^-  ^.flfeft- 

M 

It 


if 


1 1 
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i i 


n I 


a®"t  34 

as  m 


1220  . 
1221. 
1222. 
1223» 

1224. 

1225. 

1226. 

1227. 

1228. 

1229. 

1230. 

1231. 

1232. 

1233. 

1234. 

1235. 

1236. 

1237. 

1238. 

1239. 

1240. 

1241. 

1242. 

1243. 

1244. 

1245. 

1246. 

1247. 

1248. 

1249. 

1250. 

1251. 

1252. 

1253. 

1254. 

1255. 

1256. 

1257. 

1258. 

1259. 

1260. 
1261. 
1262. 

1263. 

1264. 

1265. 

1266. 
1267. 
1263. 

1269. 

1270. 

1271. 

1272. 

1273. 

1274. 

1275. 

1276. 

1277. 

1278. 

1279. 

1280. 


SUBROUTINE  Cfl'J 

COMMON  /ARRAYS/  TGT(250,13),  ATT ( 50 » 1 1 ) » Av  .,20,2),  T0<250,2), 

X 1 ZONE  ( 50 , 2 } , NH  IT  ( 250)  ,MHIT(20>  ,Hl  T (20,3 ,25  ) ,NRW(  5 ) ,Hl  T<U  5 , 3, 250 ) 

X » P(  250 , 3)  ,COV  ( ? 50) , M TYPE  ! 10 ) , NAME  ( 2 50 , 2 ) , HPNF  EL  ( 1C)  ,NCBU(  250) 
COMMON /!  NT /NT,  NA  ,ND,  NTM , KT  EST  , KCR  , MC  N,MCl. , MODE  ,NPR  INT  , NAM , NST  , MT  T 
COMMON/S i 4TS/NTRI4L, I TR l AL , NST AT, ST ATI  250 , 8 ) , STA T2 (20 , 5) ,MSTAT(B> 
COMMON  /CBUHIT/  C8UHT ( 200,2 ) , IR<50),  KCBU,  KPT  1 
OIMC-NSI  IN  NCOV(  16,161,  ICrVC16, 16)  ,PSCDV(  16,16)  , 

X X ( 4 ) , Y(4),  XX(4|,  VY(4),  WC(10) 

COMPUTE  PATTERN  DIAGONAL  DIMENSION 
00  2 NWPN  = 1,  10 
WD(NWPN)  = 0 
LTH  = - AMD ( Nw  PN , 1 , 1 ) 

IF  (LTH  . LE . 0)  GO  TO  2 
WID  = AM0 (NWPN, 2*1) 

WD(NWPN)  = ( t LTH»LT  H ♦ W I D* W 1 0 ) ** ( . 5 ) ) /2 . 

CONTINUE 

AT  THIS  POINT  THE  PROGRAM  NOW  CHECKS,  TARGET  BY  TARGET , FOR 
WHATEVER  C3U  * S «AY  HAVE  CCV1PED  AMY  P4RT  OF  THE  TARGET. 

'IC  200  L = 1 , NT 

«»*  ^ ,*#*:.#* ***«  «<*  c POP  c ACH  TARGET 

INIT1  = 0 

COMPUTE  T A° Gc T CcNT£k 

TCX=(TGT(L,l)«-T'jT(L*5)  ) *0.5 
TC  Y=  ( TGT  ( L , 2 ) *-T  GT  ( L . 6)  )*  0 .5 
TARGET  DIAGONAL 
XA=T^T( L , 1 ) -TGT  t L ,5 ) 

YA=TGT(L,2)-TGT(L,6) 

TO  = 0.5«((X4«X4  * Y\*YA)**(0.5) ) 

NCB'JHT  * 0 

DO  40  1=1,  NA 

«**  FOP  FACh  ATTACK 

IF  UR(I  ) .EO.  1)  GO  TT  40 
IN  I T 2 = 0 
NWPN  = ATT( 1,9) 

LTH  = -4M0(NWPN.1,1) 

IF  (LTH  .LF.  0)  GO  T~  40 
rtIO  = AMO (NWPN, 2,1) 

TOT  = Wr>(  NWPN ) + TO 
NS  = ATT ( 1,7) 

IN  IT  3 = 0 

00  20  M = 1,  NS 

•>■** ♦**»** *********  fgr  fach  weapon 

NCB'JHT  = NCRUHT  ♦ 1 
XB  = Cf?UHT(NC9UHT,  1) 

IF  (XB  .EO.  -10000.)  GO  TO  20 
Y8  = CBUHTI NCBUHT , 21 

DISTANCE  9FTWEEN  TARGET  AND  PATTERN  CENTERS 

0 = ( (XB-TCX)*(XB-TCX>  «•  (YB-TCY) *< YB-TCY ))**(. 5 ) 

TARGET  CANNOT  BE  HIT  IF  0 GREATER  THAN  COT 
IF  (D  .GT.  TOT)  GO  TO  20 
IF  ( INIT3  .GT.  0)  GO  TO  16 
1NIT3  = 1 

PHI  = ATT (1,1)/ 57. 3 
S = S IN ! PHI  ) 

C = COS (PHI  ) 

SL  = S*LTH 
SW  = S*WID 
CL  = C*LTH 
CW  ^ C*WIP 
00  5 J = 1,  16 


FGk  FACH  WEAPON 


GO  TO  20 


1 


^ >-,  jt*85-"iSsi ;-  4 V - , 


1281. 

00  5 K = 1,  16 

1282. 

5 

NC0V(J,K)  = 0 

1283. 

16 

CONTINUE 

1284. 

XX(i;  = X8  - ( $L<-Cto  1 /2  • 

1285. 

xx<2>  = xxm  + SL 

1286. 

X X ( 3 ) = 1X(2)  ♦ cw 

1287. 

xx { 4 1 = xxm  *•  cw 

1288. 

YY  C 1 1 = YP  ♦•  (SW-CU/2. 

128?. 

VY  ( 2 1 * YYU)  ♦ CL 

1290  t 

YY ( 3 ) = YY(2)  - SW 

1291 . 

YY ! 4 > = YY( 1 J - SW 

1292. 

IF  ( INIT2  .GT.  0)  GO  TC  18 

1293. 

INIT?  - 1 

1294. 

C 

FIND  WE STCf  LY  C JRNSR 

1295. 

ILX  = 1 

1296. 

DO  IS  = 1,4 

129  7. 

15 

IF  (XX(NN)  . L T . xxtuxn  ILX  = NN 

12' 8. 

IF  IS  .£0.  1.0)  ILX  = 4 

1299. 

r 

RENUMBER  CORNERS  So  TH1Y  CLRNEP  «1  IS  THE  MOST  WESTERN 

130C. 

ID I F = ILX  - l 

1301  . 

r 

v 

CDMPijT  c AND  ADJUST  r&N  AND  COT  AN  AS  RfOUIfttO 

1302. 

IF  US  .EO.  0.01  .r>P.  (C  .50.  0.0)1  GO  TO  28 

1303. 

IF  ((ILX  .LO.  2)  .CP.  (IIX  .80.  41)  GO  TO  26 

1304. 

T = S/C 

1305. 

GO  n ?7 

1306  . 

26 

T = -C/S 

1307. 

27 

CT  = l./T 

1308. 

28 

CONTINUE 

1309. 

IS 

DO  22  NN  - 1.4 

1310. 

NEW  * NN  - I 01 F 

1311. 

IF  (NEW  .LT.  1)  NEW  = 6 EW  ♦ 4 

1312. 

X(NcW)  = XX(NN) 

1313. 

Y(NtW)  = YY(NN) 

1314. 

22 

IF  (KTFST  .GT.  4)  WP IT E ( 6, 1003 (NSW, * < NEW  I , Y( NEW) 

1315. 

IF  ( INIT1  ,GT.  0)  GG  70  31 

1316. 

C 

CHEATS  A 16-PQI NT  GRID  ON  TARGET  - US  £ MOKE  POINTS  FOR 

1317. 

IN  IT  1 = 1 

1318. 

NX  0 = 8 

1319. 

NYO  = S 

1320. 

IF  ( T GT ( L « 1 2 ) .GT.  2S0.)  NYO  = 18 

1321. 

IF  ( TGT ( L » 12)  .GT.  1000.)  NYC  * 32 

1322. 

IF  I T GT  ( L « 1 3 1 .GT.  250.)  NXO  = 1(. 

1323. 

IF  1 TGT I L * 13 ) .GT.  1000.)  NXO  = 32 

1324. 

NXT  = NXO/2 

1325. 

NYT  = NY0/2 

1326. 

DO  29  J * It  NXT 

1327. 

DO  29  K = 1,  NYT 

1328. 

ICCVIJ.K)  - 0 

1329. 

29 

PSCOViJ.K)  - 1.0 

1330. 

NX1  = NXO  - 1 

1331. 

NY  1 = NYO  - 1 

1332. 

XO  = NXO 

1333. 

YO  = NYO 

1334. 

AO  = TGTIL.l) 

1335. 

A 1 = (TGT { L » 7 )-TGT  t L » 1 ) ) /XO 

1336. 

A 2 = ITGT<L,3)-TGT(L,1))/Y0 

1337. 

A3  = TGTIL.2) 

1338. 

A4  = (TGT(L,4)-TGT(L,2> J/YO 

1339. 

AS  = (TGT (L , 8) -TGT (L*2) )/X0 

1340. 

31 

CONTINUE 

1341. 

C 

TEST  TO  SEE  IF  TARGET  CORNERS  COVERED  6Y  PATTERN 

tfeiteo  J-^^&^’*mw^*i^>mxXsiaxtoW&iiSas^si!i.'Mt3tyi 


sSlt 


1383. 

1364. 

1385. 

1386. 

1387. 

1388. 

1389. 

1390. 

1391. 

1392. 

1393. 

1394. 

1395. 

1396. 

1397. 

1398. 
1395. 

1400. 

1401. 

1402. 


IF  ( KTEST  ...V.  4)  WRIT E ( 6, 1001 » XT  ,YT 
IF  <(v'  t '<(  1)  ).0R.  (XT  .GT.  X ( 3)  ) ) GO  TO  10 

IF  ('  .•  M4H.0R.IYT  .GT.  Y t 2 » J > GO  TO  10 

IF  • ■'  . _.)  .OP.  (C  .EQ.  0.) ) GO  TO  7 

IF  . . . |Y(1KTMXT-Xll))ll  G3  TO  10 

IF  (,.  ,-.(Y(2!-TMXT-X(2>m  GO  TO  10 

IF  (YT.LT.(Y(4)«-CT*(XT-X(4))))  GO  TO  10 
IF  (YT.LT.IYI l)-T*(XT-XCl))))  GO  TO  10 
MM  * N I N ♦ 1 
CONTINUE 

IF  (KTEST  .GT.  3)  WRITE  (6, 1002)  NIN 
IF  (NIN  . LT . 41  GO  TC  34 

IT  ALL  CORNERS  COVERED  BY  PATTERN,  T4RGET  FOLLY  COVERED 
00  32  J = 1,  NXT 
DO  32  X - 1 , NYT 
NCOV  < J , K ) = NCOVtJ.K!  + 1 
GO  TO  33 

IF  PARTIALLY  CnVFR?0,  ESTIPATc  FRACTION  THAT  IS  C.OVFREO 
CONTINUE 

OP  30  NX  = 1 ,NX1,2 
DC  30  NY  = 1 ,NY1,2 
J * (NXMI/2. 

K = (NYMJ/2. 

GRID-POINT  DIMENSIONS 
XT  = AC  N NX*  A1  ♦ NY  *A2 
YT  = A3  ♦ N Y* A4  F NX*A5 

CHECK  IF  WITHIN  R EC TANGLE  ENCLOSING  PATTERN  THAT  IS  PARALLEL  TC 


1342. 

1343. 

1344. 

1345. 

1346. 

1347. 

1348. 

1349. 

1350. 

1351. 

1352. 

1353. 

1354. 
1 355. 
1*36. 

1357. 

1358. 

1359. 

1360. 

1361. 

1362. 

1363. 

1364. 

1365. 

1366. 

1367. 

1368. 

1369. 

1370. 

1371. 

1372. 

1373. 

1374. 


NIN  = 0 
NGTN  = 0 

DO  10  NC  = 1, 
XT=TG T ( L »NC ' 
YT=TOT(L,Nf H ) 


1375. 

IF 

((XT 

.LT.  X(l)). OR. (XT  .GT. 

X ( 3 1 ) ) 

-4  r* 

DU 

TO 

30 

1376. 

IF 

( ( YT 

.LT.  Y ( 4 ) ) «CR. ( YT  .GT. 

Y ( 2 ) ) ) 

GO 

TO 

30 

1377. 

IF 

( (S 

. EQ.  1.)  .OR.  (C  .EQ.  0 

.))  GO 

TO 

35 

1378. 

C 

CHECK  IF 

POINT  IS  WITHIN  ACTUAL 

CBU  PATTEFN 

1379. 

IF 

(YT  . 

GT.(Y(1)*CTMXT-X(l)m 

GO  TO 

30 

1380. 

IF 

(YT  . 

GT.(Y(2I-  T* ( XT-X( 2) ) ) ) 

GO  Tn 

30 

1381. 

IF 

(YT  . 

LT.(Y(4H-CT*(XT-X(4))I  ) 

GO  TO 

30 

1382. 

IF 

(YT  . 

LT . t V ( 1 ) -T*( XT-X( 1 ) t > ) 

GO  TO 

30 

NGJN  * NG IN  ♦ 1 
NCOV( J , K } = NCOV ( J , K I 


IF  IKTfcST  .GT.  5)  WRITEI6.10 
CONTINUE 
CONTINUE 

IF  ((NIN  «•  NGIN)  .EO.  0)  GO  TO  20 
RECORD  ANY  COVERAGE  AS  A 'HIT' 

NHIT(L)  = NHIT(L)  + 1 
NCBUU)  = NCBU(L)  ♦ l 
IF  (KTEST  .GT.  4)  WRI TE (6,10 
IF  ( (TGT( L, 11 1 .LT.  1)  .PR.  ( 

IF  (NH1T(L)  .GT.  25 > GC  TO 
00  120  J = l,  NST 
IF  (MHIT(J)  ,E0.  0)  GO  TO  130 
IF  (MHIT(J)  . Nc.  L)  GO  TC  120 
HIT  (J, 1,NHJT (L)I=X8 
HIT  ( J, 2,NHIT(L))=YB 
HIT  ( J,3,NHIT(L))=NWPN 
IF  (NHIT(L)  .EQ.  25)  WPITE 
GO  TO  130 


WRITE (6, 1003)  NX, NY, XT  »YT ,NGIN»NCOV(J,K) 


WRITE(6, 1006)  L,  NHIT(L),  NCBU(L) 
l)  .OR.  (TGT(L,10)  .EO.  1.1  I GO  TO  130 
) GC  TO  130 


WPITE  (6,1007) 


ITRIAL 


1403. 

1404. 

1405. 
1406  • 

1407. 

1408. 
140  <9. 

1410. 

1411. 

1412. 

1413. 

1414. 

1415. 

1416. 

1417. 

1418. 

1419. 

1420. 
1421  . 

1422. 

1423. 

1424. 

1425. 

1426. 

1427. 

1428. 

1429. 

1430. 

1431. 

1432. 

1433. 

1434. 

1435. 

1436. 

1437. 

1438. 

1439. 

1440. 

1441. 

1442. 

1443. 

1444. 

1445. 

1446. 

1447. 

1448. 

1449. 

1450. 

1451. 

1452. 

1453. 

1454. 

1455. 

1456. 

1457. 

1458. 

1459. 

1460. 

1461. 

1462. 

1463. 
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120  CONTINUE 
130  CONTINUE 

IF  (TGTU.,10)  .N5.  1.1  G)  TO  150 
IF  (MHITtU  .GT.  250)  GO  TO  150 
00  140  J = 1,  5 
IF  (WW(J|  .50.  0)  GO  TO 

IF  (NRW! J ) .NE.  L>  GO  TO 

HITRC J,1»NMIT(l)>  = X3 
HITR! J,2,NhIT(L ))  = Y8 
HITR!J,3,NHlT(L>)  = NWPN 


150 

140 


RECYCLE  PQR  MORE  WEAPONS 


Rfc'CYCLS  FOR  MOKE  ATTACKS 


IF  (NHIT(L)  • cO.  250)  KRI TF(6, 1009)  L,  ITRIAL 
140  CONTINUE 
150  CONTINUE 

r i .**,•*** • ,**.*«* 

20  CONTINUE 

IF  ( INIT2  .EO . 0)  GO  TC  40 
PSP  = 1.  - AMOINWPN.TGTU  ,10), 2) 

DO  160  I = l,  NXT 

DO  160  F = 1,  NYT 

ICOVU.K)  = ICOV!J,K)  4-  NCOVI J »K) 

160  PSClV(J,K)  = °SCOV( J,K)*FSP**NCOV(J,K) 

40  CONTINUE 

IF  ( t lNtTl»l\‘lT2)  .GT.  0)  GO  TO  165 
P ( L * 2 ) = 0.0 
COV(L)  =0.0 
GO  TO  ?00 
TCOV  = 0.0 
PST  = 0.0 

00  170  J = 1,  NXT 

CO  170  K = 1,  NY" 

IF  UCOVIJ.K)  .GT.  0)  ICOV  = TCOV  * 1. 

PST  = PST  + PSC7V(J,K) 

7CTC  = NXT*NYT 
COV(l)  = TC0V/T7TC 
P(  L > 2 1 = 1.  - PST/T3TC 

IF  ( XTEST  .GT.  3)WRITg(6,  1004)L, TCOV. TOTC. PST, C0Va),P(L.2> 

• ,-*#******«****«*«**.*«*  RECYCLE  pCR  MORE  TARGETS 
CONTINUE 
RETURN 

F8.0, ’ YT  *,r8.0) 

• , 14) 

1 , 1 3, 2F10 • 0) 

• C0V*,F6.0,*  T0T',f=6.0,'  PST « 

COV ' * F 10.4, * P(L»2) ' , F 10 . 4 ) 

• NX  ’*13,*  NY  •, 13, 2F 10.0, 21 10) 

' TGT  * , 1 4 , ’ NHir’,15,'  NCBU',15) 

ONLY  25  HITS  WERE  STORED  FOP.  TARGET  #*, 
14,  • ***** ) 


165 


170 


C * 

200 


1001 

FORMAT!'  * 

,*  XT  *, 

10G2 

PJRMAT! * 

*,  * 

NIN 

1003 

FORMAT ! " 

*,  ' 

N5W 

1004 

FORMAT! • 

' , ' TGT ' , 14, 

1005 

1006 
1007 


F10.4,  * 

FORMAT!' 

FORMAT!  • 

FORMAT! *0* ,♦***»* 
13,*  DURING  TRIAL 


1009  FORMAT!  *0'  ,'***■**  ONLY  250  HITS  WERE  STORED  FOR  TARGET*, 

X 1 K » » 1 3, * DURING  TRIAL  #’,I4,»  **»«**) 

END 

SUBROUTINE  PLOT  H'f  ( NH  . NR.  L TH,  WI D) 

COMMON  /HITS/  XN( 250  I ,YNt  250 ) ,NZ( 250) 

DIMENSION  I COL  I 130) 

OATA  I8K  / 1H  /,IX/  1H*  /,IY/  1H*  /,1S/  1H-  /’,!£/  1H'  / 

C THIS  ROUTINE  PLOTS  THE  IMPACT  POINTS  (RUT  NOT  CRATERS)  FOR 

C ALL  HITS  THAT  HAVE  BEEN  STORED  FOR  A RUNWAY/ TAXI  WAY . IT 

C WILL  PLOT  ALL  HITS  THAT  AFFECT  RUNWAY  OPERATION  UP  TO  50  'FEET' 

C OF  EITHER  SIDE  OF  (UP  TO)  A 300  'FOOT*  RUNWAY.  RUNWAY  LENGTH 

C IS  LIMITED  TO  13000  * FFET  * . 


i 


7 ' ?'•'  « v»  » r- » 


If 

-7$} 

IJ 


1464. 

1465. 

1466. 

1467. 
1468  . 

1469. 

1470. 

1471. 
147?. 

1473. 

1474. 
147  5, 

1476. 

1477. 
1476. 
1479  . 

1480. 

1481. 
148?. 

1483. 

1484. 

1485. 

1486. 

1487. 

1488. 

1489. 

1490. 

1491. 

1492. 

1493. 

1494. 

1495. 

1496. 

1497. 

1498. 

1499. 
1500  , 

1501. 

1502. 

1503. 

1504. 

1505. 

1506. 

1507. 

1508. 

1509. 

1510. 

1511. 

1512. 

1513. 

1514. 

1515. 

1516. 

1517. 

1518. 

1519. 

1520. 

1521. 

1522. 

1523. 

1524. 


IV.  I D = 910/10.  * 5. 
UN  = L7H/100.  * 1, 
IF  (LC-N  .0T.  129) 
LU  = L5"/10. 

LI  = 10*t'l  6 1 


LSN  = 129 


LU  = LU 
00  40 


+ 1 

J = 1.40 
1 = 41 - J 
DC  10  N = 1,129 
ICOL(N)  = IBK 
ICCLIl)  = 15 
ICQL(LSN)  = IF 
IF  t(I  ,U£.  5)  .AND. 
PO  11  NS  - l.LFN 
ICOLtNS)  = IS 
00  12  NS  = i.LT, 
ICCL(NS)  = Ir 
CONTINUE 

00  20  N = l.NH 
NY  a YN(N)/10.  «■  5 


I W ID ) ) 


GO  TO  14 


TO  30 


GO'  TO  20 


IF  (HY  .NE.  I)  GC  TO  20 
NX  = *N  < N ) / 1 30 • 1. 

IF  ((NX  .LT.  1)  .04.  (NX  .GT.  129))  GO'  TO  20 
ICnL(NX)  = IX 

IF  <NZ(N)  .50.  1)  ICU(NX)  = IY 
0 CONTINUE 
Y = 1/5. 

LY  = Y 

IF  ( ( Y-LY)  .Nr.  3.0)  C TO  30 
LYY  = 5*LY  - 5 

WRITS  (6,1011  LYY.t  IC0L(V),M='l,U9) 

GO  TO  40 

3 rtRITC  (6,102)  (ICCL(M). *=1,129) 

0 CONTINU' 

WRITE  (5,103)  ( I,  1 = 1,  I?),  NR 

RETURN 

01  FORMAT  (•  * , I ? , 129 A l * 

02  FORMAT  (•  ' , 2X  , 1 29 A 1 ! 

3 FORMAT  (•  •,*  O',  1218X.12)/, 

X 40X , * TENi  BY  THOUSANDS  OF  LENGTH  UNITS  »,/,/, 

X 60X,  > 'M®aCT  POINTS  ON  RUNWAY  NUMBER  ',12 
X 37X , * U = POINT  IMPACT  WPN  * = C8U  CENTROID)') 

END 

SUBROUTINE  JMEMQ( NJ » PH , C ) 

COMMON/ I NT/ NT ,NA,ND  ,\’TM,'k  TEST, MCR  , MOW,  MOL .MODE, NPRINT ,NAH,NST , MTT 
DIMENSION  F ( 9 ) » 0(11) 

THIS  SU3R0UT INE  PRDVlOES  THfc  USER  'METHOD  O'  AS  OUTLINED  IN  THE 
8-74  WANG  700  USSRS'S  MANUAL  FOR  JM£m  OPEN 'END  METHODS. 

THIS  TRAJECTORY  PROGRAM  PERMITS  AIDA  ’HERS  TO  PRESCRIBE  THE 
ATTACK  DATA  AS  IN  JMJM.  THE  SUBROUTINE  LOGIC  IS  TAKEN  DIRECTLY  FROM 
THE  REFERENCED  PUBLICATION  AND  USES  NOTATION  CLOSELY  PARALLELING 
THE  ORIGINAL.  ONLY  THE  'PATTERN  RADIUS'  COMPUTATION  i USED  WITH 
ROCKEYc ) IS  OMITTED. 


;ENTROI! 


NJKEM  * NJMEM 
NCNT  = 0 
IF  I NJMEM  .60- 
READ  C5.102)  ( 

WRITS'  (6,103) 
WRITE  (6,104* 

DO  20  l = 7,9 

D(I)  = 0(1*1) 


n write  (6,101: 

( 5U>  , 1 = 1,91 
NA,  (7(1), 1=1,6) 
<£( t) , 1=1,9) 


•nSt&f&rfas  _«sa*as4S 


s&U?^  $ *"T?}  *;!-•»-’ 


1525. 
1626. 
1527. 
1528  • 
1629. 
1830. 

1531. 

1532. 

1533. 
1536. 

1535. 

1536. 

1537. 

1538. 

1539. 

1560. 

1561. 

1562. 

1563. 
1566. 

1565. 

1566. 

1567. 

1568. 
I56c  . 

1550. 

1551. 

1552. 

1553. 
1556. 

1555. 

1556. 

1557. 

1558. 

1559. 

1560. 

1561. 

1562. 

1563. 
1566. 

1565. 

1566. 

1567. 

1568. 

1569. 

1570. 

1571. 

1572. 

1573. 

1576. 
1575. 
1376. 

1577. 
I5?e. 
\579„ 
1580. 
1SC1. 

1582. 

1583. 
1586. 
1585. 


VI  = 

1.688*£f 

n 

V = 

VI 

TH  « 

t( 1) /57. 

3 

VX  * 

V*COSITM 

) 

VY  * 

-V*SIN(T 

H) 

Tft  * 

10(7)  - 

1. )*=<»)  * (e(5)-i.)*F(9> 

OX  » 

T»*VX/2. 

Y1  « 

t ( 2 ) - VY*  Ts 

YI  • 

(Yl  ♦ =.( 

2)  )/2. 

Y » 

V! 

TF  « 

O.T 

VF  * 

0.0 

THfso 

.0 

IF  ( j 

(7 i .GT. 

'■•00.) 

GO  TO  1 

Yf  * 

0.0 

IF  I: 

(I)  .co. 

0.T) 

GO  TO  ? 

TPD*r 

(7) 

GO  TO 

3 

YT  a 

FI  7 ) 

TPO* 

99. 

CONT! 

NUF 

OG  =■ 

32 . 1 7/  ( F 

(6  )*F<  ‘ 

>1) 

CONT ! 

NUC 

IF  (KTfcST  .G’ 

. 5) 

WRITE  (6,108)  TF,  y , v 

NCNT 

* NCNT  * 

1 

IF  ! NCNt  .GT. 

2000) 

GO  T7  30 

NCNT2 

* 0 

WRITFt6,107> 
GC  T7  30 


C = V/220.  - 3. 

TF  (C  .GE.  3.)  C = 3. 

IF  (C  .IT.  1. I C = 1. 

00  = C*DG*cXP(-Y/3lOOO.) 

OT  = 10. /<DG*V«V) 

IF  (OT  .GT.  0.5)  r>T  = 0.5 
GO  TO  13 
CONTINUE 
fF  = TF  3 

Y Y = V>0 

Y = YO 

IF  (KTEST  .GT.  5)  WRITFt6,107) 

NCNT2  = NCNT2  ♦ 1 

IF  (NCNT2  .GT.  100)  GC  TT  30 

CONTINUE 

TTO  = Tc 

TF  = 7 F ♦ OT 

IF  (TF  .IT.  tpo>  Go  TO  6 
TF  i TPO 
DT  « TF  ~ TFO 
CONT I NU? 

VYO  = VY 

VY  =1  VY*(  l. -V ♦OT^OD  I - 32.176DT 
YO  ' Y 

Y Y «■  OT*(VY*VY Ol/?. 

IF  ( ( Y - YT)  ,G=,  -1.)  GO  TC  6 
ll  ■=  7YG*V*00  «•  3T.17 
7 = vygwYO  « 2. 0*  { YO-YT  )*Ll 
OT  = WYO  «■  l**  1.5)  )Ul 
GO  TO  S 
CONTINUE 
VXO  = VX 

VX  = VX*U.-V*OT*OD) 

OX  = OX  t -'■'T*  (VX0+V41/2. 


^3*5? 


-'*m 

w 

I 
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•1 

•I 

'4 


1586  . 

V = (VX*VX  * VY«VY)«*{.5> 

1587. 

IF  (TF  .60.  TPO)  GC  TO  12 

1586. 

IF  ( F (7 ) .GT.  500.)  GO  TO  11 

1589. 

GO  TO  VO 

1590. 

11 

AID  = A3S ( ( Y-E( 7) ) ) 

1591. 

IF  (AID  .LT.  I.)  GC  tc  12 

1592. 

IF  (KTEST  .G£.  5)  WRITE  (6,109)  AID, 

1593. 

10 

CONTINUE 

1594. 

IF  (KTEST  . G? . 5)  Wfi I TF (6 , 108)  Y 

1595. 

IF  (Y  .GE.  1 . ) GO  to  7 

1596. 

GO  TO  15 

1597. 

12 

CONTINUE 

159  8. 

Z = -VY/V 

1599. 

ZZ  = ( l.-Z*Z )««(  .5) 

1600. 

THF=ATAN( 1/ LI ) 

1601. 

DG  = 32 . 1 7/  ( E ( 5 ) *E(  5)) 

1602. 

TPC=  99. 

1603. 

YT  = 0.0 

1604. 

IF  ( ? ( 7 ) .GE.  500.)  E ( 7 ! = 0.0 

1605. 

GO  TO  10 

1606  . 

15 

Z = -VY/V 

1607. 

22  - ( 1.-2*  Z)  **  t . 5) 

1608. 

AI  = AT  AN ( 2 / 2 Z ) 

1605. 

SR  = (OX*DX  ♦ Y1*Y1 ) '*(.5) 

1610. 

IF  (0(5)  .SQ.  0.0)  GO  TO  8 

1611. 

HP  = 0(41/1000. 

1612. 

CP  = 0(5)/1000. 

1613. 

GO  TO  9 

1614. 

8 

HP  = 0.000573*0(4) 

1615. 

CP  = BP 

1616. 

9 

CONTINUE 

1617. 

U = (fc(6CS(6CTF*VF) 

1618. 

1619. 

1620. 
1621. 
1622. 

1623. 

1624. 

1625. 

1626. 
162  T. 
1628. 

1629. 

1630. 

1631. 

1632. 

1633. 

1634. 

1635. 

1636. 
1631. 

1638. 

1639. 

1640. 

1641. 

1642. 

1643. 

1644. 

1645. 

1646. 


TF 


30 


101 


102 

103 

104 

105 

106 
101 
108 

109 

110 


0(4)*  (U+(SR*S  R*8P/Y  11**2.  >**(.  5) 
0(5)=(U*CP*CP*SR*SR  ) **(.5) 

SRT  = SRI 1000. 

0(8)  = V1*TR*SIN( AI-TH)/SIM(Al  ) 
OISP  = 1(6) 

P(  10 ) = 0.674«<M  SP*SRT 
0(6)  = 0(101/2 
AID  = 57 . 3*A! 

THFD  = 57. 3* TUP 

WRITE (6» 105)  To,  SR,  AID,  THFD 
RETURN 

WRITE  (f.HOI 
STOP 

FORMAT Cl'»20X,'JMEM 


FORMAT  ATTACK  DATA*//'  * 


X • ATT 

NO. 

HOG 

X-MPI  Y-M° I 

CEP(REP) 

1DEP) 

DIS 

XP'/' 

X' 

' , 

SPEED 

ALTITUDE  DIVE 

TERM! 

TERM2 

WIN 

XD 

TD/HF 

T TD',/,/) 

FORMAT ( 6X,  6F6.0,  3F6.3) 

FORMAT  ( 1 O’ , I 6, 7F10 .0 1 
FORMAT  (•  *,6X,  6F10.0,  3F10.3) 

FORMAT ( * • , 30X, * TF  »,F5.2,'  SEC  SR 
'IMPACT  ANGLE  '.F5.2,'  DEG  (FUZING 
FORMATC  **  • , F6.3. * SEC 

FORMAT ( • ',10X,'  OT  • ,F10.4, ' 

FORMATC  ' , • A LT  ',F10„2> 

FORMATC  • , ' A LT  OIFF ' , F 10.2, * TF  '.F10.2) 
FORMATCO'/'  ', 'LOOPING  IN  JMEMO.  CHECK  INPUTS 
• WITH  KTEST  = 6.0* ) 


’ »F7 .0, ' FEET  ', 

ANGLE  ' , F5  .2 , * DEG)') 

,F8.0, ' FEET  ' , F8. 1 ,'  FT/SEC  ) 
VERT  VEL ' , F 10 . 3 ) 


ANO/OR 


f 

I 

J 

I 

'I 

*"U 

$ 

i 

$ 

1 


84- 


1*47. 

END 

1648. 

SUBROUTINE  EXPHIT 

1640. 

COMMON  /ARPAY?/  TGT ( 250,13),  AtT( 

1650. 

XI  ZONE (50,2), NH IT (250) ,MHIT <20 > ,HIT 

1651. 

X ,P! 250, 3) ,COV ( 250) , M TYPE ( 10), NAME 

1652. 

COMMON/ 1 NT/ NT ,N  A , NO , NTM, K TEST  »MC« 

165  j. 

COMMON  / HI  TON  ( NH I TO , NR~DO 

1654. 

C 

THIS  ROUTINE  ESTIMATES  THE  EXPECT 

1655. 

C 

POINT-IMPACT  WEAPONS  ONLY. 

1656. 

c 

1657. 

c 

1653. 

c 

1650. 

NCYCLE  = C 

1660. 

PI  - 3.14159 

1661. 

00  200  I = l,N6 

1662. 

NW  = ATT ( 1,9) 

166?. 

IF  (AM0(NW,1,1»  .LT.  0.0)  GO  TO 

1664. 

REL  = WPNREL(NW) 

166  5. 

NFL AG  = 0 

1666  . 

NS  = ATT  ( 1 , 7) 

1667. 

LS  = ATT ( 1,8) 

1666  . 

PH 1 = ATT! 1 , 1 ) / 5 ? . 3 

1660. 

S = SlN(PHl) 

1670. 

C = CCS(PHF) 

1671. 

IF  ( ( A?  T ( I , l ) . S3 . 0.)  ,{R.  (ATT! 

1672. 

CT  = C/S 

1673. 

GO  TO  20 

1674. 

10 

NFL AC  = 1 

1675. 

20 

CONTINUE 

1676  . 

NP  = 10 

1677. 

OLS  = L S/9. 

1676. 

IF  (LS  .GT.  500)  G'  TO  30 

1670. 

MP  = 5 

1680. 

OLS  = LS/4. 

1681. 

IF  (LS  .GT.  50)  GO  TC  30 

1682. 

NP  s 2 

1683. 

OLS  = LS 

1684. 

30 

CONTINUE 

1685. 

X = AT  T { 1 , 2 ) - 3*1  S/2  . 

1686. 

Y = ATT ( 1,3)  - C*LS/2. 

1687. 

SIGRS  = 2.200*ATT( I » 4) *ATT( I ,4) 

1688. 

SIGOS  = 2 .200  'A  TT(  1 ,5)*ATT(  1 ,5)  ♦ 

1680. 

TSRS  = 2.  *SI GRS 

1690. 

TSCS  = 2.  * SI GOS 

1691. 

SIGR  = S I G°  S* ‘I .5) 

1692. 

SIGO  = SI  GOS** ( . 5 ) 

16)93. 

F = 1.  / (SIGR*SIGO*6.2832) 

1694. 

IF  (KTEST  .GT.  1)  WRITE (6, 1008 ) 

1695. 

X SIGO,  c 

1696. 

NFL  AG 2 = 0 

169  7. 

DO  180  L = 1 , NT 

1698. 

IF  (TGT (L ,10)  .SO.  21.)  GO  TO  1 

) 699. 

NFLAG2  = 1 

1700. 

DEN  = 0.0 

1701. 

TL1  = TGT (L , 1 2 ) 

1702. 

TL2  = TGT (L, 13) 

1703. 

TH  = TGT  ( L,  91  /5  7.3 

1704. 

CHI  = TH  - PHI 

1705. 

SC  = SIN(CHI) 

1706. 

CC  = COS (CHI) 

1707. 

SI  = ($IGRS*CC*CC  ♦ SI GOS*SC*SC ) * 

N'JMPER  5 OF  HITS  FOR  ATTACKS  WITH 


240 


.cO  • 1 BO. ) ) no  TO  10 


ATT ( 1 i 6 )*ATT( 1,6) 
ATT ( lt 1 0 )*ATT l 1,10) 
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1708. 

1709. 

1710. 

1711. 

1712. 

1713. 

1714. 

1715. 

1716. 

1717. 

1718. 
1715. 

1720. 

1721. 

1722. 

1723. 

1724. 

1725. 
1126. 
1727. 
1728  . 

1729. 

1730. 

1731. 

1732. 

1733. 

1734. 

1735. 

1736. 

1737. 

1738. 

1739. 

1740. 

1741. 

1742. 

1743. 

1744. 

1745. 

1746. 

1747. 

1748. 

1749. 

1750. 

1751. 

1752. 

1753. 

1754. 
17';  5. 

1756. 

1757. 

1758. 

1759. 

1760. 

1761. 

1762. 

1763. 
176*1. 

1765. 

1766. 

1767. 

1768. 


S2  = (S1GRS*SC*$C  + SlG0$*CC*CC)*M .5) 

Nl  = C 4 . *■  TL  1 / S 1 1 ♦ 1 
N2  = (4.  * TL  2/  S2  ! 1 

IF  IKTtST  .GT.  1)  WRITF (6, 1009)  Nl,  N2,  SI.  S2 
00  90  K = 1.7,2 
XT  = TGT(L.K) 

YT  = TGTU.K4-1) 

XX  - X - XT 

IF  (NFLAG  ,E0.  1)  GO  TO  40 
YY  = YT  - Y «-  CT*XX 
0 = -YY*  S 
R = YY*C  - XX/S 
GO  T')  50 


40 

0 = 

XX 

R = 

(YT  - Y ) * C 

60 

AI01  * 0*0/TS0S 

IF 

( ! TL If TL2 ) .GT. 

200.) 

GO  TO  60 

IF 

(Aim  .GT.  10.) 

GO 

TC 

180 

40 

IF 

(A  101  .GT.  10.) 

GO 

TC 

90 

00 

= SEXP(-AIOl) 

00 

80  M = l.NP 

AID2  = R*  R/TS  RS 

IF 

t A 102  .GT.  10.) 

GO 

TO 

70 

DR 

= SFXP (-AI 0? ) 

OEM  = OEM 

NCYCLE  = NCYCLfc  ♦ 1 

IF  (KTEST  .GT.  2)  wRlTft6»1005»>„L»K,XT,YT,0,P,DC,QR,CfcN 
70  R - R - OLS 
80  CONTINUE 
90  CONTINUE 

C IF  THE  TARGET  DIMENSIONS  Ac  E SMALL  (I.E.  LESS  THAN  CNE-CUARTER 
C THE  PROJECTION  7F  SIGMA  PARALLEL  .TO  THE  TARGET  EDGE)  THE  HIT 
C DENSITY  IS  TAKEN  AS  THE  AVERAGE  C1-’  THE  VALUES  AT  THE  FOUR  CORNERS. 

C IF  IT  IS  LARGER,  A GR 10  CF  INTERNAL  POINTS  IS  ESTABLISHED  AND 

C THE  HIT  DENSITY  IS  TAKEN  AS  THE  AVERAGE  0VC»  THE  CORNERS  AND  THE 
C INTERNAL  POINTS. 

IF  ((Nl  + N2)  .GT.  2!  GO  TO  100 
OEN  = DEN/ ( 4. *N° ) 

GO  TO  160 
100  CONTINUE 

STC  = SIN(TH) 

C7G  = COS(TH) 

DILI  = TL1/IN1+1) 

0IL2  = T(_ 2/  ( N2+1 ) 

DO  150  M = 1,N1 
DO  150  N = 1,N2 

XT  = TGTIL.U  M*niu*STG  ♦ N*0IL2*CTG 
YT  = TGTIL.2)  ♦ M*0IL1*CTG  - N*0IL2*STG 
XX  = X - XT 

IF  (NFLAG  .FO.  1)  GO  TC  110 
YY  = YT  - Y ♦ CT*XX 
D = -YY*S 
R = YY*C  - XX/S 
GO  TO  ) 20 
110  0 - XX 

k = (YT  - Y)*C 
120  A101  = D*D/TSOS 

IF  (AID1  .GT.,  10.)  GC  TO  150 
OD  = SEXPI-AID1) 

DO  140  K = l.NP 
AID2  * R*P/TSRS 


1769. 

mo. 

1771. 
1772  . 
1773. 
177A. 
1775. 
1 776  . 
1 777. 

1778. 

1779. 
178C . 

1781 . 

1782. 

1783. 
170A. 
1785. 
1786  • 

1787. 

1788. 

1789. 

1790. 

1791. 

1792. 

1793. 
179A. 

1795. 

1796. 

1797. 

1798. 

1799. 

1800. 
1801. 
1802. 
1803. 
180A  . 

1805. 

1806. 

1807. 

1808. 

1809. 

1810. 
1811. 
1812. 
1813. 

181  A. 
1215. 
1816. 

1817. 

1818. 

1819. 

1820. 
1821. 

182  2. 
1823. 
182A. 

1825. 

1826. 

1827. 

1828. 
1829. 


160 


IF  < A I D 2 .GT.  10.) 

OR  = SFXP (-41 02 ) 

OEN  = PFN  OR*OC 
NCYCLE  - NCYCIF  + 1 

IF  (kTCST  .GT.  2)  WPITEI6, 1006)1,1,  XT,  YT  , 0 R,DD,CR,DFN 

R = a - OLS 

CONTINUE 

CONTINUE 

OcKI  = DEN/  IIA»N1*N2)-»NP) 

CONTINUE 
0EN1  = DON 

DEN  = NS*F*RCl*ATTU  ,li)*OSNl 

IF  (KTrST  .GT.  1)  WRITE  (6,1007)  I , L »DEN1 ,DeN 
SMC  = 4M0(NW,TGT{L, 10) ,1) 

OLl  = TL1  ♦ R*L  Mr- 
012  = f L2  ♦ 2*FM0 

EH lT  1 = OEf!*(  DL  l^OL  2 - . 8 5rtA*f MO«EMO) 

COVtLI  = COV(L)  t-HITl 
CONTINUE 

IF  (NFLAG2  .£0.  0)  GO  TO  230 
CONTINUE 

WRITS <6 ,1000)  MCYCLF 
WRITE  (6,1002) 


210 


220 

230 

2A0 

1000 

1001 

1002 

1003 

100A 

X 

1005 

X 

1006 

1007 

> 

1008 

1009 

1010 


00  220 
NN  = 0 
TCOV  = 
9LCGHT 
00  210 


M = 1 , MTT 


MJ 


GC  TO  210 


0.0 
= 0.0 
L = 1, 

IF  (TGT(L.IO)  . 

NM  = NN  ♦ 1 

TCCV  = TCOV  + CDV(L) 

FLDGHT  = PIOGHT  «•  I 1 . - c XP ( -COVI L ) ) ) 

IF  (NN  .SO.  1)  wFITE  * 6 , 100 1 ) M 

WRITS  16,1003)  UCOVIL), \AMt(l,l) ,NAME!L,2) 

CONTINUE 

IF  ( NN  .SO.  0)  GO  to  220 

8L0GHT  = BIDGHT/NN 

WRITE  (6,1010)  TCOV,  PLCGHT 

CON  f I NU2 

CONTINUE 

RETURN 

WP.ITF  <6,  1006) 

STOP 

FORMAT ( ' 1'  * 25X,  'CYCLES', 

FORMAT  C0',10X,'«*  TARGET 
FORMAT  { 'O' , 10X, 'TARGET 


IV//) 

TYPE  13,'  *«',/) 

HITS', /,10X,»  NO  EXPECTED') 

FORMAT  I • ' ,10X,IA,6X,F6.3,6X,2AA) 

FORMAT  Cl','  COMPUTATION  STOPPED  - CBU  WEAPONS  ARE  NOT', 

• PERMITTED  WITH  SU8RQUTINE  EXPHIT') 

FORMAT!'  *,'C0  ATT  ',13,'  TOT  ',13,'  COP  ',12, 

AF6.0,  3E12.A  ) 


, "IP 


FORMAT  I 
FORMAT! • 

F8.5, • 

FORMAT  ( 

FORMA!  < • 

FORMAT! • 

END 

FUNCTION  SEXPIX) 

IF  (X  .IT.  -0.025» 
SEX°  = l.*X 


• ,I3,AF6.0,3E12.5) 
,13,'  AVG  OEN  ', 


ATT  ' , I3» • T GT 
' ,'  ATT  ',13,'  TGT 
NOR  DEN  * »F 1 A. 10, // ) 

' , ' ATT  ',13,218,  F8.2,  AF8.1,  612.5) 

' , ' N1  * , I A, * N2  ' , IA,6X,2F10.l> 

• ,20X  , ' *,/,19X,F8.3,M',F5.3,')» ) 


GO  TO  10 


I 


Wm 


1830. 

1831 . 

1832. 

1833. 
1834  . 
1835. 
1836  • 
1837. 
1638. 
1835. 

1840. 

1841. 
1842'. 

1843. 

1844. 

1845. 

1846. 

1847. 
1843. 

1849. 

1850. 

1851. 

1852. 

1853. 

1854. 

1855. 

1856. 

1857. 
1856. 

1859. 

1860. 
1861. 
1862. 

1863. 

1864. 
1065. 
1866. 
18'7. 
1868. 

1869. 

1870. 

1871. 

1872. 

1873. 

1874. 

1875. 

1876. 

1877. 
1876. 
1879. 
1380. 
1881. 
1882. 

1883. 

1884. 

1885. 

1886. 

1887. 

1888. 

1889. 

1890. 


RETU°N 

SEXP  = =XP(X) 

RE  Tim 
END 

SIBPOUTINF  GRI05N 

C9MM.QN/ARF AY S/TGT( 250,13) , ATT (50,11)  »AM0U0 ,20  ,2  > ,T0<  250 , 2»  , 

X 1 ZONE ( 50,2  ),NH IT (250) ,MKI  T ( 20 ) ,HJ  T(  2 0, 3,25 ) ,NR W(  5 ) ,H I TR(  5 , 3 , 250) 
X.P ( 250,3 > ,Cr VI 25  0) , MTYP; { iO) , NAME ( 250,2  I, WPNAEU  10 ) , NCBU ( 250) 
COMMON /I  NT /MT,NA,f!0,NrM,K  TEST,  MCA,  mCw.MCL  , MODE ,NPP INT , NAM, NST , MTT 
01 MENSION  OAT A( 17 ) 

NCYCL5  = 0 
PI  = 3.14154 
SR3 1 = 1.7724 
DO  140  L = i,N7 

IP  (TGT(UIO)  .NP.  21.)  00  TO  140 

N P L AG  2 * 0 
INC  = 250 

Ir  (TOT (L , 9 ) .OT.  0.)  INC  = TGT(L,9)/16. 

NX1  = TGTIL.l) 

NYl  = TGT  < I,  , 2 ) 


10 

20 


30 


DO  120 
DO  120 
NY  = 
DEN 
XT  = 
YT  = 


00 


NYR  = 1,17 
NX  - 1,17 
18  - NYR 
= 0 

NX  1 «■  (NX-1)* 
NY  1 + (NY-11* 
C4TA(NX)  * 0.0 
100  1=1,  MA 

NW  = ATT (I ,9) 

IF  ( AMC( NX  , 1 , 1 ) . 
RSL  = WPN* cL ( NW I 
NE  LAO  = 0 
NS  = ATT (I ,7) 

US  = ATT( i ,8) 

PHI  = ATT ( I » 1 ) 757 , 
S = SI N(DH! ) 

C = COS (PH! ) 

Ic  ((ATTU.l) 


INC 

INC 


LT.  0.) 


GO  T'T  160 


.60.  0.)  .OR.  ( ATT ( I , 1 ) .PQ.  130.))  00  TO  10 


5001  GC  TO  30 


50) 


GO  TO  30 


CT  = C/S 
GO  TO  20 
NCL AG  = 1 
CONTINUE 
NP  = 10 
OLS  = L S/9. 

IF  (IS  .GT. 

NP  = 5 
OLS  = l S/4. 

IF  (IS  .GT. 

NP  = 2 
OLS  = LS 
CONTINUE 

X = ATT (1,2)  -S*L  S/2. 

Y = ATT  (1,3)  -C*I.S/2« 

SIGN'S  a 2.200*ATT(I,4)*ATT(I,4) 
SIGCS  = 2.200*A7T( I,5)*ATT( 1,5) 
T55  S = 2.0  * SIGRS 
TSOS  = 2.0  * SIGDS 
Si  CR  S*  * ( .5  I 
SIG0S**(.5) 

/ (SIGR*SiG 0*6, 2832) 

- XT 


♦ ATT(I,6)*ATT(I,6) 

♦ ATT ( I ,10)*ATT (I , 10) 


StGR  * 
S1GD  * 
F = 1. 
XX  = X 


w 


hi 


m 


If 


-88- 


1891. 

IF  (MFLAG  .60.  11  CO  TO  *0 

1892. 

YY  * YT  - Y ♦ CT*XX 

1R93. 

- 

D * -YY«  S 

189*. 

o * YY*C  - XX/5 

1895. 

GO  T'i  50 

1896. 

*0 

n * XX 

1897. 

l * (YT  - Y)«C 

1898. 

50 

At.'l  * C*0/TSDS 

1899. 

tF  tAlDl  .GT.  10.1  Gl  TO  90 

1900. 

rn  t SfXFI-AlOU 

1901. 

I.CYCLe  ■ NCYCle  ♦ 1 

190?. 

00  80  9 * 1*  N« 

1903. 

A102  * R»R7TSRS 

190*. 

IF  {*102  .GT.  10.1  G(  T"  70 

1909. 

°F.  « 5r  X°  { -*T  32 1 

1906. 

Os‘1  * l-tv  ♦ 0R*Pr 

1907. 

NCYCL.-  * HL  YC  Li  ♦ l 

1906. 

IF  (KT!‘T  .GT.  21  VIP  1 Tc  (1005)  XT , YT  , 0,  R ,0D  ,CR,PEN 

190c  . 

70 

8 * « - OLi 

1910. 

80 

CONTI  0= 

1911. 

90 

cont  nu  ; 

191?. 

DC I,  * IOOOU»NS»F*0=L*ATT( 1 ,11>«DEN/NP 

1913. 

r*T*i‘ixi  = ri/Ti(‘:xi  ♦ Ti.N 

191*. 

100 

CON T I Mi - 

1915. 

IF  (NFLAG2  .CO.  1)  Gl  TO  110 

1916. 

MrLtO>2  * l 

1917. 

WRITE  (6,  1001) 

1518. 

no 

CONTI  NIJB 

1919. 

IF  (NX  .IT.  17)  GO  TF  120 

1920. 

NYT  = YT 

1921. 

IF  (NY  .1  T.  17)  50  TC  115 

1922. 

NX2  = NX1  ♦ 16»1NC 

1923. 

WRITS  (<.,1002)  ( MCC,  NLPC«NX1,NX2,INC) 

192*. 

115 

WOlTfc  (6,10031  NYT,  »i>ATA(l),  1*1, IT) 

1925. 

120 

CONTINUE 

1926. 

1*0 

CONTINUE 

1927. 

MS' IT:  (6,100*1  NCYCIE 

1928. 

r(ETU?N 

1929. 

160 

CONTINUE 

1530. 

WHITE (5,10061 

1931. 

STOP 

1932. 

1001 

F0RMAT(*l*,20X,*SXi>ECTH>  HIT  DENSITY  PER  10000  SO  FT 

1933. 

1002 

FORMAT!*  * ,// ,30X, • X - LOCATION*  »/, * *,lOXt17I7,/, 

193*. 

X 

• Y L3C  • ) 

1935. 

1003 

FORMAT  ( * *,/,*  *,  16,  * X,  17F7.3  1 

1936. 

100* 

FORMAT l • *,///, 20X, ‘CYCLES  IN  GRIOEN  ',181 

1937. 

1006 

FORMAT! • *,///,*  COMPUTATION  STOPPED  - C8U  WEAFONS 

1938. 

X 

• NOT  PERMITTFO  WITH  SUBROUTINE  GRIDEN* ) 

1939. 

1005 

FORMAT I ' * ,*F6.0,3E12.* ) 

19*0. 

ENC 

